Formas diferenciais e aplicações em 
eletromagnetismo e termodinâmica 


Valmir Peixôto de Souza Sobrinho 


Setembro 2022 


Sumário 


1 


2 


Revisão de cálculo vetorial 

1.1 Teorema fundamental do cálculo ...lccccccccl 

1.2 Cálculo em várias dimensões ..lcccccccccclc 
1.2.1 Produto escalar e vetorial ....ccccccccccc 
1.2.2 Operador nabla . ...ccccclc 
1.2.3 Integrais de linha... ccccccll 
1.2.4 Teoremas fundamentais ...c.cccccccccccc 
1.2.5 Integração por partes .....ccccccclcc 

1:35" "Jacóbiano” sm sim es ds SQL Ea E EA SE SG GN ÚE dado dera Lied do DD 
1.3.1 Áreas e volumes infinitesimais sob parametrização 
1.3.2 Áreas e volumes orientados ...ccccccciiicii 

1.4 Vetores contravariantes e covariantes . ....ccccclccc. 
1.4.1 Transformações contravariantes . ....cccccccc. 
1.4.2 Transformações covariantes .....ccccccclcc. 
1.4.3 Vetores sob mudança de base .....c.ccccclccc. 

1.5: BBcudovebOres:o mz ras EE dia E E a CA GUESS STA, GDE aba o E 


Espaço dos p-vetores 

2.1 Determinantes .....ccccccccccccc 

2.2 Produto externo entre p-vetores . . . cc ccccccls 

2.3 BrOduto: LNLELDO sta e drag dad ad BOA UBE cf À PA AA AS 
2.3.1 Produto entre I-vetores .....cccccccccs 
2.3.2 Produto entre p-vetores .....ccccccclc 

2.4 Operador estrela .....ccccccccclc 
2.4.1 Calculando o operador estrela. . cc clic. 

2.5 Transformações lineares ....ciccccccccclc 


3 Derivada exterior 
Sd. -Espaçodual ss quan gs dedo Es Al E SD RS a E 
32, - Espaço COLangente. aguia tita ed Sa oi e ca aa as Dali 
3.3 Álgebra de p-lormas= a atas Sao dd ÇA jaca UR NR cf ga pa 
3.4 Mapeamento e mudança de coordenadas . .....cccccc. 
3.5 Inversa do lema de Poincaré . ....cccccccllc 
3.6 Coordenadas generalizadas ....lcccccccccc 
3.7 Gradiente, Divergente, Rotacional e Laplaciano . .. cc... 
38 SUperfícies . 2 ess ur mo E eo ce a a E EN JE ES ia 


4 Formas diferenciais e integração 

Hd. Variedades sn sao ar tear e DGE BRO RR AR SE TAP e a a a 

4.1.1 Diferenciabilidade ...ccccccccl 

4.1.2 Orientabilidade ....lccclclll 

4.1.3 Vetores tangentes. .... ciclo 
4.2 Formas diferenciais . . ..ccccccccl 
AS. “SIM pIeXOS E. jar zaas e e ga 2 Coe ça DO E JR O a Rg o RÃ ar oa 
4.4 Integração de formas diferenciais ....ccccccclccccc. 
4.5 Teorema de Stokes generalizado . . ....lccccclcc 
4.6 Formas exatas e fechadas ....ccccccclccc 


5 Aplicações de formas diferenciais na Física 
5.1 Eletromagnetismo .....ccccccccccc a 
9.1.1 Equações no vácuo ....ccccclcl 
5.1.2 Teorema de Poynting. ....ccccccc 
5.2 Tetmodinâmica-sc. Gg mis nas E da Ee MIS Gon GUN) Car a E o E 
5.2.1 Relações de Maxwell . . Lc cccclcccc 
5.2.2 Propriedades dos materiais ....cccccccccc 


1 Revisão de cálculo vetorial 


1.1 Teorema fundamental do cálculo 


Teorema 1.1 (Primeiro teorema fundamental do cálculo). Seja f uma função 
integrável em um intervalo [a,b|, c um escalar tal quea < c< be F uma função 
definida pela expressão 


Flo) = [ rod Yxa<xr<b» 


Então uai F'(x) existe para todo « em um subintervalo I de (a,b) onde f é 
x 


contínua e 


F'(x)= f(x) Vrxel 


Demonstração. Pela definição de derivada de uma função, temos 


onde 


a+h a 
Flo +) = Pta) = | ra [ Ft) dt 


E [mo dt + (ro dt 
h 


= fo f(t) dé 


Considerando x um valor fixo, podemos fazer f(t) = f(x) +[f(t) — f(x)] e então 


a+h 
F(z+ h)- F(s) = Pi fa) + [160 — Ha] dt 


de forma que 


F(x + h)— F( 
h 


a+h 
Do fno+;/ 10-1od 


Como f é contínua em um intervalo T, tem-se, por definição, que Ve > 0 existe 
um ô > 0 tal que 


a-ô<t<r+ro = I|f()-Ffa<e (1) 


para todo x,t em 1. Agora escolhendo um h tal que rx <xvx+Hh<ax+o, 
ou equivalentemente O < h < 6, a proposição (1) será válida no intervalo de 
integração [x, «x + h], no que implica 


a+h a+h ah 
Há F(t) — f(x) dt <[ 0 = fHojjat< [ Ediede 


ou seja, 
E < 
— € 
h 


ac+h 
É H(0) — f(x) dt 


Mas essa é a definição de limite de h e é equivalente a dizer que 


a k A) Fledi=0 


h>0+ h 


Analogamente, para esse mesmo ô e e, podemos escolher h < O tal que O < 
—h < 6, ou equivalentemente O < |h| < ó para chegarmos à conclusão que 


Como ambos os limites laterais são iguais, temos que 


li E 
mm — 
h50 À 


ag+h 
[ to-tma-o 


e como consequência 


dE a Plz+h) — Pla) 
de ho0 h 
a+h 
«poros; f 016947 
= f(x) 
para todo x em T. E 


Teorema 1.2 (Teorema de Rolle). Seja f uma função contínua em um intervalo 
[a,b] e diferenciável em (a,b) tal que f(a) = f(b). Então existe ce (a,b) tal 
que f'(c) = 0. 


Demonstração. Suponha que não exista c e (a,b) tal que f(c) = 0, isto é, a 
função f não tem mínimos ou máximos em (a,b). Isso significa que f tem seu 
máximo e mínimo em [a,b] justamente nas extremidades. Mas como f(a) = 
f(b), segue que f é constante em [a,b], no que implica f(x) =0Vxe (a,b), o 
que contradiz a hipótese. E 


Teorema 1.3 (Valor médio). Se f é uma função contínua em [a,b| e f” existe 
em (a,b), então existe algum ce (a,b) tal que 


Demonstração. Seja g uma função definida por 


g(x) = bf(x) — af(a) — 2 f(b) + 2 fla) = Ha)(b — a) — alf(b) — F(a)] 


Então g(a) = bf(a) — af(b) = g(b). Como f é contínua em [a, b] e diferenciável 
em (a,b), então multiplicar f por uma constante e somar com uma função afim 
alf(b) — f(a)] manterá a continuidade e a diferenciabilidade de g em [a,b] e 
(a,b) respectivamente!. Aplicando o teorema de Rolle (1.2), temos que existe 
ce (a,b) tal que g'(c) = 0. Mas também 


Como f é integrável em [a, b|, segue que E 


dE 
Teorema 1.4 (Derivada nula). Se E O para todo vz em um intervalo aberto 
bi 


I=(a,b), então F é constante em T. 


' dr , o Pe 
Demonstração. Como mio F'(x) é uma função limitada em (a, b), segue que 
JJ 
F é uma função contínua. Então podemos usar o teorema do valor médio (1.3) 
para mostrar que para qualquer x em [a,b] existe algum ce (a,b) tal que 


F(x) — F(a) = F(c)(z — a) 


Mas F'(c) = 0, no que implica F(x) = F(a) para todo x em [a,b], isto é F é 
constante em [a, b]. E 


Definição 1.1 (Funções primitivas). Se F é uma função tal que 


dF 
E F(x)= f(x) Vrxel 

para alguma função f, então denominamos F uma função primitiva de f no 
intervalo T. 


Teorema 1.5 (Diferença entre funções primitivas). Se F e G são primitivas de 
f em um intervalo I, então F e G só diferem de uma constante em 1. 


“Note que | f(x) + ax — f(t) — at] < | f(x) — f(tb)l+ lal-lz — t| < ese | f(x) — f(b)l < /2 e 


lz—t|<ô paratodo xr tal que |x—t|<ô6= 2a? o que garante a continuidade de g. 
a 


Demonstração. Pela definição de primitiva, temos 


Fa) — G'(x) = Ha) — f(x) =0 


Usando a propriedade distributiva da derivada, 


Pelo teorema da derivada nula (1.4), temos que 
F(a) —- G(ax)= O 
onde C é uma constante. E 


Teorema 1.6 (Segundo teorema fundamental do cálculo). Seja f for uma 
função contínua em (a,b), ce (a,b) e F qualquer primitiva de f em (a,b). 
Então para todo x em T 


Fla) = (9 + [ sa 


Demonstração. Seja A uma função definida pela expressão 


A(x) = a f(ty)dt Vxe(a,b) 


Então pelo primeiro teorema fundamental do cálculo (1.1) temos que A'(x) = 
f(x) para todo x e (a,b) Além disso, A é uma primitiva de f. Como F também 
é uma primitiva de f, o teorema (1.5) nos garante que 

F(a) —- A(o)=€ 
onde C é uma constante. Mas quando «x = c, temos A(x) = 0, no que implica 


F(co)=C 


Como consequência, temos 


ou equivalentemente 


Fla) = (+ [ sa 


O teorema (1.6) também aparece mais comumente na forma 


[ soar=r(-ra (2) 


Teorema 1.7 (Lema fundamental do cálculo variacional). Se f é uma função 
contínua em [a,b] e para qualquer função contínua n(x) que satisfaz as condições 
de contorno n(a) = n(b) = O for válido 


b 
/ Ha)n(x) de = 0 


então f(x) = 0 Yxe [a,b]. 


Demonstração. Suponha que f seja positiva em um intervalo (11,12) € [a,b] 
com q * «92. Seja n(x) a função definida por 


(x — xy)(xa — 2º) q € [71,45] 
me) = | 0 É x É [01,45] 


Se x e (x1,x2), segue que 14 < x < x e consequentemente x —- x > 0€ 
g9— 4 > 0,0 que implica n(x) > 0 Vx e (a,b). Então 


f siemeras= [P rtomtojar+ [ rtonteras + [stone do 
- [. Hentoj as 


que é positivo, uma vez que f(x)n(x) é positivo em (71,12). Isso contradiz a 
hipótese da integral ser zero e portanto f(x) = 0 Vx e [a,b]. E 


1.2 Cálculo em várias dimensões 


Para facilitar a notação, iremos definir 


d=>— cm q=e0=% qp=e=y qg=e3=2Z 
Oqi 


e também 


1.2.1 Produto escalar e vetorial 


Definição 1.2 (Produto escalar). Sejam U = (us,uy,uz) ev = (Vr,Uy,Vz) 
elementos do espaço Rê. Definimos o produto escalar de ú com Y como sendo 


> 


Definição 1.3 (Produto vetorial). Sejam U = (uz,uy,U:) €V = (Vr,Uy,Vz) 
elementos do espaço Rê. O produto vetorial de ú com Y é dado por 
Uxv= 


DD SisrtivjSe 


3; 23 
=1j=1k=1 


) 


onde e;;r é o símbolo de Levi-Civita, definido por 


1 se(i,j,k)e ((1,2,3),(2,3,1),(3,1,2)) 
cines —L se(,9,h)€4(1,9,2),(2,1,3),(3,2,1)) 
0 sei=9 ouj=koui=k 


Permutações pares dos índices do símbolo de Levi-Civita não alteram o sinal, 
ao contrário das permutações ímpares. Quaisquer índices iguais levam a zero. 
Desse modo, o produto escalar triplo ú - (V x w) é dado por 


Permutando os índices k <— i e em seguida k <— j, temos que ep;j = €jjk € 
então 


3 3/3 
u- (v x W) = z » a Cijk UV; Wk (3) 


A 


Sobre o símbolo de Levi-Civita, note que (€p)q = dpq implica 


3 


3 
E (gx) =D, >) cimndiidjmden = Eijk 
onde ô;; é o delta de Kronecker, definido como sendo 


fue O seil 
UC)1 sei=l 


Com isso, temos que 


3 3 
a Cijk€ilm = o é; . (6; x Ao É; . (é x em) = 38 x Si(é x em); 
i=1 i=1 ) 
= (6; x &) a (é x em) 


Como os versores cartesianos são ortonormais entre si, segue que o segundo 
membro só pode ser diferente de zero se 
é; x E = é x em ou é; x 6, = —é x em 


No primeiro caso obteremos 


3 
> EijkEilm = Ôj1Ôkm 
= 


e no segundo, 
3 


> ejnéitm = —ôjmôm 


i=1 


Como & x €& = O para qualquer p (isto é, os dois casos são mutualmente 
exclusivos), podemos somar ambos para obter a forma mais geral 


3 
a Cijk€ilm = Ôj1Ôkm — djmôki (4) 
i=1 


1.2.2 Operador nabla 


Definição 1.4 (Gradiente). Seja V(x,y,z): XC Rê R uma função cujas 
derivadas de primeira ordem existam em X. O gradiente de V em coordenadas 
cartesianas é definido pela expressão 


3 
VV(2,y,2) = > 0;V é (5) 
El 


Definição 1.5 (Divergente). Seja F: X CR? Rº um campo vetorial dado 
pela expressão 
F=F(F(x,y,2), Fo(x,y,2), Fa(x,y,2)) (6) 


cujas derivadas de primeira ordem existam em X. O divergente de F em coor- 
denadas cartesianas é definido por 


3 
V:F=50F; 
i=1 


Definição 1.6 (Rotacional). Seja F: X cRºs Rº um campo vetorial dado 
pela expressão (6). Definimos o rotacional de F em coordenadas cartesianas 
pela expressão 


Definição 1.7 (Laplaciano escalar e laplaciano vetorial). Sejam V: XCRº 
ReF:XcCRº-R funções contínuas com derivadas de primeira e segunda 
ordem contínuas. Então definimos o laplaciano escalar de V em coordenadas 
cartesianas pela expressão 


3 
VV=V(VV)-5,00V 
i=1 
e o laplaciano vetorial de F' por 


3 
VºF = 5(VºFi)é 


i=1 


Teorema 1.8 (Aplicações únicas de nabla). Sejam f,g: Rê» R funções com 
derivadas de primeira ordem contínuas e U,V campos vetoriais R$ > Rº com 
primeiras derivadas contínuas. Então 


1 V(fg)=[Vg+gVf 
D=úix(VxW)rvix(Vxw+(i-V)v+ (vv) 


x (fu) =(V)xui+fV xa) 
e Vx(ixD= VD) AV D+ V)di-(d.V) 


Teorema 1.9 (Aplicações sucessivas de nabla). Seja f:Rº->R uma função 
real com derivadas primeira e segunda contínuas e U um campo vetorial Rê >» Rº 
com primeira e segunda derivadas contínuas. Então 


L.Vx(V$=0 
2V(Vxi)=0 
3a. Vx(VxD)=V(V-ú)- Vi 


Demonstração. Para o primeiro item, temos 


3 3 3 
= » > >, EijhÓiO;V Ex 


Trocando os índices mudos à por j e vice-versa, obtemos 


3 338 3 3 3 
ba E: ba Ejih0,0;V Ex = — > x A Eijh0iO;V Ex 


i=1j=1 k=1 i=1j=1k=1 
="Vx(Vf 
onde usamos o fato dos operadores 0; e d; comutarem em funções com derivadas 


de primeira e segunda ordem contínuas e que uma permutação ímpar de €;jk 
acarreta na mudança de sinal. Por consequência, temos 


x (VÍ) 


Vx(VhO=o 
Para o segundo item, temos 
3 3 3.3 
V. (V xu) ô(V x De= 04), >; erra; 
k=1 k=10 i=lj= 
333 
=5,5, >, eyrôrôiu; 
i=19=1k=1 
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Trocando de lugar os índices mudos k e à e usando o mesmo raciocínio do 
primeiro item, chegaremos à conclusão que 


Vo(VxD)=-Vo(Vxi) => V(Vxu)=0 


E quanto ao último item, temos 


EsjrOG0V x U);&k 


Mu 


Vx(Vxm)=> 


i= 


m 
= 
I 


1 


3 3 
pa bo EijhOi (Eimj O Um) Ei k 


Il 
Me 
E TM 
Du 


1=1 k=1i=Im=1 
3 
1 
SS > Eijketmj | Dirt 
1=1 k l=1Im=1 


Usando a propriedade (4) do símbolo de Levi-Civita e sabendo que e;jk = Ejhi 
€ Eimj = €jlm, temos que 


3 


— ja A 
Vx(Vxi)= 5) (ôpidim — dm) 050 rum Se 
vkl;m=1 
3 3 
Aa o e 
ba Ôp10 im 0; 0VUm Elk — pa Ô mi 0,0 NUm Sk 
vkl;m=1 vkil;m=1 


3 0; OU ae 3 0; 0; UE 


vk=1 i,k=1 


Como U é um campo vetorial contínuo, o operador 0; comuta com à; e então 


3 3 3 3 
Vx(Vxm)=5 O [> cus) &— > [> com) Sk 


k=1 =1 k=1 Ni=1 
3 3 

e Ok(V - u)éx — >v up) Ek 
k=1 k=1 

=V(V ii) - Vi 


Como consequência, também podemos escrever o laplaciano vetorial usando 
a equação 
Vi=V(V.)-Vx(Vxi) 


Definição 1.8 (Conjunto estrelado). Uma região C é denominado de estrelada 
se existir algum ponto P em C tal que qualquer reta que ligue P a qualquer 
outro ponto de C esteja inteiramente em C. 


qu 


Teorema 1.10 (Campos irrotacionais). Um campo vetorial com derivadas de 
primeira e segunda ordem contínuas em um subconjunto estrelado X = Rº, 


F:X- Rº, definido em (6) satisfaz 
VxF(r)=0 


para todo r em X se e somente se existir alguma função vetorial V tal que 
Vre X o campo vetorial F pode ser escrito na forma 


F = VV 


Demonstração. Seja Vo(r) uma função definida pela expressão 


Vo(r) = Di Fo(t-r)-rdt 
com 
Fo(t-r) = F(t(r + ro)) 


onde ry é um vetor que aponte para algum ponto P em X tal que qualquer 
qualquer reta unindo P a qualquer outro ponto em X pertence inteiramente a 
X. Por X ser um conjunto estrelado, existe pelo menos um ponto de X que 
satisfaz essa condição. 

Seja também Xç a translação da região X pelo vetor ro até a origem (de 
forma que O e X9). Dessa forma, qualquer vetor posição r € Xo obedecerá 
t(r+ro)e X Vte [0,1]. 

Calculando o gradiente de Vo(r), temos 


VVo(r) = v/ Fo(t-r)-rdt 


0 
= [ 9 (Bot-s)-sat 


uma vez que a integral é em t e os limites de integração não possuem dependência 
em x, y ou 2. 
Usando o item 2 do teorema (1.8), iremos obter 
1 
VVo(r) = / Fo(tir)x(Vxr)+rx(Vx Fo(t-r)+ 
0 
(Po(t-x): V)r+ (x: V)Fo(t-r) dt 

Do primeiro rotacional temos 


VxFo(tir)=VxF(tr+ro)) 
= 0 
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para qualquer t e [0,1] já que qualquer ponto da reta t(r + ro) nesse domínio 
sempre pertencerá a X devido a ser um conjunto estrelado. Além disso, V xr = 
O, no que implica 


VVo(r) = ! (Po(t-r)-V)r+ (rr: V)Fo(t-r) dt 
(Fo(t - r) de rs e;6j= dudu Fidist = A Fis no Fo(t é r) 


(rr: V)Fo(t-r) = du dci 
or, (te;) 
= Ads idltes) e Si 
F; E 
= a (asi) e; 


Mas também 


d 0F; dei 
Fra Ren,  Ó(tei) | po a Fe) 


i) 
e portanto 


d 
dt 


F. 
(r- V)Fo(t-r) =t5) E &=t:Fo(t-r) 
j 


Consequentemente o gradiente de Wo fica 


1 
VVo(r) = ii Fo(t-r)+t- CRo(t-r) dt 
0 


Usando integração por partes no segundo termo, temos 


|, Fo(t-r) dt = Fo(r) 


VVo(r) = |, Fo(t E r) dt + tFo(t . r) 


ou seja, 
VVo(r) = F(r + ro) 


Agora definindo uma outra função V(r) tal que V(r + ro) = Vo(r) para todo 
r € Xo, temos que 
VV(r + ro) = F(r + ro) 


13 


Isto é, para qualquer vetor posição r tal que re X teremos 


VV(r) = F(r) 
A recíproca é trivial, uma vez que 
Vx(VV)=0 
devido ao primeiro item do teorema (1.9). E 


Um campo vetorial F' tal que V x F = 0 em um subconjunto X é denominado 
de irrotacional em X. E se F(r) puder ser totalmente descrito meramente 
pelo gradiente de alguma função V(r) em um conjunto X, então ele é dito 
conservativo em X. Além disso, a família dessas funções V é denominada de 
funções potenciais de F. Se X for um conjunto estrelado no qual F é irrotacional, 
então o teorema (1.10) garante que F é conservativo nessa região. 


Teorema 1.11 (Campos solenoidais). Um campo vetorial com primeira e se- 
gunda derivadas contínuas em um subconjunto estrelado XCRº,F:X- Rº, 
definido em (6) satisfaz 

V.F(r)=0 


para todo r em um subconjunto X = Rº se e somente se existir algum campo 
vetorial A tal que Vre X o campo vetorial F pode ser escrito na forma 


F=VxA 


Demonstração. Seja Ao(r) um campo vetorial definido pela expressão 


Ao(r) -[ Fo(t-r)x (t-r) dt 


onde novamente 
Fo(t-r) = F(t(r+ r0)) 


com r restrito ao conjunto Xy que é a translação de X por um vetor ro que 
aponta para um ponto P em X no qual qualquer reta que ligue P a qualquer 
outro em X esteja totalmente em X. 

Calculando o rotacional de V(r), temos 


Vx Ao) =Vx Fo(t-r) x (t-r) dt 


- [ V x (Fo(t-r) x (t-r)) dt 


uma vez que a integral é independente das coordenadas. 
Usando o sexto item do teorema (1.8), obteremos 


V x Ao(r) - [ Folt NV (t:2)) = (E -2)(V- Folt-r)+ 
((t-7)- V)Fo(t- x) — (Fo(t-r)- V)(t-r) dt 


14 


Como V-(t-r)=3teV-Fo(t-r) = 0, então 


V x Ao(r) = / 3tFo(t-x) + t(r: V)Fo(t-x) — H(Fo(t-r) - V)r dt 


Conforme calculado na demonstração do teorema (1.10), temos 

d 

(xr: V)Fo(t-r) = taFolt 1) e (Fo(t-r)-V)r=Fo(t-x) 

e então 

E d 

V x Ao(r) = / 3tFo(t-r) + PF ot -r) — tFo(t-r) dt 
0 
dd 
= / -Fo(t-r) + 2tFo(t-r) dt 

Agora note que 


d d 
PF ot -r) +2tFo(t-r) = ar (Fo(t-r)) 
e portanto 


1 


V x Ao(r) = À e (E Fo(t-r)) dt = [PFo(t-r)]| = Fo(r) 


0 
Definindo um outro campo vetorial A(r) tal que A(r + ro) = Ao(r), temos que 
VxA(r+ro)=F(r+ro) 


Ou seja, para qualquer vetor posição r que aponte para um ponto pertencente 
a X temos 
F=VxA 


A recíproca é trivial, uma vez que 
Vo(VxA)=0 
devido ao segundo item do teorema (1.9). E 


Um campo vetorial F no qual V-F = 0 em um subconjunto X é chamado 
de solenoidal nessa região. O teorema (1.11) garante que se X for um conjunto 
estrelado, então existe algum campo vetorial A tal que F = V x A em todo o 
conjunto X. 
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1.2.3 Integrais de linha 


Teorema 1.12. Se F é uma função vetorial integrável em X tal que para toda 
curva fechada C contida em X for verdade 


pro 

Cc 

! pdr= [ F.dr 
4 Cs 


para quaisquer curvas C4 e Cy em X que comecem em um mesmo ponto e 
terminem em outro ponto. A recíproca também é verdadeira. 


então 


Demonstração. Seja Cy qualquer curva que comece em um ponto p e termine 
em q, ambos em X. Seja também uma outra curva qualquer C5 que comece em 
q e termine em p. Então 


À Pdr+ [ F-dr=0 P-dr=0 
CG Cs, 6) 


Invertendo o sentido de €5 e obtendo uma nova curva Cs, invertemos o sinal da 


integral, no que resulta 
F pdr- [ F.dr=0 
for Co, 


j Pdr= [ F.dr 
C C 


Usando o mesmo raciocínio de trás para frente, pode-se provar a recíproca. HE 


1.2.4 Teoremas fundamentais 


Teorema 1.13 (Teorema do gradiente). Seja V: XCRº>R uma função 
diferenciável e C uma curva em X tal que comece em um ponto p, termine em 
aq e seja suave por partes. Então 


Í; VV(r)-dr = V(q) — V(p) 


Demonstração. Sejam C; curvas retas no R? alinhadas de modo que C; começa 
no término C;., e todos os pontos de começo e fim (denotados por x;, y; € Zi) 
estão em C. Então podemos aproximar a integral de linha por 


Lva ar VV(r)- dr 
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onde dr = dgr&k, de modo que 


vv(r)- de= “o dr + dy + dz = dV 
VÁ 


e então a 
VV(r)-dr= ) dV 
/. 3 Ci; 


Mas no limite infinitesimal 


(xiyizi) 
/ dV = / 1IdV = Vw; vi, Zi) = Vlgio,Yi-1,Zi-1) 
Ci; (x;-1,Yi-1,2:-1) 
Segue disso que 
n 
K VV(r)-dr= a Vai, Yi 2) — V(mia, Via, 201) 

C i 

Mas isso é uma soma telescópica, no que 


j VV(r)-dra V(g) — V(p) 
E 
No limite em que n tende ao infinito (refinar o conjunto das curvas C;), temos 


[vio ar= via vt) 


Corolário. Note que não foi preciso refinar o conjunto de curvas C; para obter 
o resultado correto, ou seja a integral de linha poderia ter sido calculada em 
uma linha reta que une p e q. Além disso a integral de linha não depende do 
caminho mas exclusivamente dos pontos inicial e final. Como consequência do 
teorema (1.12) segue que 
VV(rx)-dr=0 
6) 
para qualquer curva fechada C onde VV(x) está definido. 


Os seguintes teoremas serão provados apenas para superfícies ou volumes 
retangulares, porém possuem validade para formas mais gerais. 


Teorema 1.14 (Teorema de Green). Sejam P(x,y) e Q(x,y) funções definidas 
em uma região fechada 3 com fronteira OX e que possuem derivadas parciais de 
primeira ordem contínuas nessa superfície. Então 


“mtoq op 
PP) do + Olou) dy = | (SS de dy 


com a integral de linha orientada no sentido anti-horário. 
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Demonstração. Dada uma superfície retangular X no plano xy dada por [xo, £n|x 
[yo; Ym], podemos dividi-la em regiões retangulares menores X;; contidas em D. 
A soma de todas as integrais de linha que percorrem cada região 5;; devem se 
cancelar exceto nas bordas de X, no que implica 


Ô P(x, 8) da + Q(x,1) dy = 515) o) Plsujde O(A dy 
o TT Ei 


onde 
Li Yj 
P Pede + Ole )dy= [ Plempado + Oleo) dy+ 
Eij Ti—1 Yj-1 
Ti—1 Yj-1 
[o Pleupars / Orma)ay 
Di Yj 


com a ordem das integrais começando de x; 4 até x; com y = y;-1 e percorrendo 
no sentido anti-horário, ou seja z;-1| < Z; CYj-1 <Yj- 

Agora suponha que a grelha seja refinada o suficiente para que os valores de 
x; sejam muito próximos de x;-1 e o mesmo para y; e y;-1. Então na primeira 
integral podemos fazer uma aproximação de P(x,y;-1) por série de Taylor? até 
a primeira ordem, isto é 


oP 
P(x,y;-1) = Pl(x;oa,yj-1) + EE (x — x;-1) + O(x?) 


Ti—l 


onde O(x2) são os termos que dependem de x? ou potências maiores. 
Denotando Ax; = x; — x;-1, temos que 


E Ax;)2 OP 
K P(x,y;-1) de = AxiP(a;,y;-1) + o de 


i— 


+ O((Ax;)*) 


Til 


Analogamente temos que 


vi Ay)? O 
Qin) dy = Am O(ossuj-1) + CUL + o(a) 
Ui =1 Yj-1 
Lia. A i 2 oP 
o P(z,y;) de = —Ax;P(a; 1,9) ( o a O((Ax;)?) 
Yj-1 A N2 Fa) 
! Q(xi-,y) dy = —Ay;Q(xi-,Yj-1) um e O((Ay;)*) 
Yj Yj—1 


2Em toda a demonstração, todas as expansões serão centradas no ponto (xi, yj-1)- 
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e então 


Ê Pl; 1) de + Q(e,9) dy = Sai Plecas tg) = Plre volt 


ij 


2 | | + As [Ola vp=a)— Qlra U4=4)] + 
(Ay;)? q oq 3 3 
Ag; Ay; 
Yj-—1 Yj-1 


Ê P(e;1) de + Q(e,9) dy = Sei P(erastg-) = Pla, yo] + 


Ay; [Q(ei, ty;-1) — QUxio 1, y-1)] + O(A), (Ay)? 
Mas note que no limite infinitesimal, temos que 


2 Pl =P =) 
Oy Ay 


no que implica 
oP 
Oy 


onde a derivada parcial é calculada no ponto (x;-1,Yyi;-1). Analogamente temos 


Az; [P(x;,y;-1) — Pl(x;-1,9;)]) = ——— Az; Ay; 


oq 
Ay; (Q(x, y;-) — Qlxi-,y;-1)] = dy ATA; 
Ou seja, a integral de linha em cada retângulo menor 5;; é dado por 
oq) oP 
À Plaado+ O(oa) dy = (SED) dsãy; + Oo)? (05) 
Dij L y 


No limite em que Ax; > 0 e Ay; — 0 (e consequentemente n e m tendem ao 
infinito), o termo O((Ax;)?, (Ay;)*) tende a zero rapidamente. Agora seja dx, 
ôy definidos por 


dx = max(Ax;, Vi) óy = max(Ay;, Vj) 


Então temos 
Gm (00 OP 
P = li ) ) — — — | Ax;Ay; 
Ê. (x,y) de + Q(x,y) dy lim 2 (5 =) ay 


Porém o membro direito é exatamente a definição de integral dupla na superfície 


X,ou seja 
oQ OP 
Ê. Pass Odo IA ( - —) dd 
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Teorema 1.15 (Teorema de Gauss da divergência). Seja V E Rº um conjunto 
compacto que possui uma fronteira suave e definida por partes OV. Se F for 
uma função continuamente diferenciável em V e suas vizinhanças, então 


dh P-da— [) V-.Fdr 
ov Vv 


onde o vetor da é normal à superfície OV e orientado para fora do volume. 


Demonstração. Seja V um volume retangular dado por [x0,%n] x [yo;Ym] x 
[20,2]. De maneira análoga à demonstração do teorema de Green, pode- 
mos repartir esse volume V em vários paralelepípedos menores Vi;; dados por 
[u;-1,2;] x [yj-1,y5] x [2h-1, 2%). Às integrais de superfície desses volumes me- 
nores se anulam em diversas faces devido à orientação do vetor área, exceto as 
da fronteira OV. Então 


nm 
dp Fda-3555db Foda 
av à jd ok HVijk 
Cada volume V;;, é um paralelepípedo de 6 faces, então 


dh pda [| P-da+ |) F-da+ |) F. da+ 
Vijk Ses Sica s 


í Zk—1 


IA P-da+ |) F-da+ |) F.da 
Sa; Sy; Son 


Agora denotando as componentes de F por PQ e R, 
F=F(P(x,y,2),Q(x,9,2), R(x,y,2)) (7) 
então na primeira integral temos da = — dx di 
A P-da=— // P(x;-1,y,2) dydz 
Se; Se; 


: DR duo ar ds dA, o 
Agora sejam y; e z,. os pontos médios que delimitam Ay; = y; — yj-1 € Az, = 
Zk — 2Zk-1 respectivamente, ou seja 


Iremos expandir P para o centro do quadrado [y;-1,y;] X [2;-1, 2»). Primeira- 
mente, na variável y temos 


A 


0 
P(x;-1,9,2) = P(g;-1, 4,52) + oy (y ai Ui) + O((y E y;)?) 


1 


Yj 
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E agora expandindo em z, obteremos 


oP oP 
P(£ta1,4),2) — Pla 1,9. 2%) + Ôz « 2h) os Oy ; (y Y;)+ 
Yj Yj 
Pp 
poa UT UNE) + Ol abs (e— 4) 


Ao integrar P(x;.1,y,2) de yj-1 até y; e de z; 1 até 2», teremos que 


Yj Zk 
A (y-y;)dy=0 e / (2 — 2) dy = 0 
7 Zk— 


j—1 


devido à imparidade do integrando em relação ao centro do intervalo de inte- 


gração. 
No limite infinitesimal onde Ay; e Azy tendem a zero, o termo O(Ayi, Azz) 


tende a zero rapidamente. Nessas condições, isso implica 


A F.da= —P(x,1,Yj, 2) Ay; A zy 
Sejm 


A integral da superfície oposta, de maneira análoga, é dada por 


A F.da = P(z,, Yi» 2%) Ay; Ã za 
So; 


de modo que 


A F.da+ A F.da= (Pla; y), do) — P(ax;a, Yi» Zk)) Ay;Ã£Zk 
Sai Serio 


Mas no limite infinitesimal, temos 


oP 
Plax, y;, 2h) — P(xi, 9; 2) = dg Ati 


. . Z / 1 
onde a derivada parcial é calculada no ponto (x;,y), 24). Uma vez que Ay; e 
Az, tendem a zero e P tem derivadas contínuas, o valor da derivada parcial 


também tenderá à derivada no ponto (x;,y;, Zk). Isso significa que 
oP 


A F-das |) F.da- — Azx;Ay;Ã£Zy; 
So; Se; da 


De maneira totalmente análoga para as demais quatro superfícies, temos que 


A peda+ |) F.da= DNIT A 
Suja Syj oy 
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IAN F. da+ fo F. da Aridi Ãa 


oP q 
Pp. F.da= (+ a E =) Ax;Ay;Ãz,; = (V- F)Az;Ay;Ãz 


e então 


Novamente denotando dx, dy e dz os maiores valores de cada Ax;, Ay; e Azy 
respectivamente, temos 


pb. F.da= 0) (V- F)Az;Ay; Azy 


Porém o membro direito é justamente a definição de integral tripla em todo o 
volume V, no que acarreta 


dh P-da- [) V-Fdr 
ovV Vv 


Teorema 1.16 (Teorema de Stokes). Seja D uma superfície suave orientada no 
Rê com fronteira 05. Se F(r) for um campo vetorial dado por (7) definido em 
X e possuir derivadas parciais contínuas em uma região que contém 3, então 


) VxF-da- 6 F.dr 
õ FD) 


Demonstração. Podemos repartir a superfície » em várias outras menores 5, 
que pode ser descrita por uma parametrização r(u,v), onde 


r=r(ax(u,v),y(u,v),2(u,v)) 


com (w,v) pertencendo a uma região IL,g no plano uv. 
Para parametrizar a fronteira de X,q, suponha que u = u(t) e v = v(t), de 
forma que descrevam OII e também 


e a Emo das | dy; dz 
dr = dai + dy) + dzk (is Sis 7 k) ar 


die ôvdu | dx du E ôydu | ôydv 1 dz du dz du Ela: 
R du dt Ov dt E du dt Ov dt J du dt Ou dt 


e então temos que 


F dr = Pla(u(t), o(b), y(u(t), v(t)), 2(u(t), v(t))) fr e o 2] gi 
O(a(u(t), v(t), ylu(t), v(t)), 2(u(t), v(t)) (5 » da) E 


Real (O) luto (0) luto voto) (SE + ES) de 
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dx 


ê oz) d 
F.dr = (Pon + Ae,y,2) 5 o Re) e d+ 


dx To) Õz 
(PionoZ + ey, 2) 5 + Ren) 


F.dr=P(u,v) du + O(u,v) dv 


onde 
P(u,v) = Pig; y, e +” Q(a, Yy;, o E R(a, y, e 
UU UU UU 
o foh Oy z 
O(u,v) = Elo; y, ão Eur Q(a, 1) ão Es R(a, 12) 50 


com xr =a(u,v),y=y(u,v)ez=az(u,v). Segue disso que 


Ê. 


“A pq 


F.dr= , P(u,v) du + O(u, v) dv 
Opa 


dv 
E dt 


Por hipótese, as funções P, Q e R são funções deriváveis. Idem para as para- 
metrizações de x, y e z em função de u e v, assim como estas em função de 
t. Portanto podemos usar o teorema de Green (1.14) na integral de linha do 
membro direito com a orientação tal que o vetor de área infinitesimal da seja 


du x dv e então 


o) pdr= |) (E) duas 
0Lpa Hp Eu du 


Agora note que podemos reescrever P e Q na forma 


Puv)=F. e Au)=F.5 
de forma que 

2 Ç pr sed es 

du Ou Ou Ou Ou OuOv 


2 
oP OP é q dr 


Ou dv Ou “Ovôu 
e então 
009 OP OF dr E 02r oF dr 02r 
ou Ou Ou Ov Oudv dv Ou Ovou 
Na st ; 02r 02r 
Como a parametrização de r em u e v é suave, temos —— = —— 
Ovou OuOv 
quentemente 


q oP OF Or oF Or 
ou Ou Odu dv dv Odu 
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(8) 


e conse- 


Mas também, usando a regra da cadeia e denotando F, = P(u,v), F) = Q(u,v) 
e E; = R(u,v), obteremos 


oF E OF; 5 ê 0F; Oq; se 
ou 7 du É 2 e 


Ou Og; Ou 
e também 
or OG 5 Õr E ÔMh 
Ou du Ou E dv E 


Portanto temos 
oF Or 0F; Oq E Ogk 6 OF; Oq; Oqk E: 
du o Oq; Ou ) Ra ) - 3 Og; Ou Ou SAERR 


oF Or OF; Oqi Cah p OF; Oq; Oq 
Ou Ou -> Og; Ou ao -> Oq; Ou Ov 


De maneira totalmente análoga, obteremos 


0F or =D Oqi Oq; 


Ou Og; dv Ou 


de forma que 


oF Or OF 
Odu Ov Ou 


r -> OF; Oq; Oq 3 OF; Oq; Oq; 


0 
ou Og; Ou Ov Og; Ov Ou 


3 
E nm [Ôn0g Oq; Og; 
= 2d) (sa dv dv Ou 


ij 
A equação (8) fica? 


a Oqi dg; Oq Og; 
Po F.dr= JA sy 0;F';) (5 O ARCA du dv (9) 


Hpa ij 


Por outro lado, o vetor de área infinitesimal da na superfície X,y é dado 
pelo produto vetorial da variação do vetor posição r sob um acréscimo du com 
a variação de r sob um acréscimo dv, com a ordem do produto vetorial escolhida 


30 termo entre parênteses com as derivadas parciais em u e v é o jacobiano das coordenadas 
u e v em relação às coordenadas cartesianas q; e q;. Isso será visto na seção (1.3). 
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com base na orientação escolhida da superfície devido ao teorema de Green. Isso 
significa que 


or or 
or Or É or or 
= Ao A. dv = mn | 5 SA 'n d 
da ( E x =) du (D) as ae du dv 
Portanto 
E A E or or A 
VxF-da= du sra : A Eimn | 5a Em Mm du dv 
É or or R 
= > EijkEimn (0;F';) | — — |) &- eEndudo 
= Odu), NOV), 
igkimn 
3 
= e EijkEimn (0:F;) (=) (5) ôkn du dv 
É ou v 
ijkimn À m 
ro) 


3 A 
= pa EijkElmk (0;F';) (=) (5) du dv 
l m 


ijkim 
Usando a equação (4), temos 


3 3 
Ba EijkElmk — Bs EkijEkim — Ôjldjm — Ôimôjl 


k=1 k=1 


e portanto 


3 2 
or or 
VxF.da= 5. (Dida —OsmOGt) (0:F;) (5), (57) dude 


ijim 


3 A 
) m 


(0) 


3 
or or 
x» dimôji (0; F;) (=) (5) dudo 


ijim 


igim 


3 A, A 
VxF-da=5 (05) (=) | (=) dudo— 
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3 So , 3 . . 
V x F.da=5(0,F;) E qudv — (Gif) GR quer 


3 Ee] 


VxF-da=5 (05) (5 dg Oq; E) du dy 


E Ou Ou Ou Ou 


Integrando ambos os lados em u e v na superfície Ig, obteremos 


A VxFida- |) Sar) Cm da DUO a 0) 
E dl Ea du Ov dv Ou 


Pq ij 


Agora compare o membro da direita da equação (9) com o da equação (10). 


Portanto 
, F.dr = A VxF.da 
Fab) Hpa 


“A pq 
Mas integrar u e v na superfície IL, do plano uv é totalmente equivalente a 
integrar na superfície 5,y no espaço das coordenadas cartesianas, ou seja 


) Pdr= |) VxF.da (11) 
oz Diga 


“A pq 
Para obtermos a integração em toda a superfície D, basta somarmos (11) em 
peq. À integração nas fronteiras de L,y se anularão exceto nas bordas de » 
enquanto que a soma das integrais das superfícies », será a integral de toda a 


superfície X, resultando 
2 pdr= |) vxE-da 
o zL 


1.2.5 Integração por partes 


Teorema 1.17. Sejam A, B e f funções contínuas com derivadas contínuas. 
Então é verdade que 


1. hf tis nar = db sacas ff acc nar 
2. Lica | tax(vn-da+ db s4-ar 
E Bv war = ff acyxBjar+ 9) (ax) da 


Demonstração. Do terceiro item do teorema (1.8), temos que 


V=(FA) = AV A) + As(VI) 
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Integrando ambos os membros em um volume V obteremos 


hs radar = ff rev -apar+ ff accvnar 


Usando o teorema de Gauss (1.15) da divergência no primeiro membro, ficamos 


Co dEisacãa fl nvcazars fl acconas 
hf tv nar = db sacas ff accvnar 


o que demonstra o primeiro item. 
Do quinto item do teorema (1.8), temos 


e então 


Vx(JA)=(VOxXxA+Sf(VxA) 


Integrando ambos os lados em uma superfície 3, obteremos 


fiveua-aa= fl (vu a das ff tv xa) da 


Usando o teorema de Stokes (1.16) no primeiro membro, chegaremos na ex- 


pressão 
Potacar= ff (vixA-das fl 17x A)-da 


Segue disso que 


fl tio x a) -da -— fl (VÍ XA- das do fA dr 
-| [A x (VA) da + É fA-dr 
õ oz 


o que prova o segundo item. 
Do quarto item do teorema (1.8), temos 


V(AxB)-=B(VxA)J-A(VxB) 


Integrando ambos os lados em um volume V e usando o teorema de Gauss (1.15) 
da divergência no primeiro membro, obteremos 


db ta xB)-aa= ff] B.(v x A) ar= ff a (V xB)d 


e então 


hs (Ux Adr fl A (VxBjar+ 9) (AxB)-da 


o que demonstra o terceiro item. E 
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1.3 Jacobiano 
1.3.1 Áreas e volumes infinitesimais sob parametrização 


Sejam u e v parametrizações suaves? das coordenadas cartesianas «x e 4 no plano 
ay. Podemos expandir x e y em série de Taylor em termos de ue v a partir 
de um ponto (u;-1,v;-1) para saber o comportamento dessas coordenadas sob 
uma pequena variação de u e v. Para a coordenada x temos 


ox ox 
eu v)=e(wva)+ 5) (u-u)+ so (ova)+ 
Ui-1 Ui—1 
Vj-1 vVj-1 
2x Ao 
NOR (u— u-)(v — v;-1) + O(uê, ví) 
e 
e para a coordenada y, 
9) 9) 
y(u,v) = y(u;-, Vj-1) + Em 1 = Ui) + ao no — vj-1)+ 
Vj-1 Vj-l 
0? 
os (u— ua)(v — v;-1) + O(u2, v2) 
Ui—1 


Vj-1 


Agora fixemos v — v;-1. Variando u de u;-1 para u;, as coordenadas x e y 
também irão variar. No caso de x, temos 


a(u;, Vj-1) E gui, vj-1) = Zu (ui; = Ui) + O(u2,v?) 
Uluioa 
Vj-1 
ox 23 
a(u;,vj-1) = 2(u-,vj-1) + Zu Au + O(ul,v”) 
Ului- 
Vj—-1 
de modo que 
ox 923 
Ay, = (ui, vj-1) E e(u;m,Vj-1) = — Au + O(u UV ) 
ou Ui—-1 
Vj-1 
Similarmente, para y temos 
Õ 2,2 
Ayu = ylu va) — y(u-,vja)= |  Au+ Ou", vê) 
OU lusa 


Vj-1 


4Parametrizações suaves são aquelas cujas derivadas parciais em u e v existam. 
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No limite em que Au é infinitesimal, o termo O(u2, v2) rapidamente tenderá 
a zero. Então o vetor ray que liga os pontos em que houve uma variação 
infinitesimal Au é dado por 


De modo totalmente análogo, o vetor TA» que liga os pontos em que houve uma 
variação infinitesimal Av é dado por 


Tau = (= 5) Au 


dv Ou 


A área À de um paralelogramo formado por dois vetores à e b é numericamente 
igual ao volume de um paralelepípedo cuja base é o paralelogramo formado pelos 
vetores à e b com altura dada pelo versor k. Isto é o produto escalar triplo 


; a dr ay O es a 
A=[k-ãx5|=|det| be by O - [aee | q 5º | 
0 01 = “y 


com o módulo aparecendo devido que a área A é não-negativa. 
Portanto a área infinitesimal da de uma região formada no plano xy devido 
à uma região retangular infinitesimal no plano uv é dada por” 


du Ou 

da = |det du dv 
do Oy 
dv Ou 


Mas como da = da dy no plano cartesiano, temos 


ca Oy 
du Ou 

dx dy = |det du dv 
da O 
dv Ou 


Agora no caso tridimensional suponha uma parametrização suave do tipo 
e(u, v,w), y(u,v,w) e z(u,v,w). Realizando o mesmo procedimento de expansão 
em série de Taylor, os vetores que ligam cada variação infinitesimal Au, Av ou 


5Os diferenciais aparecem fora do determinante devido à propriedade multiplicativa por 
escalar. 
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Aw são dados por 


O volume infinitesimal dV formado por esses três vetores é dado pelo produto 

escalar triplo 

dx Oy Oz 

du Ou Ou 

E 3; » dx Oy Oz 

dV = |Fau Tao xXxTaw)= |det| > <L EL du dv dw 
dv  Ov Ou 


dx Oy Oz 


dw dw dw 
Ás duas matrizes que apareceram nos motiva a fazermos a seguinte definição: 
Definição 1.9 (Matriz jacobiana). Sejam u ev parametrizações das coordena- 


das cartesianas x ey cujas derivadas parciais existam. Então denominamos o 
determinante da seguinte matriz de jacobiano: 


os Oy 
du O 
x,y) dd u cu Fobi êy Oy Ox (12) 
Ou, v) dx Oy dudv dudu 
dv Ou 


No caso tridimensional em que há parâmetros u, v e w, temos 


du Odu Ou 

Eee) = det ox Oy oz 

Ou, v, tw) dv dv Ou 
Go oi e 
dw Odw O 

Dessa forma temos que 
da = drdy = d(x, 1) dudv 
Ou, v) 
ou então no caso tridimensional 

dV = dxdydz = A 9,2) du dv dw 

(1,0, wW) 


1.3.2 Áreas e volumes orientados 


A área A e volume V de uma região parametrizada 3 no espaço cartesiano (ou 
a região II no espaço uvw) pode ser calculada integrando as diferenciais dx e 


dy, ou seja q j som f E 
V= fl arvar= ff [5 


Entretanto podemos atribuir sinais negativos as essas quantidades se levarmos 
em conta a orientação (ou ordem de escolha) dos vetores Tau, TAv € TAw, € Isso 
pode ser feito ao não aplicar o valor absoluto no jacobiano: 


A= |) aray= ff SE V) qud 

DD) un Ou,v 

Vu ficanas= | 5 Mm, 19,2) idas 
n O(u,v,wW) 


Denominamos À e V respectivamente de área e volume orientados. 


du dv 


Oa,y,2) 
Ou, v,wW) 


du dv dw 


Exemplo 1.1 (Troca das coordenadas cartesianas). Suponha a(u,v) = ve 
yu(u,v) = u. Então 


x,y) dxdy dydz 
u,v) dudu dudu 


A= |) dedy=— |) quav=— |) duas 


Isso significa que inverter a ordem de integração altera o sinal da área orientada. 


0-0-1.1=—1 


Portanto 


Exemplo 1.2 (Área de uma linha). Suponha a(u,v) = u ey(u,v) = u. Então 
x,y) dxdy dydz 


u,v) dudu dudu 


A= [| drdy= | odudo= |) 0aras =o 


Exemplo 1.3 (Área em coordenadas polares). Suponha x(r,0) = rcos(b) e 
y(r,0) = rsin(0). Então 


=1:0-1-0=0 


E então 


x,y) Odxdy dydz 


2 DB es 
27,0) * Dr 00 Dr 00 rcos"(0) + rsin“(0) =r 


Logo a área em coordenadas polares é dada por 


A= |) aray= |) rardo 
z m 
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1.4 Vetores contravariantes e covariantes 
1.4.1 Transformações contravariantes 


Sejam € = fu,ua,...sUnje 8 = fui,va,..., vn) duas bases de um espaço veto- 
rial V cujo corpo é K. Podemos escrever cada vetor de 8 como uma combinação 
linear dos elementos de é: 


vj= Dat (13) 


onde a é o elemento da i-ésima linha e j-ésima coluna de uma matriz 4. Em 
outras palavras, 4 é a matriz que leva a base é para 5. 
Seja y um elemento de V cujas coordenadas na base 8 são (yl,y2,...,y"), 


ou seja 
de . 
v= yo; 
J 


Segue disso que 


nm n n n 
es yiaju; = ) las Ui; = ) vu; 
j i à Nj F 


n 


onde xº = > y'ai é a i-ésima coordenada de y na base é. 
j - 
Usando a notação matricial, temos que 
x= Ay 


onde x e y são os vetores coluna contendo as coordenadas de y na base é e 
8 respectivamente. Sabemos? que det A + 0 e portanto A possui uma inversa 
A-1, de forma que 


ou equivalentemente 


onde à) é o elemento da i-ésima linha e j-ésima coluna da matriz 47. 

Agora note que a transformação das coordenadas do elemento y de E para 
8 se faz usando a matriz 4-1, apesar da base £ se transformar para 8 usando 
a matriz 4. Com isso, denominamos que as coordenadas do elemento %y se 
transformam contravariantemente sob uma mudança de base. 


êSe det A fosse nulo, então O seria um autovalor de A e portanto existiria algum autovetor 
coluna k, cujos elementos pertencem a K e não todos nulos, tal que Ak = 0. Ou seja, existiria 
alguma n-upla não nula na base 8 que seria o vetor nulo na base É (e consequentemente o 
elemento nulo de V), o que é falso. 
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1.4.2 Transformações covariantes 


Seja a:V > K um funcional linear que leva um elemento de V ao corpo K. 
No contexto da álgebra linear, esse funcional é uma transformação linear cujo 
contradomínio é o espaço vetorial formado pelo corpo K. Podemos representar 
a univocamente por meio de coordenadas em relação a uma base, de forma 
semelhante à coordenada de um vetor. Sejam p; as coordenadas de a na base 
É, isto é 
pi = a(u;) 
Então a aplicação de a no elemento q resultará 


a(y) = a [> cu) = Dyt'o(u;) = D'pi 


Por outro lado, substituindo xº pela sua representação na base 8, temos que 


a(Y) = XX gap; = >;y [> cs 
ij j i 


Se denotarmos q; = a(v;), temos também 


alg)=a Ê o) - Diyia(o;) = via 


Como a() é o mesmo independentemente da base escolhida, podemos igualar 
as duas expressões para obtermos 
il ls & ig 
So Bum) 0 
j i 
Essa equação é verdadeira para qualquer y de V (ou seja, é válida para qualquer 
escolha dos valores y'). Em particular, podemos escolher y' = d;; e então 
n n n n n 
o a = DE y' = E (a = Sin) djj! = Gp — Djaji =0 
j i j i i 


Trocando o índice mudo j' para j, obteremos 


n 
i 
q = e 4;Di 
i 
Na notação matricial, essa expressão é equivalente a 
q = pÁ 
onde p e q são vetores linha contendo as coordenadas de a na base € e 8 
respectivamente. 
Ágora note que a transformação das coordenadas do funcional linear a da 
base é para 8 se dá aplicando a matriz 4, exatamente a mesma matriz usada 


para levar a base é para 8. Nesse caso, denominamos que as coordenadas de q 
se transformam covariantemente sob uma mudança de base. 
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1.4.3 Vetores sob mudança de base 


Podemos separar os vetores em três grupos: um nos quais se transformam con- 
travariantemente, outro nos que se transformam covariantemente e um terceiro 
no qual são invariantes sob mudança de base. Isso nos motiva a fazer as seguintes 
definições. 


Definição 1.10 (Vetores contravariantes). Dada uma matriz A que leva uma 
base E para 8, um vetor cuja n-upla se transforme da base É para B sob a 
aplicação de A! é denominado de contravariante. As coordenadas de um vetor 
contravariante são denotadas com um índice superior. 


Exemplo 1.4 (Vetor posição). Seja y = (q!,q2,qº) um vetor qualquer do Rê, 
onde q” são as coordenadas na base canônica É formada por 


6 = ((1,0,0);(0,1,0);(0,0,1)) = téi, 62,635 
Então temos que à 
= bj q'& 
Agora seja a nova base 5 formada pelos elementos 
8 = ((2,0,0):(0,2,0):(0,0,2)) = (261, 262, 263) = (fi, fo, f5) 


isto é, dobramos o módulo de cada elemento de &. Consequentemente 


e portanto 
3 q ' 
= o gh 
KA 
com agora f* = 7 coordenadas de y na base 5. 


Note que aumentar o tamanho dos vetores da base canônica E implica em 
reduzir na mesma proporção as coordenadas para obtermos as novas coordenadas 
na base 8. Esse comportamento invertido se mostra pela matriz A que leva E 
em B: 


enquanto que a matriz que leva as coordenadas de y na base É para 5 é 


Eos sd cf 
Ad =| O 2 
O a, sigo 
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ou seja, 


f! q! 
f= Ala com f=| f2 q= | -gº 
Fe q” 


e consequentemente y é um vetor contravariante. 


Definição 1.11 (Vetores covariantes). Dada uma matriz A que leva uma base 
é para 8, um vetor cuja n-upla se transforme da base É para 8 sob a aplicação 
de A é denominado de covariante. As coordenadas de um vetor covariante são 
denotadas com um índice inferior. 


oV ov ov 
Oq!" Oq?" Oq? 
uma função V na base canônica E cujas derivadas existam, onde q' são as 
coordenadas nessa base. Usando a mesma base 8 do exemplo (1.4), a regra da 
cadeia nos diz que 


Exemplo 1.5 (Gradiente). Seja (VV): = ( ) o gradiente de 


oV.  OV of “10V 
dq Of dq 20f 
no que implica 
oV » 9%V 
of' oq 


Por outro lado, o gradiente de V na base B é 


oV oV ovV OV .OV .OV 


Oq!” Oq?” Og 
Note que aumentar o tamanho dos vetores da base canônica E implica em 
aumentar na mesma proporção as coordenadas para obtermos as novas coorde- 
nadas do gradiente na base 8. Na notação matricial, isso significa 


(VV)o = (VA 


com 


oV OV OV | 


e portanto VV é um vetor covariante. 
Usando a notação com índice inferior, temos 


yi = 2x; onde vy= 


Definição 1.12 (Escalares). Um vetor cuja n-upla é invariante sob qualquer 
mudança de base é denominado de escalar. 
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Exemplo 1.6 (Módulo de um vetor). Dado y e R$, o módulo quadrático de 
na base canônica É é dada por 


3 3 3 3 
Ill ="7= [i0a) [Des =D 1d é 
i j à jd 
3 3 


3 
7] = MATEO; = at) 


) 


Já o módulo de y numa outra base 8 qualquer é dado por 
: São fi E Se my dy 
Ialla="8" = [7a] [:78) =D, 1 É -fi 

j ij 


Mas também os elementos da base 8 podem ser escrito na base canônica É, 


de forma que 


3 E 
6-8 = (Dates) (Darei) = Datas cc 
k km 


o» 

eta» 

Il 
=D 
3 [Ie 

Q 

Sm 

[eu 

= 
3 

Il 
[Du 

Q 

Sa 

Q 

Sm 


e então 
3.3 800 E 6 E 
Ilê =D fifiatas =) [> rá) | (5 pos 
id k k i j 


Mas os coeficientes ak são exatamente os que transformam as coordenadas de 
y na base 8 para € (transformação contravariante), ou seja, 


e consequentemente 


3 


3 
Ilza = 50º = 50? = Ing 
k 


) 


Isso significa que a quantidade ||? é invariante sob a mudança de base da 

canônica para outra qualquer e consequentemente também é invariante sob mu- 
24 

dança entre quaisquer bases. Logo ||y|É é um escalar. O vetor (h1,...,hn) onde 

hi; são funções que dependem apenas de ||y|| também é um escalar. 
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1.5  Pseudovetores 


Suponha um operador R : Rº > Rº que reflete um vetor em relação à coorde- 
nada y, ou seja, cuja representação matricial é 


1 0 0 
R=|0 -1 0 
0 0 1 


Se à e Rº, então temos 
33 
Rã = bos Rij4;6i 
UR ai) 


onde Ri = (—D)'H6;;. 
Agora sejam q,b e R3. Segue que 


3 
(Rã) x (Rb) = 3. esje(Ra)i(RD) jk 


igk 
3 3 3 

1 1 ! A 

= er) RauAi Rjm Bm 
igk 7 m 


3 
= + en) oa(—1)H16;m A Bm ek 
ijkim 
3 . . 
= DD en AB; 
ijk 


isto é, o vetor (Ra) x (Rb) possui as componentes x e z invertidas em relação 
ao vetor à x b. Consequentemente temos 


. 100 , : 
(Ra) x (RD) = à a (1x5) =-r(ax5) 


que é diferente de R (à x b). Esse comportamento anômalo também fica nítido 
ao inverter o sentido de todas as componentes. Seja agora K a matriz 


—1. 0 0 
K 


I 
[am 
| 
H 
[am 


Temos que Kã = —d e portanto 


(Ka) x (KD) = (-a) x ( D-axbaK (5x5) 


ou seja, o elemento à x b foi invariante sob a inversão dos eixos. Por causa desse 
comportamento anômalo, denominamos à x b de pseudovetor (ou vetor axial). 
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2 Espaço dos p-vetores 


Definição 2.1 (Espaços 0-vetor e I-vetor). Seja V um espaço vetorial de di- 
mensão n cujo corpo é K. Define-se o espaço O-vetor A! V=K como sendo o 
corpo K. O espaço Í-vetor A V=Vé definido como sendo o próprio espaço 
vetorial V. 


Definição 2.2 (Espaço 2-vetor). Dado o espaço vetorial em (2.1), define-se 
AC V o espaço vetorial 2-vetor formado por todos os elementos u, denominados 
de bi-vetores, tais que podem ser escritos como 


U=UAY para algum v,yeV 


com a operação binária A : |V x V] > A? V, supondo x, y e z elementos de V 
eabeK, satisfazendo os dois axiomas: 


Axioma 1 (Linearidade à direita). (ax +by) na z=a(x an 2) +b(y A 2) 
Axioma 2 (Anticomutatividade). va y=-—yA £ 
Denomina-se « A y de produto exterior (ou externo) dos vetores x e y. 
Teorema 2.1 (Linearidade à esquerda). z A (ax + by) = a(z n x) +b(z A y) 
Demonstração. 
za (azx+by)=-—(ax+by)a z (Axioma 2) 

=-a(xnz)—b(ya 2) (Axioma 1) 
=a(zns)+b(z Ay) (Axioma 2) 


E 
Teorema 2.2. 2 A x = 0, com O sendo o elemento nulo de A” V. 
Demonstração. Devido ao axioma da anticomutatividade (2), temos 
LTALS-TIAL 
e então 
2(rn1)=0 => 2rAzx=0 
E 


Ágora sejam y e w elementos de V cujas representações na base é são 


n n 
Y =5,r'%; e w = Suwtz; 
F F 


O produto exterior de y com w é 


one (8) (0) -Bfirêno 


toi) 
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Devido ao teorema (2.2) e ao axioma (2), temos que os termos em que à = y 
se anulam e os que aparecem com à e j trocados podem ser agrupados com um 
sinal de menos, resultando 


n i—l 


yaw=5 5 (Pu! = iu; 25) (14) 


i=2j=1 


onde a soma é feita de forma que cada combinação de à e j seja feita apenas 
B É . - 2 

uma vez (1<j<i<n). Consequentemente a dimensão do espaço A“V, que 

é o número de elementos (x; A x;) linearmente independentes, é dada por 


dim AV = 5 5 1= 56-1)= MOD (1) 


Podemos generalizar criando um novo espaço vetorial AS V formado pelos 
elementos que são a aplicação do operador A entre um bi-vetor de AS Ve um 
vetor de V, resultando em um tri-vetor. 


Definição 2.3 (Espaço dos p-vetores). Seja V um espaço vetorial de dimensão 
n com corpo K. O espaço p-vetor APV, com 1 <p <n, é formado pela 
totalidade dos elementos u que podem ser escritos na forma 


E (a Apae ro e Aa] (15) 


com os p índices q; satisfazendo IL <q <q <...<qpanxmeVetal 
operação entre vetores de V satisfazendo 


ly nyrA... AYp troca de sinal se dois elementos y; forem trocados; 
2. WAY A... AY é o elemento nulo de APV sey; = y; para algum i £ j; 


3. (ax+by)nyoA...nyp=a(zaya... ayp)td(yagoA... Ayp),ea 
mesma propriedade é válida se y; for substituído por ax + by, com a,be K 
ex,yev. 


A dimensão do espaço A? V é o número de elementos linearmente indepen- 
dentes dados pelas combinações dos índices q; em (15), ou seja 


ve v-(0) 


Sep >n, então A? V é um espaço vetorial de um só elemento, o nulo, já que 
haverão n + 1 elementos linearmente dependentes e consequentemente o item 2 
da definição (2.3) se aplica. 

Se u for um elemento de A? V, então u pode ser escrito em termos de uma 
base & = fx1,...,4n) de V na forma 


Emi hiho...hn 
u= x a ER PR 
hiha...hp 
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onde os índices h; satisfazem 1 < hy <ho<...<hp<n. Ou seja, para 
cada conjunto H = (m,ho,...,hp) E (1,2,...,n) de p elementos temos uma 
coordenada at12:--Pm associada a Xh A Tha A... A Th, Podemos compactar 
tudo isso usando a notação 
Uu= 5 alzy 
H 


onde u é a soma sobre todas as possíveis combinações de índices H, com a 


afihadhm e rg = Tha A Tha A A Th. 


H — 


2.1 Determinantes 


Dado um espaço vetorial V de dimensão n e corpo K e uma transformação linear 
A:V + V, iremos definir uma função ga de n variáveis como 


DA Vs AV 
ga(ii,...,Un) = Axy A Axo A... A Ar, 


ondex;e Ve x”V éo conjunto das n-uplas contendo elementos de V (produto 
cartesiano). Equivalentemente podemos definir uma função fa de uma variável 


em A”V tal que 
tA:AVS AV 
fala nx2 A... An tn)= Axya Axo A... A Ara = ga(xi,..., Un) (16) 


A dimensão de A”V é 1, pois há somente uma única combinação de n 
índices em n elementos de qualquer base de V. Então suponha u,v e A”V, 
P = (6) uma base desse espaço no qual 


u=mê e v=hB& aBiekK q 0 


e que T(u) = v para alguma transformação linear T: A"V »> A”V. Temos 
então que 
Bi 


T(a16i) = a T(E) = Bão — T(6) = ma 


Isso significa que qualquer transformação linear em um espaço vetorial de uma 
dimensão para ele próprio é meramente uma multiplicação por escalar. Então 
podemos reescrever (16) como sendo 


Az A AxoA ...s Ar =9YA(L ALA... A Xn) (17) 


onde y4 € K. Agora suponha que os elementos x; formem uma base é de V e 
[az] seja a representação matricial de A nessa base, ou seja, 


n 
Az; = > io 
j 
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então 


Ar, A Axo A... A Áx, = (Beto) A [lote] ALA [Bota 
k q 


Í gb I (x; 
x ara cccal(L;ALEA A Lg) 


Como A”V é um espaço unidimensional, podemos escrever todos os termos 


TA Xp A... A tg em função de 7, A 72 A ... A Zn. Os índices repetidos 
resultarão em termos nulos e os sinais dependerão se a permutação dos índices 
9,k,...,q for par ou ímpar em relação à sequência 1,2,...,n. Podemos resumir 


tudo isso pela equação 


Ax, A ÁxoA...AÁZ, = > Ejt.g) -aR . eat) (24 A 29 A...A Zn) (18) 
Jk...q 


onde e;p....q é à generalização do símbolo de Levi-Civita para n índices. Compa- 
rando as expressões (17) e (18), temos que 


n 
ok 
YA = se Ejk..g0] 09 * 0. 08 (19) 
Jk...q 


A equação (19) coincide com uma das várias definições equivalentes de deter- 
minante da matriz 4. Se |A| é o determinante da matriz que representa a 
transformação linear A na base é, então y4 = [A]. Isso nos motiva a fazer a 
seguinte definição: 


Definição 2.4 (Determinante de transformações lineares). Seja V um espaço 
vetorial de dimensão n, corpo K ex;e V. Seja também A uma transformação 
linear de V para V. Denominamos |A| e K de determinante de A como sendo 
o escalar que relaciona 


Az, a AxsA ...n Ara =|Al(r Axo A... A 2) 
Usando a definição (2.4), podemos mostrar que 


IAB| (x, A 2324 ...n 2n) = ABr, An ABao A ...n ABax, 
= A(Br1) A A(Bas) A... A A(Ba,) 
=|A|(Br, A Bars A ...n Ban) 


=|AL-IB|(xy ng2A...A &n) 


Como (x, A 2x2 A ...A Zn) £ O se os elementos x; forem linearmente indepen- 
dentes, concluímos que 
|A4B| = |4/-1B] 
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Além disso, se 4, é a matriz que representa a transformação linear A na 
base 1 e 45 é a matriz que representa esse mesmo operador na base 2, temos 
que 

Ao = [o 4 [1] 


onde [7]) é a matriz mudança de base de i para j. Mas como ([1]5)-" = [7]%, 
segue que |49| = [134,01] = |Ai;|. Isso significa que o determinante de A é 
um escalar, pois é invariante sob mudança de base. 
2.2 Produto externo entre p-vetores 


Podemos relacionar a operação A vista na definição (2.2) com as operações A 
vistas na definição (2.3) como sendo o mesmo operador, onde agora será definido 


como 
p+q 


A Av X AVES AV 
e que satisfaz 
(xi nx A...AL)Aln AA AY) LADA ALLA AVDA AYg 
Segue imediatamente disso que 
1. O operador A é associativo, ou seja, A (yAZ)=2AyAZ=(ZAY)AZ. 


2. O operador A é comutativo se pg for par e é anticomutativo se pg for 
ímpar, ou seja, 2 A y=(—1)PIy a x. 


3. O operador A continua distributivo, ou seja, A (ytZ)=2Ay+HXAZ. 


O segundo item pode ser rapidamente provado ao mover cada y; para à esquerda, 
começando com y até yq, introduzindo um fator (—1)? a cada operação: 


(x AXoA...AL)A(NAgDA.. AY) = 


(IDP AZIAZIA ATA VIA A Ya 


(LADA AL)A(N AGA. AY) 


(IDP) AyATIATIA ATP AYIA A Va 


(x AXA..AL)A(NAgDA.. AY) = 


(=) A dy Ars Do A A Bo AE 


(x AXºA..AL)A(NADA.. AY) = 


(DI (ynanya. Ayg)A(z atos... A Lp) 
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2.3 Produto interno 
2.3.1 Produto entre 1-vetores 


Ágora suponha que o espaço vetorial V esteja munido de um produto interno, 
operação binária V x V > K que satisfaça os seguintes axiomas: 


Axioma 3 (Linearidade à direita). (z,ax + by) = alz,u) + b(z,y) Va,y,2 E 
VoabeK. 


Axioma 4 (Simetria). (x,y) = (ya) Ve,y e V. Isso em conjunto com o 
axioma (3) implica na linearidade à esquerda. 


Axioma 5 (Não degenerado). Dado x e V, se (x,y) = 0 para todo ye V, 
então x = 0. 


Definição 2.5 (Matriz de Gram). Seja & = 
elementos do espaço vetorial V. Definimos [b; 


de É: 


fx, x2,...,Xn) um conjunto de 
5) = [Cx;,2:;)] de matriz de Gram 


Cx1,11) Lago) ... (11%) 
Abas) AM, Ma) re AB Ea) 


[xi 25)] = 
(Dart) REnsDa Susi MBas Dn) 
Teorema 2.3. Se £ = [11,192,...,2n) é uma base do espaço vetorial V, então 


o determinante da matriz de Gram de É é não nulo. 


Demonstração. Se det[(x;, w;)] fosse nulo, então existiria alguma dependência 
linear entre duas ou mais linhas ou colunas de [(x;, 1;)]. Suponha que a primeira 
linha seja uma combinação linear das outras, ou seja, existe (a2,a2,...,a") 0 
com aí e K tal que 


te ota) 


11 


Subtraindo (x4,a;) de ambos os lados e fazendo a! = —1, teremos 


Gu = O az; , 5) =0 


Mas isso é válido para qualquer x;, e consequentemente para qualquer elemento 
de V. Logo o elemento 5/a'x; é o elemento nulo de V devido ao axioma (5). 


Mas como (al!,a2,a3,...,a”) + 0, isso também significa que o conjunto E é 
linearmente dependente, o que contradiz a hipótese de £ ser uma base. E 
Denominamos 8 = (y,y>,...,Yyk) de subconjunto ortonormal de V se 


(Voy) = tô; VysyeB 


Teorema 2.4. Todo espaço vetorial V munido de produto interno possui uma 
base ortonormal. 
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Demonstração. Seja & = fx1,...,Ym) O conjunto com maior número de elemen- 
tos de V tais que 
Ci 2j) = Oi 


Suponha m escalares a” tais que 


m 
Se = 0 
F 


Aplicando o produto interno em ambos os lados em relação a x;, obteremos 
ts + 
3a (Li, 25) = 40,25) = 0 
i 


m 
0 = bs +a'ôi; = +aí 
i 


e consequentemente É é um conjunto linearmente independente. Logo a di- 
mensão do subespaço vetorial gerado por é, que denotaremos de M, é m. 

Sem = dmV = n, então € é uma base ortonormal de V. Sem < n, 
suponha K o conjunto de todos os elementos y de V que são ortogonais a todos 
os elementos de é, ou seja, 


(r,y)=0 Vreé VyekK 


O conjunto K também é um subespaço vetorial, pois a soma e multiplicação por 
escalar de dois elementos de K continua ortogonal a £, K inclui o elemento nulo. 
Se dim V — dim M = n—m é a dimensão do subespaço vetorial MO complemen- 
tara M, segue que MO = K (pois todos os elementos que não são combinações 
lineares de é acabam sendo ortogonais a É, uma vez que suas coordenadas nesses 
elementos, dadas pelo produto interno, serão nulas) e consequentemente 


dmkK > dmV-dmM=n—-m 


n <m + dimkK 


além de que a totalidade dos elementos que pertencem a M e K formam o 
espaço V. 

Se existir algum elemento de V que pertença tanto a M quanto a K, então ele 
é ortogonal a todos os elementos de & (por pertencer a K) assim como ortogonal 
a todos os elementos de K (por pertencer a M) e consequentemente é ortogonal 
a todos os elementos de V. Esse elemento só pode ser O devido ao axioma (5). 
Isso implica que M e K são mutualmente exclusivos (exceto pelo 0) e então 


n = m + dimkK 


Além disso, como O é o único elemento que é ortogonal a todos de K, então 
também podemos definir o produto interno dentro do subespaço K idêntico ao 
de V. 
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Agora seja ue K tal que (u,u) £ 0. A existência de u é garantida pelo fato 
de que se (u,u) = O para todo ue K, então 


Cutvutrv=Cuu+2lu,v)+Cv,v)=2Cu,v)=0 


para todo u,v e K. Como u e v também são ortogonais aos elementos de M, 
então são ortogonais a todos de V e consequentemente u = v = 0, o que viola o 
axioma (5) da não-degenerescência. 


Fazendo tm+1 = —, obteremos 


Su, u) 
(Li, Em+1) = +dim+1 Vi E (1,2,..,m,;m+ 1) 


uma vez que u é ortogonal a todos de é por pertencer a K e está normalizado. 
Então conseguimos construir um novo conjunto é = fx1,...,Ym, Um+1+ tal que 


Ci 25) = Oi 


que tem mais elementos do que £, que supomos ser o máximo. A única saída 
para esse absurdo é que m = n e então é é uma base ortonormal de V. E 


Teorema 2.5 (Produto interno como funcional). Seja a um funcional linear 
em V. Então existe um único ue V tal que 


a(v) = (u,v) VveV 


Demonstração. Seja € = (x1,...,1n) uma base ortonormal de V. Podemos 
escrever v em termos da base É: 


n 
v= õ v'g; 
F 
e então 


a(v)=a [3 a) - Brad a( 


KA 


Agora seja u um elemento de V cuja representação na base é seja u; = +a(x;), 


ou seja, 
= Con) o; (20) 


Segue disso que 


(u,v) = 6 ta(x;)t;, do = Do +via(x;) (ti; 2;) 


(u,vV) = »Divia a(x;)Oi; -3iva(ai = a(v) 


Como só há uma representação de u em uma dada base (ou equivalentemente 
só há um elemento u no qual as suas coordenadas na base É são u; = +a(x;)), 
segue que u é único. E 
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Teorema 2.6 (Caso particular da Lei da inércia de Sylvester). Dada uma base 
ortonormal qualquer € = (a1,...,2n) de um espaço vetorial V munido de um 
produto interno específico, temos que se 


Dimom)=r-s 


KA 


onde r é a quantidade de elementos x; de É tais que (x;,x;) = 1 es a quantidade 
de elementos xx de É tais que Cxp, xx) = —1, então para qualquer outra base 
ortonormal 8 = (y1,...,Yn) é válido 


DU) =1—8 


Denominamos o valort=r — s de assinatura do produto interno. 


Demonstração. Seja p o número de elementos y; de 8 tais que Cy, y)= lego 
número de elementos y; de 5 tais que <y;,y;) = —1. Como há n elementos da 
nas bases £ e 5, temos que 


r+s=p+qg=n (21) 


Suponha que r  p e que, sem perda de generalidade, r < p. Suponha também 
que a ordem das bases £É e 5 são tais que os positivos definidos apareçam 
primeiro, ou seja, Vi E f1L,...,r), <x,x;) = 1 e também Vi e f1,...,p), 
(Ui, Vi) = 

Seja f:V > K"*"—? uma função dada por 


F(z) = (C2,21), 2,220), LE, Dr), CE Yp4 DD 3-5 XD Un) VzeV 


Como r < p, temos n +r — p < n e consequentemente essa transformação (que 
é linear já que K"+"-? é um espaço vetorial e o produto interno é linear) passa 
de um espaço de dimensão n para outro de dimensão menor n+7r — p=9. À 
matriz F que representa essa transformação linear é de ordem y x n, com y < n. 
Então a equação matricial Fz = 0 tem soluções além da trivial”. Agora seja 
u € V um elemento não-nulo que satisfaça Fu = 0, ou seja, f(u) = (0,0,...,0). 
Isso é equivalente a 


Cutz)=0 Vi<xr e Luy=0 Vi: p<i<n 


Suponha que u possa ser escrito nas formas 


n n 

Dodns o mes 
u=),du;= >, Uy; 
j=1 j=1 
Para todo à < r, temos 
n n 
(u, Ti) = E q" (Lj, Ti) = = 005; Di, Ti) = 0º = O 

j=1 j=1 


TO posto da matriz F obrigatoriamente será menor do que o número de colunas, no que 
implica na nulidade de f maior que zero. 
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assim como b' = 0 Vi: p<i<mn. Podemos calcular (u,u> de duas formas, 
uma usando a base é e a outra usando a base 5: 


n 


CU, U) EA e 0º (x; 2h) a at0jn (Lj, 0) = PA 


“o, 


j=1 k 
a . 
= o (a)2(x;,2;) <0 
j=r41 
Por outro lado, temos 


n p 


(u,u) = At Pu) = (0 yu) > 0 


Logo temos <u,u) < 0 e <u,u) > 0, que é um absurdo. Portanto p = 
consequentemente da equação (21) temos q= se 


Smm=p-qg=r-—s 
% 


2.3.2 Produto entre p-vetores 


Agora sejam u,ve A”? V p-vetores. Iremos definir o produto interno no espaço 
A? V como sendo uma operação dada por 


GAVX AVE 


<u,v) = det[<u;, v;)] (22) 


onde 
u=UNUZA... A Up e V=V NAVAL. A Up 


com u;,v; EV. 


Suponha que v,w,ze A?V,wueV,a,be Ke que 


V=V AVIA.A Up 
f=(aw+bz)nuZA... Au 
PE MA A Tg 


h=2ZAUZA..AUp 
de forma que 


f=(aw)nruZa...nup+(bz) Aus A... Au, 


=a(wnus A... aup)+b(ZAUZA... A Up) 
= ag + bh 
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Segue então que 


Cv, (aw + bz) Lu,uz) ... <ur,up) 
Sua, (aw + bz)) Cuz,uz) ... Sua,up) 
(v, D = det 
Cup (aw + bz)) <upjuo) ... <UVpUp) 
alvi, wW) + o 2) Lui, uz) CUL, Up) 
rede a<vo, wW) E (u2, 2) a Wi tip) 
alvp,w) + bXvp,2) CVpyU2) CU Up) 
alvi,w) Cus,uo) +... <v,Up) 
RR sa: u) a ad Cuas tp) 
alVp,W) CUpuz) +... CVpyUp) 
bXvi,2) Cur,ua) (Ui, Up) 
bóva,2) Cva;uo) (U2, Up) 
det , 
don) (pt) co COpstp) 
Cu w) uu) ... <vi,Up) 
aa e u) o a Goa, tp) 
Cop W) CUVpua) +. CVpy Up) 
(1,2) CU u2) (vi, Up) 
Rir o 2) a uz) Goa, Up) 
Com) Cs) co (Ur) 


tu, |) E atv,9) + b<v,h) 


Mas como f = ag + bh, temos que 


(v,ag + bh) = aXv,9) + bLv,h) 


O mesmo é válido se substituirmos u; por (aw + bz), w e z em f, ge h respecti- 
vamente ou então introduzir uma soma de três ou mais elementos de V em u;. 


Portanto esse produto satisfaz o axioma (3) da linearidade. 
Além disso, o produto continua simétrico: 


<v,u> = det[<v;, u;>] = det[(u;, v;)] = det[Xu;, v;)] = <u, v) 
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uma vez que o determinante é invariante sob a transposta. 

Para verificar a não-degenerescência desse produto, suponha que € = fx1,...,%n) 
seja uma base ortonormal de V e o conjunto H seja uma das combinações de p 
índices h;, com h; < h;,1. Segue que todos os elementos xy para todos os H 
formam uma base de A? V. Sendo K uma outra combinação de p índices ks, 
temos 

(xH, TE) = det[ (zh, Xk;)] 
com TH =T A. ATh CXK STA. ATk SeHeK forem a mesma 
combinação de índices, teremos 


(Lhi Lk;) = (Uh; Yh;) = EOij 


o que significa que a matriz de Gram será diagonal com elementos +1 e portanto 
o determinante será +1. Isso significa que (xy, xx) = +1. 

E se não forem a mesma combinação, existirá pelo menos um índice q que 
pertencerá apenas a H ou a K. Logo, supondo q e H, temos (xa, tk,) = O 
para todos os índices k; e consequentemente haverá uma linha nula na matriz 
de Gramê. Segue disso que (xy, xK) = 0. 

Para quaisquer combinações H e K, temos o resultado mais geral 


(rm, TK) = LôH,K (23) 


Isso significa que os p-vetores xy formam uma base T ortonormal de A?V, 
característica herdada da ortonormalidade dos elementos x;. 

Finalmente, uma vez que garantimos a existência de uma base ortonormal 
de A? V, temos que se existir algum p-vetor u tal que seja ortogonal a todos os 
elementos de A? V e ele puder ser escrito na forma 


u=5 ul 
H 


então para todo elemento xx da base T temos 


(u,TK) = (Lutem : ox) = DEU eo = Djut(tômx) = +uf =0 
H H H 


o que implica u = 0 e isso prova a não-degenerescência do produto interno entre 
p-vetores definida em (22). 

Se € = (x1,...,4n+ for uma base ortonormal de V, então o elemento x = 
21 A... A gn forma uma base do espaço A” V. Portanto 


Cx, = Cx1,01): Cxo,09) 0... CLXn, Un) = (> 1)* = (—1)º 
Mas temos que 
s= (25) E >Oln 1) = (en — 21) = s(n+n 21) = s(n+ (r+s)-—2r) 
(n-(r-5)= 


1 
s==(n+tr+s—2r)= 
2 
8Também haverá pelo menos uma coluna nula, já que haverá um outro índice | que está 
em K mas não em H. 
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onde t=r— s é a assinatura do produto interno. Logo 


<a, a) = (—1)0-0/ (24) 


2.4 Operador estrela 
Sejaue A?Veve APV. A transformação 


VISUNV 


é linear e transforma um elemento do espaço A“ ?V para um elemento do 
A” V. Então se € = fx1,...,%7n+ é uma base ortonormal de V, u A v pode ser 
descrito como um escalar vezes um n-vetor x = x, A... A Zn normalizado, já 
que x sozinho já forma uma base do espaço A” V. Ou seja, 


unv= fu(u) 


onde f,: APV > K é um funcional baseado no elemento u. Sea,be K e 
we A”PYV, então 


fulav+ bw)z = uv A (av + bw) 
fulav+bw)z = a(u av) +b(u a w) = [afu(v) + bfu(w)] 


o que implica que o funcional f, é linear. Portanto o teorema (2.5) (aplicado 
ao espaço vetorial A”? V) garante que existe um único n — p vetor *u tal que 


Pulo) = Cru, o) 


e consequentemente 
unv= Cau, v)L (25) 


Como há apenas um elemento «ue A” ?V para cada u e A?V, podemos 
definir uma função, denominada de operador estrela (x), dada por 


E Re 
a(u)=(*u) => un v= Gu, va we À 


Naturalmente essa definição depende da base escolhida. Se y = —x, então y 
também está normalizado e também forma uma base de A” V. Consequente- 
mente 
unv=—Ceuv)y= (eu), v)y 

no que acarreta na mudança de sinal do operador estrela. Escolheremos uma 
base x = 27/A...A Zn específica que adotaremos como orientação positiva. Uma 
outra base y = y. A... Ayn de A” V terá orientação positiva se o determinante 
da matriz A que transforma a base x em y for positiva, ou seja, 


mA AY = Ara Arona ...s Ar, =|Al(r sn gx2 A... A Ln) 


50 


com |A| > 0. 
Podemos verificar que o operador estrela é linear, pois dado u,w e APV, 
abeKeve APV, temos 


(au+bw) A v=Ca(au+bw),v)r 
mas também 


(autbw)nv=(unav)+(wabv)= Ceu, av)ar+ Cew,bua 
= <a(xu),v>a + <b(sw),v>z = Ca(xu) + b(sw), va 


e portanto 
Ce(au + bw),vDaz = Ca(xu) + b(xw),v> x 


(Ce(au + bw),vD — Ca(xu) + b(xsw),v)) É = O 
Como «x é uma base, temos x 0 e 
Ce(au + bw),v) = Ca(xu) + b(xw),v) 
Ce(au + bw) — (a(xu) + b(xw)),v) = O 


Como isso é válido para todo ve A”? V, usamos o axioma (5) de não degene- 
rescência para concluir que 


«(au + bw) = a(xu) + b(*w) (26) 


2.4.1 Calculando o operador estrela 


Agora suponha u = 1x4 A 2x2 A... A xp, onde É = fx1,...,2n) é uma base 
ortonormal V com orientação positiva. Seja v = vx, onde K é uma combinação 
dos n — p índices de 1 a n. Então 


UNV=UATES (HU TR)L (27) 
mas também 
TA UNA e A A ES A A 


que pertence ao espaço A” V. Isso significa que se houver algum índice k; 
pertencendo ao conjunto (1,2,...,p), então haverá índice repetido e u A v = 0. 
O elemento u A v só pode ser diferente do nulo se K for exatamente a combinação 
fp+1,p+2,...,n). Por outro lado, isso também significa (pelo lado direito da 
equação (27)) que *u é ortogonal a todo xx exceto se K = (p+1,p+2,...,n). 
Isso significa que 


*U=CTk=C(Lo41ATp42 A. AL) 


onde c e K é algum escalar. Para determinar c, basta lembrar que x = 71 A 
Zo9 A... A Zn. Então 


unv=r=CGuzo)r=elar, rr = (tor 


o1 


Segue disso quec=1 se Xrr,zr)=1lec=-1se (xr,zxr)=-1. Ou seja, 
c= (Tr,ZK) 
e então 
su = (LE,TRK)LK 


Esse é um caso particular onde u = xp, A ...zhn, com h; = à. Se H é uma 
combinação qualquer de p índices, então K será o conjunto complementar de H 
em relação ao conjunto de índices 1 = [1,2,...,n), deforma que HU K = Te 
HnkK=% para queua vo). Logo 


*(2H) = (LR, ZK)LK (28) 


Sabendo como que o operador estrela age em um elemento da base de A? V, 
então também sabemos operar em qualquer p-vetor devido à linearidade (26). 

Além disso, devido à relação de comutação do produto exterior entre um 
p-vetor e um q-vetor, temos 


trata = (DA) Peg a tr) = (IPC (rg A ax) 
e consequentemente 
TKAXHS (—I1PC-Pg = (e(t2K),2H) 


Novamente é necessário que «(x x) = d-xy para algum de K para que zK nto £ 
O e consequentemente 


(—PM-Pg = d(zH, LH) 


Segue disso que 
(ni je = d(xy, 2H) 


Multiplicando ambos os membros por (x, «1H e sabendo que (Cxm, vH))? = 1, 
obtemos 
d= (=P (rm, ZH) 


no que implica 
*(210) = (PCP (rg, em) 2H (29) 
Note que há uma diferença entre aplicar o operador estrela a um p-vetor (equação 28) 
ou a um n-p vetor (equação 29). 
Agora note também que aplicando o operador estrela à equação (28) obte- 
remos 


*(+2n) = (gr, 21) (xx) = (IPO (rr, xr) (tm, 2H) ZH 
Como os vetores xn, são ortogonais, assim como x,, segue que as matrizes de 


Gram [Cxn,, 2n,;)] e [Cx 2%;)] são diagonais. Além disso, H e K são conjuntos 
complementares entre si e consequentemente a multiplicação de (xx, xx) com 
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(xr, xH) é totalmente equivalente ao determinante da matriz de Gram do n- 
vetor x =Í%71A...A Zn. Isso significa que 


(xr,ZK) Cx, ZH) = (1,9) (30) 


Mas da equação (24) temos (x,x) = (—1)("-9/2, onde t é a assinatura do 
produto interno e portanto 


*(*2H) = (PCP (0 DLrg = (IPC 
Devido à linearidade do operador estrela, segue que para qualquer p-vetor u é 
válido 
a (20) =: (= Pts 
Por último, seja ue A?V. Então 


ua (s2H)= Cu aza) 


O único elemento da base de p-vetores que gera u que não anula essa equação 
é xy e portanto 


gg A (exn)= tg A (Certo) vK) = Cer,vro(zg A TR) = CLR,TK)L 
Mas multiplicando (30) por (xy, xH)> em ambos os lados, obtemos 
(xr, TR) = (zu, tn) (2,0) = (um, e) (10-02 


e então 
ze A (s2g) = (DD rg,an) 


Agora se u,v forem dois p-vetores quaisquer de forma que 


! 
= atra e v=5 buy 
H H' 


ua (sv) = ba av (gm A *(xH')) = » ab (zm A *(xH))dHH 
HH! HH! 


ua (xv)= Da bE (1) (XE, ZH) 
H 


ua (uv) = (—1)M—0/2 p apr Cencen)| x 
H 


Por outro lado, temos 


tuv)= 5, atbE Lam, Tm) = o a! rm, 20)ônH! 


HH! HH” 
tu, v) = x» ad id (LH; TH) 
H 
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e portanto 

ua (ev) =(-D CD av) 
Note que trocar u com v e vice-versa não altera o resultado, já que o produto 
interno é simétrico. Consequentemente 


ua (a0)=0vA (au) (31) 


Exemplo 2.1 (Produto vetorial). Seja R$ o espaço euclidiano real com a base 
canônica ortonormal É = (61, €2, E3) que convencionamos como orientação po- 
sitiva e com produto interno definido por 


(Lp = ay + ay + ty? 
onde 


Ú 


Il 
o o di co 
8, 
D, 
& 
fas 
« 
Il 
Iv 
€ 
D 


Definimos o produto vetorial 
Ú x Es =*(T Ay) 


Segue disso que 


Ex = «((ely? - acyDé né) + *(elyê — vêyDé A 63)+ 
* ((x2y? — vy?)é» A €3) 


x y= (ly? - uy") «(64 &) + (ay = 1!) * (6 A 63)+ 


au 


(22yº — ay?) x (62 A &) 
Mas também 
«(& A €&2) = (63,63)63 = 63 
«(& A 63) = —(6,62)62 = —E 
«(6 n 6) = (GG, EE = & 
Note o sinal trocado em *(€& A €3) devido à orientação da base, pois sabendo 
que «(€j A €3) = yé2 com Yy = +1, então da equação (25) temos 
(61 n 63) An 6 = CG(& A 63),6>) (61 A 63 A 63) 
(En é A 6) = Gés,6>) (6 A 624 63) 
—(& n €& A 63) = 9(65,6>)(E1 A E2 A 63) 
—y(é n & A 63) = (€& A 624 63) 


Segue disso que Y precisa ser —1. E então 


Exemplo 2.2 (Produto escalar). Dado Y = qi + xo + «sk um elemento do 
Rº ewy o 2-vetor dado por 


Wy = nr) +y(kad)+ys(inj) 
com y = ni Er 2) + yak e Rº, então 
TA wy= (ray ay) é AjA k) 


Logo 
“(TA wy)= Ty +Hiyo+Hisys= Ty 


Por outro lado, temos que *(y) = wy, já que 


MSTNÊ 
«))=kni 
RR 


Então 


Note que de (31) temos 
(End) =+PaeB)=0:] 
o que implica na conhecida propriedade simétrica X-y=y-L. 
Exemplo 2.3 (Produto escalar triplo). Seja Rê o mesmo espaço euclidiano real 
do exemplo (2.1). Se 7 = >, 2'é;, então 
TAyAZ= Dixiyi (é AGA &)= Degnziyi (e A 62 A 63) 
ijk ijk 


Logo 


a(CAyAZ)= Demnniyi at «(& A 624 63) = > eipntyl 
ijk ijk 


que é exatamente o produto escalar triplo, definido pela equação (3). Logo 
Tlyx)=s(TAgAZ) 
2.5 Transformações lineares 


Vimos na seção (2.1) o caso particular das transformações lineares que levam 
um conjunto a ele próprio. Agora consideremos A uma transformação linear 
que leva um espaço vetorial M de dimensão m para um outro espaço vetorial N 


do 


de dimensão n. Dessa forma, dado x = x, A w2A...n xp, com x; E M, podemos 
construir uma nova transformação linear A? A dada por 


[Ã) (o): ÁM- ÀN 
[ÃA) (x) = Axy A Axo A... A Az, 


Denotamos A? A de p-ésima potência exterior de 4. A linearidade de 4 e a 
multilincaridade do produto exterior garantem que A? 4 também é linear e 
portanto é uma transformação linear. 

Considerando agora B uma transformação linear que leva o espaço vetorial 
V para M, temos que 4B leva V para N. E então, dado v = vy A V2 A ... A Up; 
com v; E V, a p-ésima potência exterior de AB fica 


DE anã 
[Ã as) (0) = [ÃA) REED 
[Ã as) (0) = [ÃA) (Ás) Area 
á)o-(ádhoe 


Como isso é válido para todo v e A?YV, segue que as duas transformações a 


seguir são equivalentes: 
p p p 


Dadosu=uAuUZA...AUpew=uwnwA.. aw,comue A?M 
ewe AM, temos que un we A?Í!M e portanto podemos definir uma 
transformação linear A??? A usando 4: M » N pela expressão 


p+q 
(AA) (un w) = Au, A Aug A... A Aup A Aw A Au A... A Aug 


Como o produto exterior é associativo, podemos reescrever a expressão acima 
como 


p+q 
[AA] (un w)= (Au A Aus A... A Aup) A (Aw A Awo A... A Awg) 
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3  Derivada exterior 


3.1 Espaço dual 


Definição 3.1 (Espaço dual). Seja V um espaço vetorial de dimensão n sobre 
um corpo K. Definimos o espaço dual V* como sendo o conjunto de todos os 
funcionais linearesV >» K. Sea,6BeV* ue V, ceK, a soma e multiplicação 
por escalar ficam definidos respectivamente por 


(a + 8)(u) = a(u) + B(u) 
(co)(u) = e(a(u)) 


Pelo fato dos funcionais serem lineares e da definição de soma e multiplicação 
por escalar, V* forma um espaço vetorial. 
Sejam a*'e V* ea;e K. Suponha que 


Duo =0 (33) 


onde 0 é o elemento nulo de V*, ou seja, é o funcional que leva qualquer elemento 
de V para o elemento nulo de K. O teorema (2.5) garante que cada funcional 
pode ser escrito em termos de produtos internos (desde que V esteja munido de 
m). Suponha que 
a" = (lu, > wevV 


Cada elemento u; pode ser descrito em função de uma base é = (x1,...,2n) de 
V: 
o “ 
Ui; = a bi Lj 
j=1 


de forma que 
Logo 


0= ns abilr;, )= E (Bout) os ) 


Isso significa que esse produto interno com qualquer elemento de V sempre será 
zero. Pelo axioma da não-degenerescência, temos 


Como é é uma base, 


o8 


Essa é uma equação matricial do tipo 
aB = 0 
com b! sendo o elemento da i-ésima linha e j-ésima coluna da matriz B e 
a= [ q... Gg | 


Essa equação matricial tem n equações e q incógnitas. Se q > n, obrigatoria- 
mente teremos soluções não triviais a O. Se q < n, só teremos a solução trivial 
a = 0 se o determinante de B for diferente de zero (ou equivalentemente se os 
elementos u; que definem os funcionais o” forem linearmente independentes). 
Isso implica, da definição de independência linear, que a dimensão do espaço 
dual V” é igual à dimensão do espaço V, pois q = n é o maior valor para o 
qual a equação (33) possa ser satisfeita apenas com a solução trivial a; = 0 e 
consequentemente fa!,...,a?) seja linearmente independente. 

Nesse caso, satisfeita a equação (33) com q = n, o conjunto fal,...,a 
forma uma base do espaço dual V*. 


Definição 3.2 (Base dual). Sejam € = fri,...,un) eT = fal,...,a”) bases 
de V e V* respectivamente. Denominamos P de base dual de E se 

o(z;)=6; Vij 
Teorema 3.1. Toda base de um espaço vetorial V munido de produto interno 
possui uma base dual em V* e é única. 


Demonstração. Suponha é = fx1,...,%n) uma base ortonormal de V. Se defi- 
nirmos os funcionais 
ot = (ti, ) 


segue imediatamente que o“(x;) = d; e portanto (a!,...,a”) forma uma base 
dual de €. Se fal, ..., ai) também for uma base dual de £, então 
(a* — ag)(x;) = 0(x;) — ab(x;) = 0; —- 0; =0 (34) 


para qualquer j. Segue disso que o funcional (a* — af) leva qualquer elemento 
de V a 0. Logo (a* — af) é o elemento nulo de V* e 


o = ab 
Como a equação (34) também é válida para qualquer 1, segue que a base dual é 
única. E 
3.2 Espaço cotangente 


Seja M um espaço topológico (como uma superfície ou um volume) de dimensão 
m. Definiremos um mapa de coordenadas y : U > R” com U sendo um 
subconjunto aberto de M e y uma função bijetora, contínua e com inversa 
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pl: p(U)CR" SU contínua. Cada ponto p que pertence a U pode ser 


levado univocamente a um vetor do espaço R” pelo mapa ( e vice-versa. 
Fixemos p um ponto em UC M. Sejam também y(t) e Y5(t) duas curvas 

em U definidas em um domínio [—1,1] tais que 44(0) = 95(0) = p. Definiremos 

“1 € Y2 como sendo curvas equivalentes se (p 0 7,h1)'(0) = (p 0 y5)'(0), ou seja, 


(aretnto)) | = (retrato) 


Em outras palavras, Yj e Y2 são equivalentes se o vetor tangente a M no ponto p 
pela curva 94 for igual ao de y, no mesmo ponto p. Definiremos [] como sendo o 
conjunto de todas as curvas equivalentes a y, também chamado de classe. Segue 
disso que [7] está associado a um vetor (po y)'(t = 0) que denominaremos de 
vetor tangente de M em p. Cada classe [7] fornece um vetor tangente distinto e 
a totalidade desses vetores formam um espaço vetorial T,M que denominaremos 
de espaço tangente. 


ER” (35) 


t=0 


Exemplo 3.1. Seja M uma esfera de raio r centrado na origem eU CM a 
superfície dessa esfera excluindo os polos norte e sul. Em coordenadas esféricas, 
o mapa py: UR? é dado por 


EN A 
º 


P(0,4) = rcos(4) sin(0)i + r sin(4) sin(0)j + r cos(0)k 


T 
Agora considere o ponto p de U em (o = 5 = 0) e duas curvas 91 e 92 que 


passem pelo ponto p em t = 0: 


Aplicando o mapa y nessas curvas obteremos 


P(n(t)) = r cos(0) sim(5 + Di + sin(0) sin(5 + Ji + reos(5 + Dk 


= rcos(t)i — r sin(b)k 


P(yo(t)) = r cos(t) sim(5 )ã + rsin(t) sin(5)j — roos(5 )k 


EN 
º 


= rcos(t)i + rsin(Dj 


e então d 
nl) = —rsin(t)i — rcos(Dk 


EN 
É 


Sep (t) = —r sin(Di + r cos(Dj 
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Em t = 0, obtemos dois vetores tangentes distintos de U em p: 


A 


(0 7Á1)(0) = —rk 
(p072)(0) = 1) 
que é o esperado, já que o plano tangente da esfera no ponto p (em que ié 


normal à superfície) é justamente o plano yz. Portanto uma base ortogonal de 
ToM é fr), —rk). Note também que 


(po n)(0) = rô (5.0) e (poy)(0) = r$ (5.0) 


onde 
(0, 4) = (cos(8) cos(d)i + (cos(0) sin($))j — sin(0)k 
99,6) = — sin(g)i + cos(d)j 


Isso não é uma coincidência, conforme veremos a seguir. 


Como estamos lidando com curvas equivalentes, iremos escolher a mais sim- 
ples de [y;]. Dado um sistema de coordenadas ortogonal, iremos fazer 


valt) = (pi + tô, po + tôi,...,Pm + tôim) 


no qual passa pelo ponto p = (p1,p2,...,Pm) quando t = 0. Esse foi o caso 
que usamos no exemplo (3.1), onde 44 e y> estão relacionados a uma mudança 
na coordenada q e O respectivamente. Denotando u; = p; + tô;; como sendo 
a coordenada generalizada u; parametrizada pela curva y; e xy a coordenada 
cartesiana no R”, temos que 


(port) = >, ze(u;(t)) e 
k=1 


Derivando em relação a t e usando a regra da cadeia, obteremos 


k=1j=1 Ou; dt 
n m 
ÔTk ca 
=D, 7 di 
k=1j=1 tj 
3 0% a 
=— Ck 
k=1 Ou; 
ep 
Mas da definição de versores em coordenadas generalizadas temos 
ro POR 0P 
Ú=— onde = 
o A; Ou ' Ou; 


Portanto, em t = O ficamos com 


e= A g= prio) (26) 


ou seja, os versores ú; no ponto p são os vetores tangentes (po y;)'(0) normali- 
zados. No caso do exemplo (3.1), fazendo új = ôeús =ó,temosM = A =r. 

Agora seja f: M > R uma função infinitamente diferenciável, ou seja, f E 
CY(M) (pequenos deslocamentos em M causa pequenas variações na imagem 
em R). A função f parametrizada pela curva %; é dada por 


(fot) = Hunt), volt)... tum(t)) 


Derivando em relação a t e usando a regra da cadeia, 


(Pen = om = Do daçto = da 
Agora vamos definir uma função 9, : TM —» R dada por 
Ds([n:]) = (fo w)'(0) 
Isso significa que 
dif) = 5] 1 (37) 


p 
ou seja, Dy([7:]) nos diz a variação da função f em p na direção relacionada à 
coordenada u;. 

Agora note das equações (36) e (37) que para cada curva y; parametrizada 
numa direção específica ú; no ponto p temos a aplicação de uma derivada parcial 


a . Nesse sentido, podemos denotar o vetor tangente como 
Ui 
p 


0 


a Ou; 


(9 0 w)'(0) 


p 


deixando implícito que esse elemento está definido em um mapa y sob uma 
parametrização [;] na coordenada u;. 

Como presumimos que as coordenadas u; são ortogonais, então os versores 
ú; o são e consequentemente, da equação (36), assim também são os vetores 
tangentes (poy;)'(0). Garantida a independência linear desses vetores, podemos 
formar uma base do espaço tangente T, M escrito na forma [or] das sos lp) 

Por outro lado, sabemos que &([;]) é um funcional linear, uma vez que a 
sua imagem é o corpo dos reais e 


Ô 
Oui 


padÃi 9 = ab;(bi)) + bold) 


p 


Daprio(bi)) = 5a] (af +09) =a 
p 
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para todo a,b e R. Logo os elementos b; satisfazem os requisitos da de- 
finição (3.1) e portanto formam um espaço dual de T;M, que denominaremos 
de espaço cotangente, denotado por T;“M. 

Os elementos de TÊM são denominados de vetores cotangentes (ou de cove- 


tores tangentes). 
)- 0 
p ou | 


Além disso, também temos que 
Segue disso que (9,1, 9,2,..., Dum) é à base dual de [| 


0 


Oui 


Du ([75]) — Pu ul = 6; (38) 


b Mais 
ainda, verificamos que 9; satisfaz a regra de Leibniz do produto de derivadas, 
pois 


ur | 
poco dum lp 


A 


og 
Ou” 


Byolbi)) = 55] U9) = 55) 9+ 


Dra(Ini)) = Sr) + Po(lu))f 
Isso nos motiva a reescrever 9; como sendo (df),, denotando a diferencial de f 
no ponto p. Portanto ((du!)p, (du?)p,..., (du )p) é a base dual de [5 E os TÊ o 


Como tudo isso é válido para qualquer ponto p em U, iremos omitir p nas 
notações dos vetores tangentes e cotangentes. 


j 


0 
Papi 
Uu 
É p 


3.3 Álgebra de p-formas 


Definição 3.3 (P-forma). Dado o espaço vetorial cotangente T*M, definire- 
mos os elementos do espaço AM(T*M) de p-formas. O espaço das O-formas 
AT*M) será o corpo CPM). 


As p-formas aparecem nas integrais, com a 1-forma aparecendo em inte- 
grais de linha, 2-formas em integrais de superfície e 3-formas em integrais vo- 
lumétricas: 


Ih Pa de roa db dn(o) dúdio ) He 2)deayas 


Porém essa relação será vista mais adiante. 
Sendo & = (da! ,da2,..., da”) uma base de A (TEM), então df, com f € 
Cº(M), pode ser escrito como uma combinação linear desses elementos: 


n 
df = > a; dy” 
i=1 
Lembrando que as 1-formas dx* são funcionais lineares, podemos aplicá-los a 


um vetor tangente -——: 
0x) 


(ap (5) - Doido (55) 
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, A A A ass 
Mas como (da'!, dx?, ..., da”) é a base dual de (5, 55,..., 3é= + (equação 38), 
temos 


ô Da 
(df) (55) = adj=a; 
1=1 


Mas da equação (37) sabemos que 


(d$) (55) = Oy([y;]) = SL 


Ra a of 
no que implica a; = ——. Portanto ficamos com 


Fon 


df = - da” 39 
f x aa 4º (39) 
que é a regra da cadeia para as diferenciais. 

Atentando à equação (39) podemos observar uma transformação d : 09 +» 
T*M tal que 


3 


a 
(0) 


d= - 
Ox 


i=1 


da' (40) 


ou equivalentemente d: A(T*EM) » A(T*M). Iremos generalizar essa ideia 
para uma transformação linear d que leva uma p-forma em uma (p + 1)-forma, 


p p+l 
do ADM) A(T*M) 
que satisfaz quatro axiomas: 
Axioma 6 (Linearidade). Vu,ve A(T*M), d(u+v) = du+ dv 
d(au) = adu VaeR 
Axioma 7 (Distribuição). Vue AM((TEM), Vve A(T*M) 
d(u nv) = (du A v) + (-DP(u a dv) 

Axioma 8 (Lema de Poincaré). Vue AM(T*M), d(du) = 0 


of 
Ox 


n 
Axioma 9 (Regra da cadeia). Vfe CP(M), df= > da” 
i=1 
Denominamos d de derivada exterior. 
O axioma 6 generaliza a linearidade para p-formas além das O-formas en- 
quanto que o axioma 9 foi herdado da equação (40). 
No axioma 7, colocamos um fator (—1)? no segundo termo para evitar in- 
consistências com o axioma 6. Para verificarmos isso, suponha ue A(T*M), 
ve A(T*M) e que o axioma seja modificado de forma que 


d(u nv) = dus v+a(u a dv) (41) 
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onde a é um escalar indeterminado. Trocando u e v de lugar, ganhamos um 
fator (— 1)?4, e então 


d(u nv) = d(-DPI(v À u)) = (—DPId(v A u) 
=(-1)(do sn u+b(v a du)) = (—1)PI(do A u) + (—1)PIb(v À du) 


onde b também é um escalar indeterminado. Sabendo que du é uma (p+1)-forma 
e que dv é uma (q + 1)-forma, então 


do a u=(-)AMD(u a do) = (—1)PHP(u À do) 


va du = (—1)P+*Da(du A 0) = (—1)PI*I(du A 0) 


e então 


d(u A v) = (—1)PI(—1)PI*P(u À dv) + (—1)PIb(—1)PI*I(du A v) 
= (—DP(u a dv) + b(—-D(du A v) 


Igualando essa expressão à equação (41), obtemos 
dunv+Ha(u a dv) = (-DP(u a dv) + b(—1)(du A v) 


[1 —b(=1)º|(du À v) + [a— (—1)º](u A dv) = 0 


Essa equação é válida para qualquer p-vetor u e q-vetor v, incluindo os casos 
em que du A v é linearmente independente de u A dv. Isso implica que 


1-b(-1)"=0 e a-(-1?=0 


Consequentemente a e b não podem ser quaisquer valores, eles são b = (—1)? e 
a=(—1)?, que é exatamente o axioma 7. 

Um exemplo em que (—1)? = —1 no axioma 7 é o divergente do produto 
vetorial (item 4 do teorema 1.8), dado por 


V(FxG)=G(VxF)-F(VxG) 


conforme será mostrado mais adiante (exemplo 3.5). 
Considerando x* e C(M), o axioma 8 implica que 


d(x! n dr? n dr a... n dr?) = d(al a (de? n dr A... n da?) 


d(x! n dr? À dr? A... a dz?) = 


dr! a (dr? À dei A... sn da?) + (—1) d(dr? À dê A... A da?) 
O último termo se anula para qualquer distribuição que fizermos e então 


d(x! n dr? na dr n...ndx?)=drl a dr ndria...n da? (42) 
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Segue disso que 
d(d(x! À da? A dai a... nda?) =0 


O axioma 8 tem como consequência as propriedades 


V(VxF)=0 


Vx(Vf)=0 
devido a aplicações sucessivas da derivada exterior, conforme veremos na seção 3.5. 


Teorema 3.2 (Existência e unicidade da derivada exterior). Existe apenas uma 
única transformação d que satisfaz simultaneamente os axiomas 6, T, 8 e 9. 


Demonstração. Suponha inicialmente a existência de uma transformação d que 
satisfaça os quatro axiomas e seja w uma p-forma dada por 


w = Say do! 
H 


com dr! = der na da2 A ...n date e H = (h,ho,...,hp) combinações de 
p índices. Iremos compactar os quatro axiomas numa só expressão aplicando 
todos eles em w: 


dw=d [ar uu!) 
H 
dy = sa (an da”) (Axioma 6) 
H 


dy = S'dag A da” (Axioma 7 e 8) 
H 


dy = x > gui dx' À dr! (Axioma 9) 


Essa transformação d então determina uma única (p + 1)-forma para cada p- 
forma w. Então se d existir, ele será único. Para provar sua existência, suponha 
agora que d é uma nova transformação definida estritamente pela expressão 


w = 5 am do => du = 5531 28 dx” À def (43) 
H 


H i=1 
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Então dados À = » by de” uma p-forma e p,q€ R, temos 


H H 


d(po + qA) = d a da! + qbg “) =d [tr +qbH) aut) 


d(w + A) = a +qdA 
o que mostra que d satisfaz o axioma 6. Além disso, agora considerando que À 
seja uma q-forma, temos 


WA Ã= (Bona!) A [tart) = Es agby def À def 
K 


H H,K 


n 


d(w A À) 5 sal (arbr) da A def à daé 
dx 
H,K 


Aus A) - DE Sera 
H,K 


o San PK ao A dr n daÉ 


H,Ki=l 


ue) 


ae 

Ob : 

“ drF à de” n de 
Egg 


di 


d(w a À) = [255 
H i=1 


K 


ia 


e 


H,Ki 


Ai A À) = da A A +(-1P [> am , NA DE 
H K 1=1 


d(w a) = doa ÀÃ+(-1D)P(w a dA) 


e portanto d também satisfaz o axioma 7. 
Aplicando d duas vezes em w, temos 


À 


2 
d(dw) = » | Dm da? À dai À daf 


Os índices são mudos, então trocando 1 por 3, obtemos 


(da) = DD DO joagas de a dei a dot 


Mas como ay € C(M), as derivadas parciais comutam e portanto 


E ar : E 
As as apro A A 
1=19= 


O único elemento de um espaço vetorial que é simétrico dele próprio é o elemento 
nulo e consequentemente 
d(dw) = 0 


Logo d satisfaz o axioma 8. E finalmente, o axioma 9 é satisfeito de imediato. 
Isso prova a existência de uma transformação d que satisfaz os quatro axiomas. 
Usando a primeira parte da demonstração, concluímos que a derivada exterior 
d é única. E 


Exemplo 3.2 (Gradiente). Ao aplicar a derivada exterior a uma O-forma f 
infinitamente diferenciável no espaço R$, obtemos 


ô 2 ô 
df= Fido + A 


dy + dz 


As coordenadas de df na base fdx, dy, dz) são justamente as coordenadas do 
gradiente de f em coordenadas cartesianas. 
Agora definamos uma transformação linear DP : A (TER) > Rê tal que 


T(da)=1 e I(dy)=j e I(d)=k 
Note que essa transformação leva uma base de ATER?) a uma base de Rê e 
ambos os espaços vetoriais possuem a mesma dimensão, o que implica que T é 
bijetora e então existe [1 :R3 A!(T*R?) tal que 

Dlii=de e ps =dy e Dk = dz 
Usando essa transformação em df, obtemos 

Pdf) = V$ 

Exemplo 3.3 (Rotacional). Aplicando d em uma 1-forma w dada por 


w= P(x,y,2) de + Q(x,9,2) dy + R(x,y,2) dz 
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com P,Q,Re Cº(Rº), obteremos 


dy = d(P de + Q dy + Rdz) 
do = d(P da) + d(Q dy) + d(R dz) 
dv =dPA de+dQ A dy+dRa dz 


Foh o 
EC qu eo a a A dy+ 
Õ Oy z 
E aa q — dz | A dz 
dx Oy 


Os termos com dx A dx, dy A dy e dz A dz são zero e portanto 


P P 
dw = E nr dada E dinda 
Oy z Fopi 
a RN NR RO CRE 
0z fopi Oy 


Uma vez que definimos dx A dy A dz como uma base de orientação positiva 

Ê a 3 AS Rd qe 
(sistema dezxtrogiro) de A“(T*Rº), então iremos colocar em evidência as 2- 
formas com orientação positiva: 


dw = El do dy a dz + aa dz a da + EE or dx a dy 
dy Oz dz Ox dx Oy 


Identificamos que as coordenadas de dw na base (dy A dz, dz A dx, dx A dy) são 
justamente as do rotacional em coordenadas cartesianas de um campo vetorial 


F(r) = P(x,y, Ji + Q(x,y, O) + R(a,y, 2)k 
Como (dx) = dy A dz e idem para dy e dz, temos que 
«(dv) =(V x F),de+(V x F), dy+(V x F), dz 
Usando a transformação linear T do exemplo (3.2), podemos escrever w como 
=P) 
e também 
Ped) =(VxPit(VxP+(VxF)k=VxF 
Finalmente temos que 


VxF=T[:d(T!(F))] 
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Exemplo 3.4 (Divergente e Laplaciano). Dado uma 2-forma w 
w=P(x,y,2z)dya dz + Q(x,y,2)dz A de + R(x,y,2) de A dy 
com P,Q,Re Cº(Rº), temos que 


do =d(Pdya dz)+ d(Qdzn de) + d(Rdz A dy) 
dy = dP A dy Aa dz+ dQdz An de+ dRda a dy 


dw = (So dr + a Se de A dy A dz+ 
Õz 


a e ap a dz A da+ 

Fok Oy õz 
dit CE Ri o A dx A dy 
Fohi Oy õz 


Excluindo os termos com diferenciais repetidos e sabendo que 
dy n dzn de= dz A de A dy= de A dy A dz 


concluímos que 


— (oP 00 OR 
a (+ 


A coordenada de dy na base (dz A dy A dz! é justamente o divergente em 
coordenadas cartesianas de um campo vetorial 


Jar a dy sa: 


F=P(ryoi+Q(a,y + R(x,y,2)k 


Usando a mesma transformação linear T dos exemplos passados, devemos 
notar que 
w = À 


com 
A=Pdr+Qdy+ Rdz=T- HF) 
e também notando que 
oP O oR 


Ad o 


V-F 


segue que 


V.F=+d(«!(F)) 


Já para o laplaciano, que é o divergente do gradiente, temos que se f é uma 
função real escalar, então 


Vºt=V.Vf=V.T(df) =+d(«0HT(df))) 


V2f =+d(*df) 


TO 


Resumindo, temos que se f é uma função real escalar e F é um campo vetorial 
em R$, então podemos definir o gradiente, divergente, rotacional e laplaciano 
em coordenadas cartesianas pelas equações 


Vf=T(df) (44) 
V.F=+d(«["!(F)) (45) 
VxF=T[=d(T-!(F))] (46) 
V2f = «d(+df) (47) 


Exemplo 3.5. Se A ew forem uma 1-forma e 2-forma respectivamente de um 
espaço 15 M de dimensão 3, então vimos anteriormente que calcular dA é equi- 
valente a calcular o rotacional do campo vetorial associado a À e e dy corres- 
ponde ao divergente do campo associado a w. Também vimos no exemplo (2.1) 
que o produto vetorial de dois campos está associado ao produto exterior entre 
eles (produto exterior entre dois I-vetores). Então seja y uma 1-forma e 


wWw=AÃAY 
Segue do axioma (7) da distribuição que 
do = (day) — (AA dy) 


Associando À com F ey com G, então w está associado a F x G, da V x F 
edyaVxG. 

Mas também vimos no exemplo (2.2) que o produto escalar está associado 
ao produto exterior de um 2-vetor com um 1-vetor. Isso significa que dA A Y 
está associado a (VW x F)-G enquanto que À A dy está associado a(V x G)-F. 
E daí surge a propriedade 


V(FxG)=G(VxF)-F(VxG) 


3.4 Mapeamento e mudança de coordenadas 


Sejam U e V subconjuntos de R” e R” respectivamente. Um mapa do relaciona 
coordenadas em U para coordenadas em V, 


do:Uo Vo dolx)=y 

PEÃO quiçã”) (48) 
com y' e CºU para cada coordenada. Dado também f:V > R uma função 
infinitamente diferenciável, podemos compor do e f em uma única função qg f : 
UsR 

do! = fodo 

de forma que dj f agora tem domínio em U e leva para R. Ou equivalente- 
mente &y leva uma função de domínio V para uma função com domínio em 


a 


U. Conforme vimos na seção (3.3), a função escalar f configura como uma 0- 
forma pertencente a AME: *V) enquanto que dg f é uma 0-forma pertencente a 
A“(T*U). Então podemos afirmar que 


0 


0 
g8: A(TV) o A(TU) 


Analogamente, podemos definir uma outra função 97 que leva uma 1-forma 
em V para uma 1-forma em U apenas realizando a substituição (48) e usando 
a regra da cadeia (axioma 9): 


di dy' = - da? (49) 


Podemos ver que 4% é linear, pois dy! + dy? = d(y! + y?) e as derivadas parciais 
são lineares. Então podemos afirmar que 47 configura como uma transformação 
linear. 

Para uma 1-forma qualquer w em V, 


w= E ay", uy) dy' = >. ai(y) dy” 
; i=1 


to = 55 aiy(o))55 do? 


ox) 


Essa é a aplicação do mapeamento para 1-formas. Para generalizarmos q* para 
p-formas, usaremos a propriedade (32) das transformações lineares, 


(AA) (u A qw) = [(Ã A) | á [A] (1) (50) 


Por exemplo, para o caso das 2-formas fazemos p = q = 1 na expressão acima e 
teremos 


o3 (dy! n dy?) = di(dy") A ói(dy”) 
m oy! , m oy? ; 
* 1 2 — - i 2 dy 
ps (dy A dy”) p ari 1) A É 217 ea 


m 1 2; 
dildy! a dy?) = 5) ELE TE (51) 
1, 


= dx! Oui 


Como à e j são índices mudos, podemos trocá-los de lugar: 


Al 2 m Anal 2 

OU OU ca ; Oy Oy 

a es -— da J p) 
2800) aqi AT A da A dad dqi A da (52) 


ó3 (dy! n dy”) = 
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Somando dz (dy! A dy?) com ele próprio com as representações (51) e (52) ob- 
temos 


* 1 2 * 1 AN ue! t e + J 
do (dy A dy”) +65 (dy A dyê) = EM EB dx' a dz A 25 dm dx' À dz 


263 (dy! À dy”) = 5, 


ij=1 


(E oy Oy! Oy? 


() j 
Oxº dxi Ox) E) Adndo 


e por fim 


2 ei Ox' dai dxi dg? 


1 A [oyo oyloy : 
dE(dyt À dy?) = e ( Sa 404 agi n dei 
O termo entre parênteses é justamente o jacobiano de y! e y? em relação a x! e 
x), dado pela equação (12). Logo podemos reescrever a expressão anterior para 
obtermos 


E RUC a | 
* 1 QN , 1 ) 
gs (dy A dy”) 5 2 dai, 7) da" À dx (53) 


Teorema 3.3. Dados 6; : A(T*V) — AP(T*U), w ey p-formas e À uma 
g-forma, então 


1. dolw + 7) = GE(w) + d5(7) 

2. Goralwv A À) = do(w) A di(A) 

3. d(ójw) = Gp (do) 

4 Seg:U»Vey:VW, então (pod) = dio 


Demonstração. O item 1 para p = O é trivial. Sabemos que 97 é linear. Para 
verificarmos se 43 é linear, 


( 
(dy") A di(dy? + dy”) 
(dy') À (gi(dy”) + di (dy?) 
(dy 

( 


dy! À dy?) + 63 (dy! À dy”) 


e portanto 43 também é linear. Suponha que dj é linear para um certo p e que 
w seja uma p-forma, então 


dpi (tw A GR = dp lv) A di (dy? +!) 


Podemos escrever w como a soma de duas p-formas. Devido à linearidade do 
produto exterior e de dj, segue que 47, , também é linear e portanto provamos 
por indução o item 1. 
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O item 2 advém diretamente da propriedade das transformações lineares (50), 
conforme visto anteriormente. 
Para o terceiro item, vemos que se f é uma 0-forma então 


Aplicando dj e usando sua linearidade, obtemos 


exam = 5) 5, ON am 5 (5 UOOD E a 
Mas da regra da cadeia temos 


3 Of(y(x)) dy" Oflyle))  Aóo(f)) 
dy Oui Ox Oxi 


i=1 


Logo 


m 


ót(as) = a CD aus = a(o8(1)) 


Pela linearidade, essa propriedade é válida para quaisquer 0-formas f e portanto 
o par de funções dj e dj satisfaz o item 3. Agora suponha que essa propriedade 
seja válida para algum p. Usando o item 2, temos 


pelos dyPO) = dj (w) n di (dy?!) 
Aplicando a derivada exterior em ambos os lados, obteremos 


Aópy(iw n dy”) = d(gi(w) A di(dyP*!)) 
Aóka(o A dy?) = d(óf(w)) A SF(dy”HT) + (1 (ói(0o) À A(óf (dy) 


Uma vez que 


dói (dyP*)) = 43 (d(dyP*")) = 630) = 
e por hipótese 
d(dj (w)) = Apra(do) 
então 
dp (w A dy?!) — pai(do) A di(dyP*") 
Lópr(w n dy?t!)) = dP a(du À dyPt!) 


Mas também 
d(w a dyP*!) = dy a dy?! + (DP (w a d(dy?* 1) = dy A dy? tl 


e portanto 
Mópr(w n dyP*!)) = dp a(d(w À dyPt)) 


T4 


o que prova por indução o item 3 de 4% para qualquer p. 

Dada uma 0-forma g que possui domínio em W, gg é uma O-forma que 
possui domínio em V. Por sua vez, d% (15 9) é uma 0-forma com domínio em U. 
Então 


[614 0 09] GO) = g(y(ó)) = [409] (660) = do | [bo 9] 69) = [(60 0 40)9] (x) 


Para p = 1, verificamos que se z denota as coordenadas em W, então 


(Bro Wt) dz = qt [Wi dz] = dt 1-6 [5 av 


(6toyf)de= 55) : | 


Pe A DN 

i=1j=1 a) ottay 

(droyd = o om AU aut 
i=1j=1k=1 


Retraindo a soma com o índice à com a regra da cadeia, obtemos 


(ot o wpde = 51 5) SLGD E aut 


j=lk= 
m m Fa) 
(ot oude =D A jd 
dE 
* * o Ôz y(x)) 
(ot o vt) do — 3 ESA au 
(PT o Wi) dz = (Wo di dz 


e portanto o item 4 é válido para p = 1. 
Agora suponha que o item 4 seja válido para algum p e que w seja uma 
p-forma em T*W. Então 


[(40 Ajralio n de] (x) = [lb 0 SP (w)] (x) A [(1y 0 d)i(dz?*")] (0) 
= [(97 0 45)(w)] (0) 4 [(G7 0 bi) (da? *1)] (0) 
= [(p41 0 Wpra)*(to A da? *)] (x) 
o que completa por indução a demonstração. E 
O item 4 do teorema (3.3) diz que mapear de W diretamente para U é equivalente 
a mapear de W para V e depois deV para U. 


Exemplo 3.6. Sejam d er constantes reais positivas e & o mapeamento Rº >» 
Rº dado por 


aq = dcos(p) + r cos(8) cos() 
y = dsin(y) + r cos(0) sin(p) 


2 =rsin(0) 


To 


Figura 1: Toro comd=2er=1. 


ou seja, v = x(0,9), y = y(0,9) ez = 2(0,9). Considerando py o ângulo da 
projeção de um ponto q no plano xy em relação ao eixo x (ângulo azimutal) e 
9 o ângulo entre q e o plano xy (a latitude) no sentido de fora para dentro da 
superfície, então todos os valores possíveis de py e O formam uma superfície no 
R3 que é denominada de toro (figura 1). Podemos verificar isso de forma mais 
nítida se escrevermos 


x cos(p) —sin(y) O d+rcos(0) 
y [=| sin(y) cos(p) O |- 0 (54) 
z 0 0 1 rsin(0) 


que é a matriz de rotação em torno do eixo z aplicada a uma circunferência de 
raio r no plano xz cujo centro está no eixo x mas deslocado em d. 
Dada uma 1-forma w = (x? + y?) dz, temos que 


diw = [[d+ rcos(0)]? cos2(49)) + [d + r cos(0)]? sin?(9))] . [5 do + — dy 


djw = [d + rcos(0)] «r-cos(0) do 


No caso particular em que d = 0, o toro se transforma em uma esfera (que 
é uma circunferência girando em torno dela própria na equação 54) e então 
dtw = r3 cos?(0) do. 


T6 


Agora se w é uma 2-forma dada por w = dz A dy, então 


d2w = (di da) A (di dy) 


00 Õ 08 
x dx dy Odyda 
= 0 
Pau E EE da Rá 


djw = | — rdsin(6) cos?(y) — r? cos(0) sin(0) cos? ((p) — 
rdsin(0) sin2() — r2 cos(0) sin(0) sin?(() | dó À dy 
dSw = [(d + rcos(0))r sin(0)| dp A do 


3.5 Inversa do lema de Poincaré 


Usando as equações (44) e (46), podemos ver que 


Vx(V$)=Vx(Ddf)) == a(D(D(df))] 
=T[*d(df)| = D(x0) 
=0 
ou seja, o lema de Poincaré d(df) = O implica que o rotacional do gradiente 
(definidos pelas equações 44 e 46) de uma função real suave f é sempre o vetor 


nulo. 
Usando as equações (45) e (46), também temos que 


Vo(VxF)=V-T[=d(D-!(F))] 
= «d(«DMT[*d(T!(F))]) 
= + d(*(+d(T!(F)))) 
= + d(d(D!(F)) 
=0 


onde usamos o fato de que 
“(DT (E) = (IP-P+O- Dra (p) 


comp=1,n=3et=3, com t sendo a assinatura do produto interno (que é o 
produto escalar usual, no qual é positivo definido). Então o lema de Poincaré 
também implica que o divergente do rotacional (definidos pelas equações 45 
e 46) de um campo vetorial real suave é sempre zero. 

Por outro lado, vimos na seção 1.2.2 que sob certas circunstâncias, se Vx F = 
O então F = VV para alguma função escalar Ve quese V-F = 0, então 
F=Vx A paraalgum campo vetorial A. Em outras palavras, nessas mesmas 
circunstâncias, se dv = 0, então há algum a tal que da = w, de forma que 
d(da) = 0. Iremos verificar no próximo teorema quais são os requisitos para 
isso ser verdade. 


qq 


Teorema 3.4 (Inversa do lema de Poincaré). Se w é uma p-forma (p > 1) em 
um conjunto estrelado U E R” e dy = 0, então existe uma (p — 1)-forma a tal 
que 

da =w 


Demonstração. Seja 1 o intervalo dos reais [0,1]. Podemos montar um objeto 
(t,u) que está associado simultaneamente a um te T e a um elemento ue U, 
de forma que a totalidade desses objetos formam um conjunto 1 x U. 

Dada uma função escalar bem comportada f = f(ax',t), onde xº são as 
coordenadas relacionadas a U, temos que a 1-forma 


df = 3 Lar pela 


Og? ot 


i=1 


pertence ao espaço A (T*(I x U)). 
Agora seja jj o mapa que leva ue U para (1,u), dado por 


Mm Us (1 x U) 
jiu) = (1,u) 

De maneira similar, definimos jo: 
Jo Us (1 x U) 
jo(u) = (0,u) 


Conforme vimos na seção 3.4, podemos definir funções ii, e Jo, que levam uma 
p-forma com domínio em 1 x U para uma p-forma com domínio em U: 


E: NTUxU)o A(TU) 
E: NT xU) o A(TU) 


Os mapeamentos 71 e jo são dados por 


Então no caso em que p = 0, por exemplo, temos que 


diets) = fio, 1) 
Jd 00',8) Re F(a*,0) 
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enquanto que com p= 1, 


is of * i of * 
= p ari Ji, (da ] + DE Ji, (dt) 
i=1 t=1 t=1 
= gi E ap da 
a Ox? q) ot A Ox) 
ice qu t=1 t=1j=1 


ot : 
Como t independe das coordenadas, ao O. Além disso, uma vez que as 
x 
ami 


Es Sa E : Ux ; 
coordenadas x” são linearmente independentes entre si, temos - = dj. Logo 


Ox) 


si (df) = Dnles dy” 


De maneira totalmente análoga, 


jé (df) = E: dy 


Em geral, dada uma p-forma w, It, w faz a 1-forma dt desaparecer e fixar 
t=1emw. Já jo,w também limpa os termos com dt mas faz t = 0. Isso é 
garantido pela propriedade distributiva dos mapas (item 2 do teorema 3.3), já 
que qualquer termo com dt resultará em um jy (dt) = 0 ou jo, (dt) = 0. 


Agora iremos definir uma transformação linear K : AMT*(I x U)) » 
AP(T*U) da seguinte forma: K sempre levará ao elemento nulo qualquer argu- 
mento que não envolva dt, isto é, se dr” e APR U)eageR, então 


K (an(t,u) de”) = 0 


E quando há dt no argumento, se dr” e A(T*U) ea; E R, definimos K pela 
equação 


K (as(t,u)dt a dr”) = [/ as(t,u) at dx” (55) 


Iremos primeiramente provar que 
K(dw) + A(K(w)) = 1, 4 — Jo, (56) 


Para isso, suponha À = ay(t,u)dt A da” uma (p + 1)-forma em 1 x U que 
contenha dt. Então 
fas da as 
dA = + >2 dtadta de” 
se E | a adta de 


i=1 
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=5, (-5e) dt À dx À de” 
= Ox 


Aplicando K a dA, obteremos 
[au a da? 


KO) = |V as(t,u) at da” 


podemos aplicar a derivada exterior a esse resultado para obter 


a? 
o Fá osttou)dt au? A da” 


i=1 


d(K()) = = FE o at dr' À da” 


i1=1 


DITA das 


Agora sabendo que 


que é K(dA) mas com sinal trocado. Então 
K(dN) + d(K(A)) = 0 


Mas também sabemos que À é um termo só e ele possui dt, o que leva IL À =0 
€ J0,41A = 0. Esse resultado mostra concordância com a equação (56). 

Agora vamos para o caso da (p + 1)-forma w = ag(t,u) dr” que não tem 
dt. Então 


du = [Dita | 4 Ear ant 
i=1 


Logo K(dw) fará o primeiro termo com a somatória desaparecer, sobrando ape- 


nas 
=| =[f. o Vara 1! = [ay(1,u) — an(0,u)]do” 


Por outro lado, K(w) = 0 e portanto d(K (w)) = 0. E também sabemos que 


an(Lu)dr? = tw e an(0,u)dr? = gw 


e portanto 
K(dw) + A(K(w)) = 1, 4 — do, (57) 
Como K é uma transformação linear, essa equação se estende para qualquer 
(p+1)-forma em 1 x U, seja ela dependente de dt ou não. 
Agora seja 9 um mapa dado por 


puderes 
d(1,v) = v 
(O, v) = vo 
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onde vo = (vl,v2,...,v”) aponta para um ponto no conjunto U no qual U é 
estrelado. Segue que & é dado por 


dote, tt) (uliaê, seu) (58) 
onde u' = u'(x?,t) são coordenadas tais que satisfazem as condições de contorno 
u(x),1) = 2º e u' (13,0) = vº 


i i 
a existam para todo 


além de que u' e U e que as derivadas parciais 


te [0,1], Yi,j e em todo o domínio U. Exemplos de bobidenada desse tipo são 
u =v+it(a' — ow!) e ul = vt + (x! — vu?) sin(5t) 
Dado que o conjunto é estrelado em vo, então usaremos 
u = v' + (1º — v)h(t) 


dh 
onde h(0) = 0, h(1) = 1 e que — exista Vte [0,1]. Então 


dt 
Ou 5 
= 0 -h(0 
Dada uma função escalar g(u!,u2,...,u”), segue de (58) que di (9) leva a 
O-forma g às coordenadas x? e t: 
o0(9) = g(u'(a?,t)) (59) 
E para uma 1-forma dg dada por 
d = : KA 
g x Zu du 
então 


dg du" |, Og Ou! 
ido) = 2, E das as Pan 


= h() > DO des x Gta; (60) 


Aplicando jf e jo, à 1-forma óf(dg) em (60), iremos obter 
i=1 


Buds 
-da'=d 
Ox á 9 


Ji. (ot(a9)) = (1), E 


t=1 


jê,(ot(A9)) = (0) D) E 


Então para qualquer p-forma w e CPU, com p > 1 temos 
(1, Pp)u =w (61) 
Co, º Pp)w =0 (62) 


Isso pode ser verificado por indução. Se as equações (61) e (62) são válidas para 
um certo p, então dada uma 1-forma À e uma p-forma w temos 


aa (Ppailio n A) = ia (Pp (tw) A PT()) 
= (dp (w)) A jiló (A) 
=WA 


Da mesma forma, 


opa (Palio A AD) = do, (Pp (0) A do, (GTM) = 0 


Agora seja w a p-forma do enunciado, ou seja, sabe-se que dy = 0. Então 
usando 47, 1w como argumento na equação (57) obteremos 


K(d(ó%)) + A(K (Go) 
K(óf(dio)) + A(K(óo)) 
K(6%0) + A(K(óo)) 
A(K(ójw)) = «o 


A (ópw) — do, (Ppvo) 
dO Sp)u — (do, º Pp Jw 


( 
w— O 


Fazendo a = K(djw), concluímos que 
da =w 
o que completa a demonstração. E 


Note que se duas p-formas 8 e a produzirem uma mesma (p+1)-forma w 
pelo operador d, então 


d8 = da =w d(8- a) =0 
e portanto da inversa do lema de Poincaré temos que 
Borda 


ou equivalentemente 
B=a+dy 


Isso significa que 8 e a podem diferir de apenas um dy, com y sendo uma 
(p-1)-forma arbitrária. 
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Exemplo 3.7 (Divergente nulo e o rotacional). Suponha w uma 2-forma dada 
por 
w=P(x,y,2z)dyn dz + Q(x,y,2)dz A de + R(x,y,2) de A dy 
de forma que em um subconjunto U de Rº estrelado na origem é verdade 
oP O oR 
dw = ( + a 


Ee 4 SE) den dynde=0 


ou equivalentemente 


E a 
dv Oy Oz 
Podemos construir um mapa & : ([0,1] x U) +» U dado por 


dl, 9,2; t) = (tz, ty, tz) 


cuja imagem é a reta que conecta a origem até (x,y,2). Pelo fato do conjunto 
ser estrelado na origem, toda a reta pertence a U. 


Usando a equação (53), temos que &3 leva uma 2-forma em U para uma 
2-forma em [0,1] x U dada por 


4 
do (Pdyn dz) = Pts tinto) a) dz” À da? 
ij=1 


com q! =t. Os únicos jacobianos não nulos (além das comutações de xº e xi) 
serão 
O(ty, tz) a t2 
Oy, 2) 
OCtystz) 
tz) 
o(ty, tz) 
Ay, t) 
O(ty, tz) 
> = yz 
o) 


= tz 


Porém o primeiro resulta em um fator de dy A dz enquanto que o último resulta 
em dt a dt = 0. Ao aplicarmos K a d5(Pdy A dz), o primeiro desaparecerá, 
e então os únicos relevantes serão o segundo e o terceiro. Com as comutações 
de x' e x), aparecerão termos iguais porém negativos, o que pode ser mudado 
ao comutar dx” com dx), resultando em um fator 2 que cancelará com 1/2, 
resultando 


dE(P dy a dz) = Pltx, ty, tz) -[ytdt n dz +tzdy A dt+t? dy A dz] 


Para deixarmos tudo com orientação positiva (primeiro dt e depois as coordena- 
das x, y ez) e então usar a relação do operador K dada por (55), comutaremos 
dy com dt. Fazendo isso e colocando t em evidência, ficamos com 


dS(Pdyn da)=t-Plte,ty,tz)- [ydt À dz — zdt A dy] + (termos sem dt) 
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De maneira análoga, encontraremos 


Go(Qdza dr) =t- Q(te,ty, tz) - |zdt A da — a dt A dz] + (termos sem dt) 
d(Rdx a dy) =t- Rlte,ty,tz)- [xdtn dy — ydt A dx] + (termos sem dt) 


Aplicando K a G5w, obteremos uma 1-forma dada por 


K(giiw) = | | “E Plte tuto) at (dep dd 
|V t- Q(te, ty, t2) at (2 de — « do)+ 
[/ E) at dis 


Para calcularmos dK(G3w), vamos aplicar a derivada exterior para cada termo: 


1 3 1 
Õ 
d E t-Plte,ty,tz at de) = — | t-Plte,ty,tz at dx" À dz 
(|/ (ta, ty, tz) dt | y Em 5! (ta, ty, tz) 


1 a 1 
= Ed | yt: P(tx, ty, tz) at dradz+ + | yt: P(tx,ty,tz) at dyndz 
ox (Jo Oy Lo 


1 1 
-|/ ut: So de [de a de + | [ EP + utço de | dy de 
[o dx 0 oy 


Mas 
oP oP ta) oP 
de tr) de (ta) 
e também =t E e então 
fa) o(ty) 


1 1 
| q at [do de + | [ ED at| dy n ds 
0 Ata) 0 ty) 


De maneira análoga, temos 


a([[ trai) av) E 


: op : ar 
Dr ARE . 2 
| zt E) at) do n ay + | t.P+rzt (82) at[ de n dy 
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fimose 


1 1 
| so de] dy a dae + | [ 1-Q+ ate A dede a do 


a([/ tai] var) - 


1 A ij A 
2. 0Q 2 00 
| E sao) | dy a de + | [ t-Q+at 209 dt] de a de 


a([[ + nai) av) E 


: oR : oR 
teu dt d . 2 dtld 
| z 09 deJden v+|[ te R+at 65 “| x A dy 


' OR : oR 
E ão dt de ndo + | [ t.R+ ea da A dz 
| 7) j “ at) |O 


E então juntando os termos de dK(G5w) que estão relacionados a dz A dy 
ficamos com 


dK(djw) = ...+ 
1 
oP oQ oR oR 
+ at. A +. R+yo— +et = at[de nd 
/ ôlta) ty) É da) "É ata) E 
Mas da equação (63) sabemos que 


oR oP oQ 
tz) Ota) ty) 


é verdade para todo t em [0,1] (pois o conjunto é estrelado) e portanto 


dK(djw) = ...+ 
i oR oR oR 
Pepe at Rea = db 
+| z att) + R+gyt Fc a(tx) | x A dy 
Agora note que 
dR om OR ta) A OR (ty) s OR o(tz) 
Ota) Ot ty) Ot o(tz) Ot 
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Logo 


E nd 
0 


Note também que 


dp 2 

sa =P 4% 

E [tR] +2t- R 

e portanto 
Rudi 
dK(djw) =...+ | = [ÉR] at dz a dy 

1 

=... +tER(te,ty,tz)| de A dy 
0 


=... + R(x,y,2) de A dy 


O raciocínio é inteiramente análogo para mostrar que isso também acontece com 
os termos relacionados a dy A dz e dz A dz, o que implica 


dK(d5w) = P(x,y,2) dy A dz + Q(x,y,2) dz A de + R(x,y, 2) de A dy 
dK(éfio) = w 


Então encontramos uma 1-forma a = K(pw) tal que da = w. Agora note que 
mostramos no exemplo (3.3) que calcular da, com a sendo uma Í-forma é equi- 
valente a calcular o rotacional de algum outro campo vetorial. Se reescrevermos 
a na forma 


= | to uma dz+ 
| er- zP) at dy + [/ tur -sqjal] Ee 


então identificamos o espaço vetorial A cujo rotacional reproduz o campo veto- 


rial (PQ, R): 


az ( K H2Q — yR) dt, | tlaR — zP) dt, f t(yP — xQ) at) 


SeF =(P,Q,R), então observamos que 
F(t-r)x(t-r)=H2Q-yRi+taeR-—zP)j+tyP — cQ)k 
e portanto 


A=[ F(t.r)x (t-r)dt 


que é exatamente a expressão que usamos na demonstração do teorema (1.11). 


86 


3.6 Coordenadas generalizadas 


Sejam &(x,y,2), E2(x,y,2) e Es(x,y,2) um conjunto ortonormal de campos 
vetoriais que variam suavemente em uma região do Rê. Eles podem ser escritos 
na forma 


é; = fala, y;, ai + fila, Yy, aj + fista, Yy, 2)k 
Segue disso que 
é&-&-fafa+ fel; + fisf;s — ij 


Aplicando a derivada exterior nessa expressão, obtemos 


d(é; - é) = (df) f; + (dfiz) f;z + (dfis) f;s + faldf;a) + finldf;o) + fis(df;s) 
d(é; - É;) = dé; . é; = é; . dé; 

= dô;; 

=0 


Portanto temos 


Cada vetor dé; pode ser descrito usando a mesma base (&, &, E3!, 
dé; = wn& + w;péa + wi;3€3 (65) 

mas com os elementos w;; sendo 1-formas. Segue disso que 

dé; Ê é; = Wij 
Por outro lado, trocando à e j de lugar obtemos 

dé; ê é; = Wii 
Usando a equação (64), chegamos à conclusão que 

Wij + Wji = 0 


ou equivalentemente 
Wji = —Wij (66) 


Isso significa que podemos montar uma matriz anti-simétrica Q com as 1-formas 
Uma variação infinitesimal no espaço Rº em coordenadas cartesianas pode 
ser descrita com um vetor dr cujas componentes são 1-formas, dado por 
dr = dai + dyj + dzk 
Essa mesma variação pode ser descrita usando a base ortonormal (&, €2, Ea): 


dr = cj€&4 + os€2 + 0363 (67) 


com o; sendo 1-formas. 
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Fazendo 
e-| & oc=[ o 09 03 | e Q= [wi] 


e considerando que a derivada exterior d aplicada a uma matriz é a matriz com 
d aplicado aos seus elementos podemos compactar as equações (65), (66) e (67) 
com as expressões 


de = Ne (68) 
Ff =-9 (69) 
dr = ce (70) 


Uma das utilidades dessa notação é a de que pode ser facilmente generalizada 
para dimensões maiores do que três. 

Como os campos é variam suavemente no espaço, as derivadas parciais 
comutam e portanto como consequência vale o lema de Poincaré: 


d(déj) =0 e d(dr)=0 
Da primeira igualdade, de (65) obtemos 
d(dé;) = (dun)éi +... + (dw;a)é3 + (1)! (wi déy +... + wi3 d63) 
Usando (65) na segunda parcela, teremos 
d(dé;) = (dw;n)&i +... + (dw;s)é3 — [ww + wizwo + Wigw31] &— 
[wiw12 + Wigtwao + wisw32] Ea — [wiw13 + Wizwo3 + wiswaa] és 


Note que os termos entre parênteses são os elementos da multiplicação de 2 com 
ela própria. Segue disso que 


0 = (dO — Q2)e 
dO = 2 


De maneira similar, temos 


Exemplo 3.8 (Coordenadas esféricas). As transformações das coordenadas 
esféricas para cartesianas são dadas por 


x = rsin(0) cos(y) 


y = rsin(0) sin(p) 
z =rcos(6) 
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Aplicando o operador diferencial às três coordenadas obtemos 


da = sin(0) cos(y) dr + r cos(0) cos() dO — r sin(0) sin() dy 
dy = sin(0) sin() dr + r cos(0) sin(p) do + r sin(0) cos(p) dy 
dz = cos(0) dr — r sin(0) dô 


e portanto o elemento infinitesimal dr fica 


dr = (sin(0) cos(y), sin(0) sin(y), cos(6)) dr+ 
(r cos(9) cos(), 7 cos(9) sin(p), —r sin(0)) do+ 
(— rsin(0) sin(y), rsin(9) cos(y), 0) dy 


Agora sabendo que 


mo 
Il 


(sin(0) cos(g), sin(0) sin(p), cos(9)) 
(cos(6) cos(y), cos(0) sin(7), — sin(0)) 
(— sin(p), cos(y), 0) 


€& Db 


podemos reescrever dr na forma 
dr = dr? + r d9ô + rsin(6) dp 


Daí identificamos 


e = 


3» 
0) 
[14] 

Il 
“D» 
fa) 
0) 

[00] 

Il 

GS 


o; = dr os=rdo e os=rsin(0)dy 


Agora note que 
01h 02403 =r2sin(0)dr A do A dy 


é o elemento de volume dx A dy A dz. 
Para encontrarmos a matriz Q, primeiramente aplicamos a derivada exterior 
ao versor £: 


df = (cos(8) cos(2) dO — sin(0) sin() dp)i+ 
(cos(0) sin(p) dO + sin(0) cos(7) dy)j — sin(0) dok 


que é ! 
dê = 0? + d90 + sin(0) dog 


Portanto identificamos 
wi = O wo =dô e wgz=sin(0)dy 
De maneira similar, para dô e dQ encontramos 


dô = (— sin(0) cos() dê — cos(6) sin(p) dp)i+ 
(— sin(0) sin() dO + cos(4) cos(7) dp)j — cos(9) dok 
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d$ = — cos(p) dçi — sin() dçj + 0k (71) 


ou equivalentemente 
dô = — dor + 08 + cos(6) do 


Para encontrarmos dQ na base (£,0,9!, devemos escrevê-lo em função de ? e O 
de forma que a coordenada em k se anule. Para isso, multiplicamos £ por sin(6) 
e O por cos(8) para então somar: 


sin(0)f + cos(0)ô = (sin2(0) cos(9), sin?(0) sin(y), sin(0) cos(0)) + 
(cos”(0) cos(7), cos?(0) sin(p), — sin(0) cos(0)) 


sin(0)? + cos(0)ô = ((sin?(0) + cos?(0)) cos(7), (sin?(0) + cos?(0)) sin(), 0) 
= (cos(y), sin(y), 0) 


que é (71) a menos de um fator — dy e portanto 


A 


dy = — sin(0) dp? — cos(0) dp6 


Com isso, a matriz Q = [w;;] fica 


0 dô sin(0) dy 
Q = — dô 0 cos(8) dy 
—sin(0) dp — cos(0) dy 0 


Note que Q é anti-simétrica, conforme vimos anteriormente. 
Cada versor da base ortonormal fé, é», €3+ pode ser escrito na forma 
é = bai+ b;o)+ bigk 


Na notação matricial, isso é 
e= Bi 


onde agora B = [b;;] e 


br) Ce) mo> 


e:-e =[é.6]=1 


onde [é; - €;] é a matriz de Gram dos versores é, com a multiplicação dos 
elementos sendo um produto interno. Mas por outro lado, temos 


e-e" = (Bi) (1 B”) = B(-1/)B” = BIB” = BB” 
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e portanto B é uma matriz ortogonal, ou seja, BBT =I —s Bl = B7. 
Da equação (70) temos que 


dr=ce=oBi 


Entretanto também temos 


dr=dri+dyj+dzk=| dz dy dz |- 


Dr) Ce) mo) 


dr =| dx dy dz Ji 


Segue disso que 
[ de dy dz Ji=oBi 


[dr dy dz |=oB 
ou equivalentemente 
3 3 3 
dz = x bio; dy = > bo;0 e dz = b3k0k 
i=1 j=1 k=1 
Logo 


3 3 3 
dz À dy Adz= n no) A [> hos) A [> is) 
j=1 k=1 


i=1 


3 
dz À dy Adz= > byibo;bak (o; A OG A Ok) 
ijk=1 


Conforme já vimos na seção (2.1), essa última expressão é equivalente a Boy A 
Bos A Bos, que por sua vez é equivalente a |B|(o1 A 02 A 03), onde |[B| é o 
determinante de B. Mas como B é uma matriz ortogonal, então 


1=detI = det(BB”!) = (det B) - (det BT) = (det BJ? 


e portanto det B = +1. Assumindo que 0; A o; A ox forma uma base com ori- 
entação positiva, conforme vimos na seção (2.4) por definição uma base positiva 
se transforma em uma outra base positiva usando uma matriz com determinante 
positivo. Logo det B = 1, implicando que 


den dyadz=01A02A 03 (72) 
O resultado (72) aparece no exemplo (3.8), onde encontramos 
dz A dy A dz = r2sin(0) dr A do A dy 
Dado e = Bi, temos também i = B-le e portanto 


de = (dBi = (dB)B te 
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Mas de (68) sabemos que de = Ne e consequentemente obtemos a matriz anti- 
simétrica 

= (dB)B 
Esse resultado pode ser estendido pelo seguinte teorema: 


Teorema 3.5. Se 4 for uma matriz ortogonal cujos elementos são funções das 
coordenadas então (dA)A”! será uma matriz anti-simétrica. 


Demonstração. Dado que A é ortogonal, então A” A = T. Aplicando a derivada 
exterior em ambos os lados obteremos” 


d(AAT) = di 
(dA)AT + A(dAT) = 0 
(dA)JAT! + A(dA)T = 0 
(dAJAT! + ((dAJAT)" =0 


((d AJA)” = (a Aja 
o que mostra que a matriz (LA)A”! é anti-simétrica. E 


Teorema 3.6. Seja 4 for uma matriz cujos elementos são funções das coor- 
denadas definidas em um conjunto U e À é ortogonal para um certo ponto uo 
em U. Então se dA = DA, onde D é uma matriz anti-simétrica, então A é 
ortogonal em todo U. 


Demonstração. Seja C = ATA. Então 


dC = d(AT A) 
= (dADJA + AT(dA) 
= (DAJTA+ AT(DA) 
= ATDTA+ ATDA 
= —ATDA+ ATDA 
=0 


e portanto a matriz C = ATA é constante em todo U. Mas como ATA = T 
para um certo ponto uy em U, então segue que para qualquer ponto em U vale 
ATA = T e consequentemente A é ortogonal em todo U. E 


9A derivada exterior se distribui num produto de matrizes semelhante ao produto de 
funções, já que se (AB);; = Dj uikdas é o elemento da i-ésima linha e j-ésima coluna de 
k 
AB, então 


A(AB) ij = >, Maskbrj) = >, A(ask)bry + Dik d(byp;) = [(dA)B o A(dB) a 
k k E 


ij 
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Dada uma transformação u = Au,, com A ortogonal, temos que uy = 4" !u 
e então a variação infinitesimal em u é dada por 


du = d(Aug) = (dA)uo = (dA)A tu 
À nova posição do vetor u devido a uma variação infinitesimal du é portanto 


u+du=u+(dAjA"!u = (1+(dAJADu 


3.7 Gradiente, Divergente, Rotacional e Laplaciano 


De agora em diante denotaremos du A dv apenas como du dv, ficando implícito 
o produto exterior entre as 1-formas e compactando a notação. 

Seja fu,v,wy um conjunto de coordenadas do R$. O vetor posição r pode 
ser descrito como 

r= a(u,v, wi + y(u,v,w) + (u,v, w)k 
Uma variação infinitesimal em uma dessas coordenadas leva r a outro ponto 
próximo. No caso da coordenada u, o vetor de deslocamento infinitesimal é 
dado por 
Or Ou, Oya Ozs 
e a) 


du ane Sm Sê 


de um or . 
Iremos definir €& como sendo o versor de Fu? ou seja, o vetor que mostra a 
u 

direção para onde r muda quando u tem uma pequena variação. Então mul- 
tiplicando essa derivada parcial por alguma função À de u, v e w obteremos 
&: 

10r 1 0% lo . lo0zs 
= ias j+ k (73) 

Adu Adu Adu A Ou 
Analogamente para v e w temos 


10r ldx, 1oy, 1Toza 
+ 


A 


ei 


= = k T4 
sa uu Ou nov. E Ou) uu Ou SE 
10 1 - 1lOdy, 1 E 
ge lo a (75) 
vôw vÔw vOw v Ow 


para certas funções u e v que respectivamente tornam é€> e é3 normais. 
Agora note que aplicando a derivada exterior a r obteremos 
ox ox dx é Oz Õz Oz R 
ar = ( du + — dv + du it + ( du + — dv + = du )R 


Ou Ou Ow Ou Ou Ow 


e portanto 


or 
= — — ué e — = vé 
ou Ou pe w a 
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e então 
dr = (Aduj& + (udv)é2 + (v dw)és 


Comparando essa expressão com dr = 0161 + 0963 + 0363, concluímos que 
o1 = Adu oo=udvo e og=vdw 


As coordenadas (fu, v,w) são ditas ortogonais se (&, é>, €&3) formarem uma 
base ortonormal do R$. Supondo que c10403 constitui uma base positiva do 
espaço das 3-formas, então segue que 010203 = dz dy dz e então poderemos usar 
o operador estrela nas 1-formas o; da maneira usual que vimos na seção (2.4). 


Exemplo 3.9 (Gradiente). Dada uma função escalar f e um conjunto de co- 
ordenadas fu, v,w), temos que 


E qua Lapa of 


9) 
Ou Ou Ow ei 


df 


e então 10f 10f 10f 
fo) 
Xu?! ii u dv? ia v ow? 


Agora iremos definir!º uma transformação linear T tal que 


df 


Segue disso que 


Iremos definir essa expressão como sendo o gradiente de f: 


fofo gilóros 4 af 


1 es 
A du pu Ou v Ow 


Vi=Idi)= 


No caso particular das coordenadas esféricas, u=7r,v=0,UW LP, À 1, 
u=rev=rsin(0), no que leva a 

of 10f, 1 Of 
=Dt+ 0 + Ô 

or r 00 rsin(0) op? 


Vf 


Exemplo 3.10 (Divergente). Seja F um campo vetorial dado por 
F = FHé& + Fé + Fes 
Aplicando T-!, obtemos 


r!(F) = Fo, + Fo509 + F303 


10A transformação T está completamente determinada, já que du, dv e dw são linearmente 
independentes (e consequentemente as 1-formas o; também são). 
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Agora aplicando o operador estrela: 


*D1(F) = F,0903 + F50301 + F30109 
= (uvF) dv dw + (AvFo) dw du + (AuFs) du do 


Calculando a derivada exterior dessa 2-forma e sabendo que apenas uma das 
derivadas parciais deve sobreviver em cada termo, ficamos com 


ô ô 
(uv Fi) du dv dw + dy (ALE) dv dw du + ER! 


E + CC ovb) + O a) du dv dw 
Odu Ou Ow 


d(+1"!(F)) = AuFs) dw du do 


2'| ob) 


1 9 9 Õ 
=— [> (uvFy) + — (Avi) + — (AuE3) | o109203 
Auv À du dv ow 


Aplicando novamente o operador estrela, obtemos 


” 1 0 0 0 
*d(«D1(F)) = E (aluno) + dy (ALE) + 5 ura)) 


Iremos definir essa expressão como sendo o divergente de F': 


V.F=+d(+[!(F)) = xo (uv) + É (ul) + lui) 


Em coordenadas esféricas, ficamos com 


E l O 2a Epa 0 
V.F= rIsin(0) (tr sin(6) F,.) + 20!" sin(0)Fg) + E) 
ST do e da. 1 0F, 
= Rg ER rsin(0) ag (sin(8) Fo) ii rsin(0) dy 


Exemplo 3.11 (Rotacional). Dado o mesmo campo vetorial F do exemplo 
anterior (3.10), temos que 


P-1(F) = Fio, + Fp02 + F303 
= AF, du + ufo do + vFs dw 


Aplicando a derivada exterior, obtemos 


d(T!(F)) = CE) dv du + CAES dw du + Eiras du dv+ 
Ou Ow ou 
EM dw dv + E du dw + in) dv dw 
w ou Ov 
ATE) = ( E ERR q )) dv dw + (Lar) É um )) du du+ 
do TE Gus Gu Pe duo 
0 0 
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Voltando para as 1-formas oi, 


AT) = (Aus) (uEo) oros+ 55 (Ei) A (UE) Jogo 


uv Ou 
1/0 9 
a (pu (UF) = (AF) )o1o2 
Usando o operador estrela, temos que 
1/0 Õ 1,0 Õ 
—1 o 
Has) uv E (Fs) ow (no) Jos a Av E (AS) du (23) )oa+ 
Fo) fa] 
xa (pa (43) 2 (AF) os 
Aplicando T nessa expressão, chegamos a 
near) = (> qr)-L qum))ari (Ar) É (um) 
uv NOv ow Avi ow ou 
1,0 Õ 
Au E (ufa) dv (NF) ) 6s 


Iremos definir esse resultado como sendo o rotacional de F': 


VxF=T(d(T(F)) 


A 


VxP= o (ler “(us )S + (DOM) É (vFy) 6o+ 


ow vkOw Ou 
1,0 Õ a 
Au E (ufa) Ov (NF) ) 6s 
Em coordenadas esféricas, o rotacional é dado por 
1 Õ : Õ E 
VxF= Tanto (597 sin(0),) — ap(rFo) + 
1 0oF, O : A 0 OF 
r sin(0) ( dp Or (r sin(0)F)) ii E (rs) 00 ) 
o 1 0. Fo 
Vxr= to) (sono) — 55)! 


a 1 0F, (9) 


. 1/0 
sin(0) dp Or (rF%)) Ha r - (ro) 20 


E 
Exemplo 3.12 (Laplaciano). Seja f uma função escalar das coordenadas u, v 
ew. Vimos no exemplo (3.9) que 


Dior lo doi: 


df Ron o du v Ow 


3 
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Aplicando o operador estrela, obtemos 


10 10 10 
«df = ERA + EO a + REP RA 
Vou 


A Ou Hu Ou 
Avôf Au Of 
9 


= — — dudu + — — dwdu + — — dudu 
u H Ou v Ow 


Agora usando a derivada exterior, apenas uma derivada parcial deve sobreviver 
em cada termo, resultando 


didni= E (25) du dv dw + = (=) dv dw du + L (5) dw du dv 


10 fvuôdf É 0 [AvoOf E 0 (Au Of dado 

du À À du dv Nu Ou dw À v Ow 
1/0 /vuoOf a dO Av Of 4 o (Mu Of Acid 
Auv | Ou À du dv À u Odu Ow À v Ow aire 
Usando novamente o operador estrela, chegamos à expressão 

1/0 /vuodf 0 [Avof o (Mu of 
Di, d =: 
e, ala (2) + (E) a (E 


no qual definiremos como laplaciano de f, denotado por V?2 f. 
Em coordenadas esféricas, o laplaciano é dado por 


= a fã my, O frsin(0)0f | O aa 
ess =Dalé (1 Sino, e) "2. ( r 2 “do (= | 


10 /5,0f 1 0(. of 1 of 
= + sin (0 + 
12 Or (r or) sasnço) 00 9 | anita) di 
Em suma, dada uma base ortonormal com orientação positiva (&, é, E) 
com coordenadas fu, v, w+ e de forma que 


dr = cj&j + 0262 + 0363 
então podemos definir uma transformação linear T' tal que 
T(o;) = é 
e portanto, para uma função escalar f e um campo vetorial F' dado por 


F=hHé+beé-+ Pés 


segue que 
Vf=INdf) (76) 
V.F=+d(+["!(F)) (77) 
VxF=T[=d(T-!(F))] (78) 
V2f = «d(+df) (79) 
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Se o = Adu, 03 = udv, 03 = vdw, então podemos reescrever essas equações 
para 


10 lo 1 
=. as Lg + OS e 


Adu Hu Ou R De (e) 


9 Fa) 9 
V.F=—— (oa luvro) + dy MALES) + Ra) (81) 


110 /vuoOf 0 [Avof o (Mu Of 
D fi 
sea 5/à (8) + (E * dy v Ow AE) 
Uma vez que T leva o espaço das 1-formas na base (01,02,03) para o espaço 


R3, se quisermos encontrar a matriz que representa [' na base fi, jk) basta 
lembrarmos que é; = bai + b;oj + bisk e portanto 


? ? di bi 
Po epa do |=| by 
E a: Ôis big 


Segue disso que a matriz desconhecida que representa T é 


ba boa ba 
bio boo bao | =B7 
bis bos bas 


que é ortogonal. 

As equações do gradiente (76), divergente (77) e laplaciano (79) são auto- 
maticamente generalizáveis para n dimensões. Para o gradiente, df retornará 
uma combinação linear de n 1-formas linearmente independentes, no qual T(df) 
será um vetor do R”. Para o divergente, «[”!(F) será uma combinação li- 
near de (,",) (n-1)-formas, no que d(*["!(F)) será uma n-forma e portanto 
* d(*P-!(F)) será um escalar. Idem para o laplaciano. 

Já para o rotacional d(P”!(F)) deve retornar uma combinação de linear de 
(5) 2-formas linearmente independentes, no que *d([”!(F)) será uma com- 
binação de (eo) (n-2)-formas. Como só definimos T para agir em 1-formas o;, 
T em geral não estará definido para agir em (n-2)-formas, necessitando genera- 
lização. Se definirmos T para agir usando a equação 


DP(o;jo;...Om)= 6 AGA ...A em 


então V x F será um (n-2)-vetor. 
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Naturalmente há outras operações disponíveis, como por exemplo D(x d(x d(T1F))) 
retornará um campo vetorial do R” se Fe R?. Já sen = 2, então «d(T-1F) 
será um escalar: 


r!(F) = Foi + F505 
= AF, du + ufa du 


Ou Ou 


d(T!F) = | : (no) & (ar du dv 


| (ur) = Sar] 0102 


e portanto 


«ATE = | À (um) (48) 


Em coordenadas polares, podemos ver que o, = dr, 02 = rdô. Logo fazendo 
u=rev=t6entãoobteremosAÃ=1,u=re 


0 0F, | 


+*d(T/F) = - [5 (r Fo) 


Essa expressão é exatamente (V x F), em coordenadas cilíndricas no R3. 


3.8 Superfícies 


Ágora iremos fixar a coordenada correspondente ao versor €3, no que resulta 
em um subdomínio do R? dada por uma superfície, na qual designaremos de XD. 
Sabendo que dr está dentro do plano tangente gerado por € e €2, então dr não 
tem componente na direção é e consequentemente 03 = 0. De dr = ce temos 


dr = c1€4 + 0962 (84) 


Uma vez que as coordenadas são ortogonais, segue que é3 será normal a toda a 
superfície 5X. Conforme r se move em 5, o versor é acompanhará a superfície. 
O conjunto dos versores normais à superfície » forma uma imagem com um for- 
mato esférico, denominado de imagem esférica de X. Por exemplo, na superfície 
cilíndrica paralela ao eixo z o versor À gerará uma circunferência no plano xy. 

Sabendo que a matriz (), relacionada à variação dos eixos é em X é anti- 
simétrica, 

QT = 0 


então podemos escrever essa matriz na forma 


0 W3 —w1 
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Da equação de = e vamos obter 


d& 0 W3 —wi & 
de> = —W3 O) —WwW9a €&2 
des W1 Wa O) €s 


que nos fornece as equações 


dé = w3€2 = w1 Es (85) 
des = —ws& = wa €3 (86) 
des = w& + wa€2 (87) 


Já da relação dO = 2, temos 


0 dws Es dw1 O) W3 —w1 0 W3 —w1 
Ea dws 0 E dws = —W3 O) —wWwga º —Ww3 0 —Ww2 
dw1 dws O) W1 Wa 0 wi Wa O) 
—lWswgs — WIjw1 —WIitwg2 —tWwgwga 
= = Waw1 —Wsws — Wawa W3w1 
—Waw3 Witws —WIiwW1i — Wawa 
0) —Wwiwo - —twgswga 
= —tWwow1 O) Wgsw1 
—tWwows WwW3 0) 
e portanto 
dw = —Wwows (88) 
dws = WIiws (89) 
dws = —Wywa (90) 
E quanto a do = 00, 
O) W3 —wi 
| do; dos do |=[ a 09 0]. —wW3 O) —tWw9 
W1 Wa 0 
= [ —Ogws Ojwsz —Ojwi — Ogwa | 
= [ W3092 —wW301 Wj01 + w209 | 
e então 
doi = W309 (91) 
dos = —W301 (92) 
O = w/o + wo09 (93) 


As equações (84) a (93) fornecem várias informações sobre a superfície D, tais 
como a sua curvatura, elemento de área, entre outras. 
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Uma vez que 0102 é uma base das 2-formas em 5, podemos escrever ww» 
em termos dessa base: 
Wiwo = Kojos (94) 


A função K é denominada de curvatura gaussiana de um determinado ponto 
na superfície 5. Essa quantidade está relacionada à curvatura intrínseca de XD, 
ou seja, referente às medidas feitas dentro e ao longo dessa superfície, sem levar 
em conta o espaço onde X reside (que é o espaço R3). Por exemplo, quando 
K > 0 numa região da superfície, podemos nos deparar com uma circunferência 
de raio R dentro dessa região que pode ter uma circunferência e área diferentes 
de 27R e TR? respectivamente quando medidos dentro da própria superfície. 
Já a 2-forma oywa — oowy também pode ser escrita na base oj03, 


Oiwo — oow1 = 2Ho105 (95) 


onde a função H é denominado de curvatura média de um ponto em D. Ela 
está relacionada à curvatura extrínseca, isto é, associada a como X está imersa 
no espaço onde reside. 

Analogamente, dado que (01,02) é uma base das 1-formas em X, podemos 
escrever w1 e w» como combinações lineares de 0, e 09: 


Wi = poi + q02 


Wo = TO + soa 
onde p, q, r e s são funções em 5. Mas devido à equação (93), 
O = wj04 + wo09 
temos que 


0 = (po + g05)01 + (roy + sos)o» 
=(—q+r)oo» 


Como 0105 é uma base das 2-formas em 5, o,09 * O e portanto r = q. Logo 


wWj = poi + qo» (96) 


Wo = q01 + 809 (97) 
Substituindo essas equações em (94), vamos obter 
Kojos = (por + go5)(goy + sos) 
= (ps — q)o10» 


e consequentemente 
K=ps-q (98) 


Agora substituindo w, e w> em (95), temos 


2H0105 = or(goy + sos) — oa(por + q09) 


(s + pjo10oa 


101 


e então 
2H =p+s (99) 


Agora note que as equações (96) e (97) podem ser escritas na forma matricial 


Cas 


Os autovalores da matriz 2 x 2 satisfazem 


Re É [-0 (100) 
go s-u 


na qual chegamos na equação do segundo grau 
2 2 — 
ut —(p+sjutps-—qê=0 


As fórmulas de Viete afirmam que as duas soluções uy e [to dessa equação 


satisfazem ) 
—(p + s 
Ha + Ho RES fino cade =p+s 


ps— 


E 
Comparando esses resultados com (98) e (99), concluímos que 


=ps-q 


Wi" Ho = 


K=muo e 2H=mn+yu 


Os autovalores 1 e |t2 são denominados de curvaturas principais de DX. 
Para encontrarmos ws3, usamos uma das equações (91) e (92), 


doy = wgos ou dos=-ws0o1 
E da equação (90), temos 
dws = —Wiwo =— —Kojos 


Isso significa que basta sabermos o, e 02 que já temos condições de encontrar a 
curvatura gaussiana K da superfície D. 


Exemplo 3.13. Seja a curva no plano xy dada por 


v(9) = cos(0) + In(tan(0/2)) 


onde n/2 < 0 <7. A revolução dessa curva em torno do eixo x gera uma 
superfície 5, na qual pode ser vista na figura 2. 
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as à 248 ; 


0.75, 09 -1.00" 


Figura 2: Superfície 3 no intervalo 7/2 < 0 < 7. 


A parametrização dessa revolução pode ser dada pelo ângulo py em relação 
ao eixo y e introduz uma nova coordenada z. Então z = ysin(y) = sin(0) sin() 
e deve aparecer um fator cos(p) em y A superfície X, portanto, é dada por 


x(0, 7) = cos(0) + In(tan(9/2)) 


ondem/2<0<me0<y<27. Conforme vimos na seção anterior, podemos 
encontrar cy e 09 referentes a O e q usando as equações 


oj=Adô e ocy=udy 


onde 
e a = E 
[0 pre õy 
r=(0,0)i+y(0, 0) + 2(0, p)k 
Portanto 
Ss = ( sin(0) + : sec?(0/2) cor(o/2)) à + cos(p) cos(9)j + sin((p) cos(0)k 


EN A 
É 


= cos(0) cot(0)i + cos() cos(0)j + sin() cos(0)k 
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cos2(8) cot?(0) + cos2(0)(cos2() + sin?()) 


/ 
= v/cos? (8) cot? (0) + cos2(0) 
N/ 


= 4/cos2(0) (cot2(0) + 1) 


= «/cos2(0) csc?(0) 


= cot(0) 


E para a coordenada y, 


EN A 
º 


ds sin(0) sin()) + sin(0) cos(y)k 


mo 
u(8,9) = 4/sin?(0) (sin?(69) + cos2(9)) 
= sin(6) 
Logo temos 
cy =cot(0)dO e oco =sin(0)dy 
Sabendo que do, = —w301, podemos encontrar! wa: 


dos = d(sin(0) dy) 
—ws01 = d(sin(6)) dy 
—ws cot(0) d9 = cos(0) do dy 


e então w3 = sin(0) dp. Portanto 
Kojos e dws 
K cot(0) sin(0) do dp = — d(sin(0) dy) 
K cos(0) dó dp = — cos(8) do dy 


Segue disso que a curvatura gaussiana da superfície DN é K = —1 
Exemplo 3.14. No exemplo (3.6) vimos que um toro tem parametrização dada 
por 
aq = dcos(p) + r cos(8) cos() 

y = dsin(y) + r cos(0) sin(y) 

2 =rsin(0) 

onde r é o raio do tubo e d é a distância entre o centro do tubo e a origem. A 


variação do vetor posição r em relação a O é dada por 
Õ N z a 
* — rsin(0) cos(q)i — rsin(0) sin(9)j + r cos(0)k (101) 


Tm 


HPoderíamos usar do; = w302, mas a única informação que obteríamos é que ou w3 = O 


ou w3 é proporcional a dy, já que do, = d(cot(0) dô) = 0. 
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Podemos verificar que 


or 
ee] ds 
00|” 
Também iremos encontrar que 
To) x A 
- = (—dsin(y) — r cos(0) sin(p))i+ (dcos(p) + rcos(0) cos(p))) (102) 
y 
To) 
u= = = d+rcos(8) 
e portanto 
oj=rdd e o=(d+rcos(0))dy 
Sabendo que do, = O, usaremos a equação dos = —w301: 


—rws dô = d((d + r cos(0)) dy) 
—rws dO = —r sin(0) do dy 


Segue disso que ws = — sin(0) dp. Sabendo que Koj0s = — dws, 


K .r(d + rcos(8)) do dp = — d(— sin(6) dy) 
cos(8) do dy 


Logo 
K .r(d + rcos(8)) = cos(8) 


K(0) = cos(8) 


= 103 
rd + r2 cos(0) quo, 


Lembrando que 6 é o ângulo entre a projeção do ponto em D e o plano xy (no 
sentido de fora para dentro do toro), temos que K > O quando —7m/2< 0 < 7/2, 
K <0 quandon/2<0<37/2eK =0 quando 6 = +7/2. Ou seja, a curvatura 
gaussiana do toro é negativa na parte interna, é positiva na parte externa e é 
zero nos pontos intermediários. 


No exemplo (2.1) vimos que o produto vetorial entre dois vetores ordinários 
uev do Rº pode ser dado por 


uxv==(UAV) 


De maneira semelhante, o produto vetorial entre vetores cujas componentes são 
p-formas é dado por 
du x dv = «(du A dv) 


onde o operador estrela * e o produto exterior A dessa equação age apenas nos 
versores é;, enquanto que a multiplicação entre as componentes é o produto 
exterior entre p-formas. 
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O produto vetorial de dr com ele próprio é 


dr x dr = «[(m& + 0962) A (o1& + 0262) | 


= *[01026,6» e 02016261 | 


= «[01026,6> sto 01026164 | 


= [201026164 | 

= 20109€3 
A 2-forma 0103 é o elemento de área de X, da mesma forma que 01/0903 era 
o elemento de volume no Rº. Já o elemento 01/0263 é o vetor de área de 5, 


que é o elemento da que aparece nas integrais de superfície. Por exemplo, em 
coordenadas esféricas com r = R constante, temos 


da = ogopf = R2 sin(0) do dyf 


Note também que o produto vetorial de um elemento com ele próprio só é 
zero quando a multiplicação das componentes comuta, o que não é o desse caso. 
De maneira semelhante, temos 


dr x dé3 = *[ (01& + 0265) A (mn + w263) | 
= (ow Ea OoWw1) 63 
- 2Hés 
e também 
dé3 x dé = *[ (w1& +uséa) A (ué + w262) | 
= Iw wa €3 
=2K6 


Agora seja f uma função bem comportada definida em uma superfície >. 
Podemos construir uma generalização do operador laplaciano que esteja restrito 
em 5, isto é, apenas os pontos de X são avaliados na variação de f. Esse 
operador recebe o nome de operador de Laplace-Beltrami, denotado por A f. 
Ele é definido da mesma maneira que o laplaciano, 


Af = «(d+df) 
mas com a diferença de que df agora é função apenas de o, e os: 
df = cjoy + c,09 
E então 


«df = Cj092 — C901 
dx df — d(c102 =— c201) 


A última equação pode ser expressa em termos da base o102, onde o fator 
multiplicativo é justamente A f: 


dx df = (Afjoios 
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4 Formas diferenciais e integração 


4.1 Variedades 


Seja o toro dado por 


aq = dcos(p) + r cos(8) cos() 
y = dsin(y) + r cos(0) sin(p) 


z = rsin(6) 


onde d e r são constantes positivas tais que d > r. A coordenada 0 é o ângulo 
entre o ponto (x,y, 2) e a sua projeção no plano xy, orientado de fora para dentro 
do toro pelo lado superior. Já y é o ângulo azimutal de (x,y,z) em relação ao 
eixo x. Iremos denotar M como o subconjunto do Rº formado pelos pontos 
que constituem o toro. Agora seja W o subconjunto do Rº dado por f(x,y) E 
R2, d-r<y224+y2<d+r). Em coordenadas polares, W pode ser descrito 
por (s,5), onde s = 12 +y2, x = scos(B) e y = ssin(8). Vendo o toro de 
cima, podemos associar cada ponto de W a um ponto da parte superior do toro 
(incluindo os equadores interno e externo) usando a projeção dessa superfície 
no plano xy. Denotaremos essa parte superior de U, que é o subconjunto de M 
com z > 0. Essa associação cria um mapa que chamaremos de dp e pode ser 
escrito como 


Pup: Wo U, 
Pup(s,B) = (a(0!(s,8),9!(s,8)), u(0!(s, 8), !(s,6)), 2(0(5,6),0!(s,8))) 


onde 


1 mu . (d-s 
0!(s,5) 5 + aresin( E ) 


comd-r<s<d+re-m<8<7. Otoroe sua projeção pode ser vista na 
figura 3. 

Da mesma forma, podemos criar um outro mapa Ggown que associa os pontos 
inferiores!2 do toro V com o mesmo conjunto U, mas agora com 


0(s,6) = 7 + arecos( E ') 
p(s,8)=8 


A união dos mapas Yup € Pdown cobrem toda a superfície do toro e portanto 
formam um atlas dele, isto é, os dois mapas são suficientes para descrever qual- 
quer ponto do toro. Os equadores interno e externo são comuns aos dois mapas, 
e como ambos usam o mesmo conjunto W, a correspondência de um para um 
dos pontos nesses equadores é trivial. 


PNote queUUV=M. 
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Figura 3: Toro com d=2er=1. A parte superior (domínio U) está colorida 
de vermelho e amarelo enquanto que a parte inferior (domínio V) está colorido 
de cinza. Em cima do toro, colorido de verde, está a projeção da parte inferior 
ou superior do toro no plano xy, que é o subconjunto W. 


Os mapas dup € Paown também admitem inversa. No caso Ria UsWe 
dinum : V +» W são dados por 


Pão (0,52) = GA y) (,y,2) e Vv+ 
E E] == (x,y) (2,y,2) eVT 


em coordenadas cartesianas. 

O fato de conseguirmos determinar relações unívocas entre quaisquer pontos 
de uma sub-região do toro (U ou V) com um subconjunto de um espaço euclidi- 
ano W € R? usando funções bijetoras contínuas Qup € Pdown Com inversas dio 
€ Pim também contínuas mostra que estabelecemos um homeomorfismo entre 
os dois conjuntos (W com U ou também W com V). 

Agora seja um ponto p qualquer do toro, isto é, pe M. Se U for um 
subconjunto aberto de M que inclua p, então dizemos que U é uma vizinhança de 
p. Se a vizinhança U de p for conexa e pequena o suficiente, podemos estabelecer 
um homeomorfismo entre U e um subconjunto de um espaço euclidiano, que no 
caso é o R2. Para os pontos fora dos equadores, isso pode ser feito usando 
os mapas descritos anteriormente. Agora note que nenhum dos dois mapas, 
Pup € Pdown Vão sozinhos cobrir completamente uma vizinhança conexa Uç de 
um ponto que fica em um dos equadores. Porém nesse caso podemos sempre 
criar um novo mapa que é a projeção lateral da superfície no qual o equador 
se encontra, de forma que possa cobrir toda a vizinhança Uç e estabelecer um 
outro homeomorfismo com um subconjunto do R?, que seria o plano yz ou xz. 
Isso significa que todos os pontos do toro possuem alguma vizinhança tal que é 
homeomórfica com algum subconjunto aberto! de um espaço euclidiano. 


8Com a adição dos mapas laterais para cada equador, podemos restringir W de forma que 
não cubra os equadores (d—- r <s<d+r) e assim W seja um conjunto aberto do R2. 


108 


Essa característica não é única dos toros, pois outros espaços topológicos 
tais como superfícies esféricas, retângulos, discos e curvas parametrizadas que 
não se cruzam do tipo (xi(t),...,zn(t)) com x;(t) contínuas também possuem 
a propriedade de que todos os pontos que os compõe possuem uma vizinhança 
que é homeomórfica com algum subconjunto aberto do R”. Isso nos motiva à 
seguinte definição: 


Definição 4.1 (Variedades). Um espaço topológico M é denominado de vari- 
edade se todos os pontos que o compõe possuem uma vizinhança que é homeo- 
morfa a um subconjunto aberto de um espaço euclidiano. 

Se em um determinado ponto pe M a menor dimensão do espaço euclidiano 
tal que é homeomórfica à sua vizinhança é n, então dizemos que a variedade M 
tem dimensão local n no ponto p. Se todos os pontos de M tiverem a mesma 
dimensão local n, então denominamos M de variedade n-dimensional, ou de 
variedade de dimensão n. 


Dada uma variedade de dimensão n, uma vez que cada ponto possui alguma 
vizinhança homeomórfica ao R”, segue que são necessários n parâmetros para 
localizar esse ponto nessa vizinhança. 

O toro M descrito anteriormente é uma variedade de dimensão dois porque 
qualquer homeomorfismo de um espaço euclidiano com uma vizinhança de qual- 
quer um dos pontos necessita que esse espaço seja pelo menos o R2, não sendo 
possível cobrir uma vizinhança usando apenas a reta real R ou apenas por um 
ponto (dimensão zero). 


4.1.1 Diferenciabilidade 


Dada M uma variedade de dimensão n, a junção dos mapas que fazem home- 
omorfismo com R” e que conjuntamente cobrem todos os pontos formam um 
atlas de M. Sejam U e V subconjuntos abertos de M que possuem pontos em 
comum (isto é, UnV 9). Segue que dois mapas homeomórficos q: U +» R” 
ey: Vi» R” pertencentes a esse atlas possuem um domínio em comum, que 
éUnvV. Uma ilustração disso é mostrado na figura 4. Os mapas também são 
denotados por (U, &) e (V,1)). 


Figura 4: Dois subconjuntos U e V da variedade M possuem uma região em 
comum U nV representado pela cor violeta. 


Um exemplo é o caso do toro, em que U pode ser a parte superior do toro 
(não incluindo os equadores) e V uma lateral da superfície externa (excluindo 
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os pontos em que z = +r). Evidentemente há uma região em comum entre U 
e V que são os pontos dessa lateral com z e (0,7). E então & leva os pontos 
superiores do toro à projeção xy enquanto que Y) leva os pontos da lateral 
externa à projeção xzz ou yz, se essa lateral for paralela a um dos eixos x ou y 
respectivamente. A região U n V é ilustrada na figura 5. 


Figura 5: Região do toro no qual Un V. O subconjunto do Rê no qual é 
homeomórfico no domínio Un V aparece ao lado esquerdo. 


Qualquer ponto na região em comum U n V pode ser aplicado tanto por q 
quanto por 1). De uma outra perspectiva, esses pontos podem ser uma ponte 
intermediária entre a imagem de q e de y. No caso q! leva a imagem de q 
dessa região até um ponto em Un V e então aplicamos 1) para levar esse ponto 
até a imagem de y. Essa passagem é representada pela composta yo q! : 
AUnNnV) » y(U nV). Como d e 1”! são homeomorfismos, y o q! também 
é. A função yo q”! é denominada de mapa de transição. 

Agora seja f: M > R uma função real sobre a variedade M. Podemos levar 
a imagem (UnV) até U nV usando q! e então aplicar f, no que resulta 
na composta f o q!. Analogamente podemos construir fowW-l. Se M for 
uma variedade de dimensão n, então as imagens (Un V) ey(U nV) podem 
ser descritas usando n parâmetros, chamados de coordenadas locais. Iremos 
denotar as coordenadas em G(U NV) de xº e as coordenadas de y(U AV) de yº. 
Com base no mapa de transição yo q”, há uma relação entre essas coordenadas 
locais, no qual pode ser expressa por 

g=r(yim,.ouy) e y=y(al,g,...,m”) 
Disso pode surgir a necessidade de saber como que f se comporta em relação 
às coordenadas locais, e isso envolve derivadas de fem «º' cem y'. Se fog 
e f ow! forem diferenciáveis pelo menos até primeira ordem, isso é possível. 
Mas mesmo assim, isso não garante que possa existir uma regra da cadeia entre 
x' e y'. Usando as propriedades de composição podemos escrever fo q”! na 
forma 


fog! =(foyDo(pog) 
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Então usar a regra da cadeia na derivada de fo q”! em xº requer que o mapa 


de transição y o q”! seja diferenciável. Por exemplo, se isso for verdade até 
Oyl Ox. 

- e — existem e a regra da 

dx" Oy) 


à primeira ordem, então as derivadas parciais 


cadeia (de primeira ordem) também. 

Com base nisso, introduzimos o conceito de atlas diferenciável de M, no 
qual é um atlas de M composto por mapas tais que as transições entre eles 
sejam todos diferenciáveis (se tal atlas existir). Isso significa que se f for dife- 
renciável nas coordenadas «x* em U n V, então pela regra da cadeia f também 
será diferenciável em y na região U n V. 

Agora suponha que uma variedade M de dimensão n suporte um atlas di- 
ferenciável. Podemos adicionar um novo mapa (Uo, do) a ele. Se o novo atlas 
com esse mapa continua sendo um atlas diferenciável (isto é, as transições dos 
mapas já existentes com (Uso, 49) são diferenciáveis), então dizemos que o mapa 
(Uo, bo) é diferencialmente compatível com tal atlas. O conjunto de todos os 
mapas de M diferencialmente compatíveis formam um atlas máximo de M. 

De agora em diante iremos trabalhar apenas em variedades que admitem um 
atlas máximo, o que nos motiva à seguinte definição: 


Definição 4.2 (Variedade diferenciável). Uma variedade diferenciável é uma 
variedade dotada de um atlas máximo. 


Uma função f: M » R é dita ser diferenciável em um ponto p de uma 
variedade diferenciável M se as derivadas de fo qd! existirem para qualquer 
mapa à do atlas máximo que cubra o ponto p. Uma função é dita suave se ela 
é diferenciável até a ordem que precisarmos. 


4.1.2  Orientabilidade 


Seja M uma variedade diferenciável e (U,4) um mapa do atlas máximo com 
coordenadas x'. Dado que as coordenadas xº são independentes entre si, segue 
que 
QRO. 5 
Ba th 


Usando a regra da cadeia, temos 


Fo oyi 
Fazendo a;; = -eCbix= obtemos uma equação matricial 
7 dy! Ot” 


Di uisby = O 
j 


ou então AB = I, onde 4 = [a;;] e B = [b;x]. Isso significa que det(AB) = 
det(A) det(B) = det(1) = 1, ou seja, det(4) £ 0 e det(B) £ 0. Mas Ae B são 
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justamente as matrizes jacobianas 


dar dar og" og" 
out yr ôxl "dam 
A = E st cê e B= e mo 
dx” da” Oy” Oy” 
Segue disso que 7 a entre as coordenadas locais dos mapas garante 
Orsa 


um jacobiano = det(A) não-nulo entre as coordenadas de dois 


Ay, ... gy”) 
mapas. Como os mapas de transição são diferenciáveis pelo menos até primeira 
ordem, então o jacobiano é uma função contínua que não tem raízes. Pelo 
teorema do valor intermediário, não é possível que o jacobiano troque de sinal, 
já que se isso acontecesse, teria que se anular em algum ponto. Portanto o 
jacobiano ou é sempre positivo ou sempre negativo. 

Sabemos que o determinante muda de sinal ao trocar uma coluna pela outra, 

então a ordem das coordenadas é importante. Quando o jacobiano é positivo, 
dizemos que o mapa de transição preserva a orientação. Uma variedade dimensi- 
onal é dita orientável quando admite um atlas composto de mapas de transição 
que preservam a orientação, também chamado de atlas orientado. 
Exemplo 4.1 (Circunferência). Seja M o subconjunto do Rº constituído pe- 
los pontos (x,y) que satisfazem q + y2 = a2. Esse conjunto forma uma 
circunferência de raio a centrado na origem. Agora iremos criar um mapa 
dup: Ur (-a,a), onde U é o subconjunto de M cujos elementos satisfazem 
y>0. Definiremos dup por 


Puplt, y)=a (y>0) 


e que cuja inversa 471 é dada por 
up 


bug (2) =(2,va2-— a?) 
Como Gup é bijetora, é contínua e tem inversa contínua, então dup forma um 
homeomorfismo entre U e um subconjunto aberto dos reais R. 
De maneira totalmente análoga, definiremos outros três mapas, rights Pdown 
€ Qleft 
Prignt(x, y) =Yy (a as 0) 
Pdown(L,Y) = « (y<0) 
dregi(t,y) =Yy (a < 0) 


com inversas dadas por 


Esses quatro mapas cobrem todo o conjunto M e portanto formam um atlas de 
M. E como todos são homeomórficos com algum subconjunto de R, segue que 
M é uma variedade de dimensão 1. 

Os mapas de transição são quatro e cada um corresponde a um quadrante 
do plano R2. Por exemplo dup O Priaht é dado por 


Pup o Prigha(U) = dup(N/a2 — 2, y) = Va? — y? 


Agora percebemos que o jacobiano 


d 4 da y 
gr (o) righe(Y) dy V2-yà 


é negativo para todo y na região de interseção y e (0,4). Como a circunferência 
possui dimensão 1, não temos o luxo de trocar coordenadas de posição, já que só 
existe uma. Isso significa que esse atlas não pode ser orientado, mas não implica 
que a circunferência seja não-orientável. Ao invés disso, iremos escolher novos 
mapas p— e b+ dados por 


y a es 
a+T e da (x,y) = a-t 


do-(x,) = 


Denotando s a coordenada local da imagem de &— et a coordenada relacionada 


a P+, temos 


a+ ax a—r 


O significado de s é a inclinação da reta (em relação ao eixo x) que passa 
pelo ponto (-a,0) e em (x,y) E M, como pode ser visto na figura 6. Ját 


(x,y) 


Figura 6: O parâmetro s é a tangente do ângulo a. 


é a inclinação da reta (em relação ao eixo x) que passa em (a,0) e em (x,y). 
Podemos observar que P— cobre todos os pontos da circunferência exceto (—a,0) 
de maneira unívoca, enquanto que 4, faz o mesmo exceto em (a,0). Então 
ambos formam um atlas da circunferência. Em seus respectivos domínios, vemos 
que ambos os mapas são contínuos. Ambos s et podem assumir quaisquer valores 
reais. O domínio em comum entre &— e, são todos os pontos da circunferência 
exceto (—a,0) e (a,0). 


113 


2 


Usando o fato de que «2 + y? = a? e usando a primeira expressão de (104), 


chegamos na equação 

q —-r = s(a+w) 
Uma das soluções da equação do segundo grau para x(s) é a(s) = —a na qual 
não nos serve por não ser uma inversa adequada. A outra solução é 


1—s? 


TE 


2 


Substituindo esse resultado em y? = a? — x?, para y > 0 obtemos 


Portanto a inversa de P- é 


1-2 2s 
— 
s)=[a- a- 
pane) ( 1+s? E) 
que é uma função contínua. Logo concluímos que &— é um homeomorfismo com 
R. De maneia totalmente análoga, podemos mostrar que o mapa 4, também é. 


O mapa de transição t(s) ou s(t) pode ser obtido diretamente de (104) ao 
multiplicar uma equação pela outra: 


4 y? O me 
(a+a)(a— x) a? — q? 
e portanto 
1 
t(s) = —— 
()=—5 
O jacobiano é 
dt 1 
ds 8? 


que é sempre positivo na região do domínio em comum (t+ 0 es 0). Com 
isso, concluímos que o atlas formado por &- e &, é orientado e portanto a 
circunferência é orientável. 


4.1.3 Vetores tangentes 


Seja f: M > R uma função suave na variedade diferenciável orientável M de 
dimensão n. Seja também U um subconjunto aberto de M no qual um ponto 
pe M lhe pertence e xº sendo suas coordenadas locais. Uma parametrização 
dessas coordenadas pode ser dada por y' em um domínio V tal que x'(y')E U e 
no qual a regra da cadeia entre essas coordenadas existe e é contínua. Podemos 
obter f pela parametrização fazendo f(x'(y?)). Com isso, temos que 


of mofo 2 Ox O 
dyi Mm dat dy [> Oyj E) / 


i=1 
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i 
Zu sede ; Eu ; 
No somatório nos deparamos com um elemento matricial a; = EN de maneira 
y 


que 


ô de cuca 
f e f 


Oyj Rn Que 


o j Sa 
Isso significa que o operador diferencial Ei pode ser construído como uma 
a 


das e 9 
combinação linear de outros operadores —: 


Ox 
Fo) 140 
d = Dai (105) 
i=1 


Com base nisso, podemos afirmar que os operadores diferenciais em xº formam 
uma base de um espaço vetorial T de dimensão n!4. 
Agora suponha que a parametrização y' tenha n variáveis independentes e 


o seu domínio V seja um subconjunto aberto de M no qual contém o ponto p 


' ENRESA ” : Ra (o, õ 
referido no início. Então os diferenciais —— também formam uma base de T e 


Oy 
a; se torna o jacobiano entre x! e y). Cera: suponha que as derivadas em q) 
sejam dadas pela equação (105). Então a jacobiana A = a; é a matriz mudança 
de base das derivadas em x” para as derivadas em yº. 

Dado um operador diferencial qualquer v e T, ele pode ser escrito nas duas 


bases: 
RR E (106) 
Fobia 1 Oy) 


i=1 


onde u” e w? são as coordenadas de v em suas respectivas bases. Usando (105), 
obtemos 


2 j 0 
v= Ww dj 
ox! 


j=1  i=1 
n n 

Ig? Ô 

= |vra] a 

i=1 Aj=1 Ee 


Comparando esse resultado com (106) e sabendo que os elementos da base são 
linearmente independentes entre si, chegamos na equação 


n 
O; 
tl > uiai 
j=1 


HA independência entre as coordenadas garante que esses elementos são linearmente in- 
dependentes entre si. No caso do R?, se em algum ponto de uma superfície as suas duas 
coordenadas forem múltiplas uma da outra, isso resulta em uma região pontiaguda e o jaco- 
biano se anula ou fica indefinido nesse ponto. 
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Se X for a matriz coluna contendo as coordenadas u' e Y a matriz com as 
coordenadas w?, concluímos que 


X = AY 


ou então!9 
AIX=Y 


- O nu j . 
Pela definição (1.10), temos que os vetores —. são contravariantes, o que jus- 
oz 


tifica usarmos o índice superior para as coordenadas u” ou wi de v. 


Considere o ponto p no qual as coordenadas locais xº são dadas por (x!,...,a”) lp 


(c!t,...,c”). Seja F(M) o conjunto de todas as funções suaves do tipo M + R. 
Na região U estão definidos os operadores diferenciais + á onde agora essas 
derivadas são calculadas no ponto p. Segue disso que aplicar esses operadores 
em alguma função f e F“(M) 


" Og? 


, retornará um número real. Essas deri- 
p 


vadas parciais em p nos lembra dos vetores tangentes da seção 3.2, o que nos 
motiva à seguinte definição: 


Definição 4.3 (Vetores tangentes). Um operador v : F(M) » R é deno- 
minado de vetor tangente em p se v satisfazer as seguintes propriedades para 
quaisquer f eg de F(M): 


f. v(af+bg) =av(f) +bv(g), VabeR; 
2 v(f:9) = 9(p)vlf) + H)víg). 


Õ ; 
Podemos observar que os operadores | de fato satisfazem ambas as pro- 
Llp 


priedades, no que podemos chamá-los de vetores tangentes. 
A totalidade dos vetores tangentes em p forma um espaço vetorial que de- 


A 


[9 
notaremos de T;M. Além disso, o conjunto (a 


Ox 


| forma uma base desse 


p 
espaço. 


Se f = c for uma função de Fº(M) onde c é uma constante, então para um 
vetor tangente v e uma outra função qualquer g e Fº(M) temos 


v(f-9) = 9(p)v(f) + flp)víg) 


Mas f é constante e então v(f- 9) = v(c- 9) = cv(9). Logo 


cv(g) = g(p)v(f) + Flp)v(g) 
= g(p)v(f) + cv(g) 


1SVimos na seção anterior que o determinante da jacobiana é diferente de zero e portanto 
AI existe. 
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Como isso é válido para qualquer função g de Fº(M), segue que 


v(/)=0 


Isso significa que todos os vetores tangentes em p possuem a propriedade das 
derivadas de que aplicá-los em uma função constante resulta em zero. 

A matriz jacobiana responsável por mudar as coordenadas de uma base de 
vetores tangentes para outra no ponto p agora é dada por 4 = [ai], onde 
Ox 
a; = ais 

Yºip 

Denotando r = (x!,...,x”) como sendo um ponto qualquer de U em termos 

das coordenadas locais x, a expansão em série de Taylor até primeira ordem de 


uma função f e Fº(M) é dada por 


of % KA 
19=10+5 a) (6) 
onde c = (cl,...,c") representa o ponto p. Aplicando um vetor tangente v em 
f obteremos 
of % (A 
v = viglc))+ =| vg —c 
(1) = víf(e)) 2 daily ( ) 


Denotando v(x' — c”) = v(x!) = v* e sabendo que v(f(c)) = 0, temos 


VE A” Og 


II 
G 
E 
Se 
S| 
io 
“—— 
Era 


Como isso é válido para qualquer f, temos a representação do vetor tangente v 


em p em termos da base 4 =——| +»: 
Og p 


n a] 
; O 
v=) 2 
Ox? 
i=1 


Sabendo que esse resultado é análogo para qualquer outro ponto p da varie- 
dade M, podemos definir um campo vetorial em M onde em cada ponto p da 
variedade está definido um vetor tangente nele e esses vetores são suaves na 
variedade, isto é, variações infinitesimais em um ponto de M são acompanhados 
por variações também infinitesimais dos vetores tangentes. Nas coordenadas 
locais de U podemos escrever o campo vetorial v(r) na forma 


p 


4.2 Formas diferenciais 


Conforme vimos na seção anterior, denotamos Fº(M) o conjunto de todas as 
funções suaves na variedade diferenciável M. Agora também dizemos que os 
elementos desse conjunto são denominados de 0-formas e que Fº(M) é o espaço 
das formas de grau zero em M. 

Seja p um ponto em M. Sejam U e V subconjuntos abertos não-vazios de 
M em que ambos são uma vizinhança de p com coordenadas locais xº e yj 
respectivamente. Definimos as 1-formas em p como sendo objetos w que podem 
ser escritos nas bases (dx') e (dy), 


n n 
w= x u; der' = s w; dy? (107) 
i=1 j=1 
onde as coordenadas u; e w; se transformam pela equação 


ui = > wjb (108) 
j=1 


Oy ê 
comb! = | = | gy). 
O dailo dailp” 
Podemos realizar a transformação inversa de (108) ao multiplicar ambos os 
lados por aí. e aplicar um somatório em i: 


n n i 
spp Oy = 
RR Ox? Oy" Oy" Ê 
i=1 i=1 yr y 
e portanto 
n 
o UM = > MW; = Wk 
i=1 j=1 
na ; ú ar 
Isso significa que a matriz A = [ai] usada!º para levar as bases an Para o é 
JL y 


a mesma usada para levar as coordenadas na base dx* para dy. Ou seja, pela 
definição (1.11) as 1-formas em p são vetores covariantes, ou covetores. 


Ágora seja e; = um vetor tangente em p. Então 


Ox p 
Ben a 
ey”) = Ox " = bj 
Õ o Õ 
16Lembrando da equação (105 — = i— 
embrando da equação (105) que 27 A 0 3a 
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Da equação (108) temos 


Substituindo em (107) na base dz”, ficamos com 


ns, [> ese) da' = 5, ww; [> eae) 


i=1 ÀAj=1 =, 


Comparando com (107) na base dy” e sabendo que as 1-formas em p são line- 
armente independentes, então 


dy) = 3  es(y!) da” (109) 
1=1 


Ágora iremos definir uma relação bilinear dada por 


(der, e) = 


(110) 


A) d 
Xp 


de maneira que a base fey,e»,...,en+ seja o dual da base (da! dx2,...,da”) e 
vice-versa. 

Além disso, aplicando essa relação na equação (109) com o vetor tangente 
em p ey teremos 


(dy, ex) = >, es(y!) (dat,ex) = >, es(y!)ôj — ex(y) 


e portanto 


n 
dy) = >. (dy), ei) da 
i=1 
Multiplicando ambos os lados por w; e aplicando uma somatória em j, iremos 
obter um resultado que vale para qualquer 1-forma w em p: 


w= 3 w; dy! = s 3 Wj ao ja = s (w, es) da 
j=1 i=1 Nj=1 i=1 


O produto exterior entre q 1-formas em p geram g-formas em p. Com isso, 
podemos definir q-formas em uma variedade diferenciável M ao associar uma 
g-forma em p para cada ponto p de M de maneira que as suas coordenadas 
variem suavemente em M. Por exemplo, para uma g-forma w em M dada por 


w=S unltd H=(h,...,ha) 
H 


necessita-se que ug (r) sejam funções suaves. O conjunto desses objetos formam 
FM), que é o espaço vetorial das g-formas em M. O produto exterior de uma 
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g-forma com uma p-forma gera uma (p + q)-forma, que é obtida ao realizar o 
produto exterior entre as formas para cada ponto de M. 
Sabendo que 


n 


dat = ) tdo E oa dy! = [Da Ee 


= 


podemos definir uma derivada exterior d dado (na base das coordenadas em xº) 


por 
d= da” Ê 
E E Og 


no qual leva uma p-forma para uma (p + 1)-forma do seguinte modo: 


Ami 
i=1 Cl, 


d(f dae”) = [> ud) Pane e a ndnÊs df a de” 
1=1 se 


Para uma p-forma mais geral w dada por 


temos que 


dw = Ps fam dx” À daf 
H 


que é a mesma expressão que aparece na demonstração do teorema 3.2 e portanto 
a derivada exterior d satisfaz os quatro axiomas (6, 7, 8 e 9) vistos na seção 3.3. 

Se M e N são variedades diferenciáveis de dimensão m e n respectivamente 
ed: M +» Né um mapa suave entre as coordenadas em M e em N, então 
podemos definir uma função q5 : FP(N) +» FP(M) que converte uma p-forma 
em N com coordenadas locais em N para uma p-forma em M com coordenadas 
locais de M. Conforme vimos na seção 3.4, essas funções satisfazem 


1 di(w +) = GE(w) + G5(7) 

2. Pprglw A À) = Gi(w) A dG(A) 

3. d(ójw) = Gp (do) 

4. Seg:iU»Vey:VoW, então (pod) = droyk 
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4.3 Simplexos 


Definição 4.4 (Envoltório convexo). Seja X um conjunto não vazio de pontos 
no R?. A interseção de todos os subconjuntos convexos do R” que contém todos 
os pontos de X formam um envoltório convexo de X. 


Um exemplo de envoltório convexo no espaço R? é a analogia de uma coleira 
elástica que se molda nos pontos mais externos”de X. Na figura 7 temos 
cinco pontos em azul no R? que formam um conjunto X. E então usamos 
três círculos (que são subconjuntos convexos) de diferentes raios e centros que 
englobem X e pintamos a interseção entre esses subconjuntos de amarelo. Já 
na figura 8 colocamos novos quatro círculos com raios maiores que englobem X 
e então a interseção entre os sete círculos nos dá uma noção de como seria o 
envoltório convexo de X. Podemos também usar outros subconjuntos convexos 
como retângulos, triângulos entre outros para que a região em amarelo fique 
ainda mais parecido com o envoltório convexo. 


Figura 7: Três subconjuntos conve- 


Figura 8: Sete subconjuntos conve- 
xos contendo os pontos de X. 


xos contendo os pontos de X. 


O envoltório convexo desses pontos é um paralelogramo cujos vértices são os 
quatro pontos mais externos. O próprio paralelogramo é um subconjunto con- 
vexo e portanto também pode ser usado para englobar X. Esse paralelogramo 
acaba tornando-se o menor subconjunto convexo do R? que contém X. 

Seja X = (py,..., pk; um conjunto de pontos no R”. Se u;; são os vetores 
que ligam o ponto p; até ps, 


ui; = r(pi) — r(p;) 


forem linearmente independentes para quaisquer à e j, 


k 
V9: Djau;=0 => a; = 0 Va 

i=1 

1%) 
então dizemos que X é um conjunto de pontos independentes. Por exemplo 
em um conjunto X com três pontos independentes, esses pontos não podem ser 
colineares. Com quatro pontos independentes, estes não podem pertencer a um 
único plano simultaneamente. 
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Definição 4.5 (Simplexos). Seja X = (po,...,pky um conjunto de k+1 pontos 
independentes do R”. O envoltório convexo de X é denominado de k-simplezxo. 
Esse k-simplexo é denotado pela expressão [po,..., Pk). 


Usando combinação convexa pode-se mostrar que um n-simplexo s” que é 
o envoltório convexo de X = (po,...,Ppn) pode ser expresso pelo conjunto de 
pontos dado por 


sr= [tomate cttum : Des t;>0 “a (111) 
i=0 


Isso significa que um 1-simplexo é um ponto; um 2-simplexo é um segmento de 
reta; 3-simplexo é um triângulo e um 4-simplexo é um tetraedro. 

Dizemos que um n-simplexo é orientado quando permutações ímpares entre 
os pontos inverte o seu sinal, enquanto que permutações pares mantém o sinal. 
Por exemplo 

[po,P1,22] = —[p1,po,72] = [p1, 2, Po] 


A fronteira de um n-simplexo orientado c é denotado por dc e é definido 
pela expressão 


olpo, .. »Pn] = > (-1/lpo, cs Pj-1,Pj41,--- +Pn] (112) 
7=0 


o que significa que a fronteira de um n-simplexo é a soma!” de vários (n — 1)- 
simplexos orientados. 


Exemplo 4.2. Seja o 2-simplexo orientado dado por c = [po,p1,p2| represen- 
tado na figura 9. A ordem po, p1 e pa indica que o triângulo está orientado no 
sentido anti-horário, mostrado pela seta circular. A fronteira desse simplexo é 
dado por 


oc = [p1,p2] — [po, p2] + [po, pi] 
= [p1,p2] + [p2,po] + [po,p1] 


Isso significa que as fronteiras desse triângulo são segmentos de reta orientados 
que ligam os vértices, denotados pelas arestas com as setas. 


Teorema 4.1. Se c é um simplexo orientado, então d0c = 0. 


Demonstração. Suponha que c é um n-simplexo dado por 


C= [gg 25:50] 


ITA soma entre simplexos deve ser entendida como a união dos subconjuntos de pontos do 
espaço euclidiano que compõem cada simplexo. 
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PO 
P1 


P2 


Figura 9: Um 2-simplexo orientado e suas fronteiras. 


Aplicando novamente O na equação (112) e sabendo que a fronteira da soma de 
simplexos é a soma das fronteiras dos simplexos, temos que 


00c = = (—1)0[vo, cs UG, UG4I, ee Un] 


j=0 
n n—1 
j i 
= > (— 1)? x (—1) [vo, dá av usg Vj-1, UVj41; ...14 Vn]i 
j=0 i=0 
onde [v0,...,Vj-1,UVj+1,---, Un]i denota que o i-ésimo ponto de [vo,...,Vj-1,UVj+1,--. 


foi suprimido!8. Sabemos que o resultado será uma soma de (n — 2)-simplexos, 
ou seja, dois pontos distintos foram suprimidos do n-simplexo original em cada 
termo da soma. Agora suponha dois pontos vp e vg. Sem perda de generalidade, 
iremos presumir p < q. Há duas formas de obter um (n — 2)-simplexo sem esses 
dois pontos do n-simplexo original: eliminando primeiro v, e depois v, ou o 
contrário, primeiro vp e então vg. Eliminando v, primeiro, segue que o p-ésimo 
ponto do (n — 1)-simplexo resultante será vp. Então esse (n — 2)-simplexo é 
gerado pelos índices ) = q ei = p e deve aparecer um termo na soma do tipo 


(—D)P[vo,... 4Up— 15 Up4 1, ces Ugo 15 Vg+ 1, Un] 


Por outro lado, se eliminarmos primeiro v,, então o q-ésimo ponto do (n — 1)- 
simplexo resultante será v;+1. Logo o (q — 1)-ésimo ponto será v, e consequen- 
temente um (n — 2)-simplexo dado por 


(—1)P+(a-D Era RR DE E ORE PRA] 


é gerado com os índices j = pe i = q-—1. Porém um tem sinal (—1)*P e o 
outro (—1)?+(9-1) = —(—1)P+9, Como os sinais são diferentes e ambos possuem 
a mesma ordenação, então eles se anulam na soma. Como isso é válido para 


qualquer par de pontos distintos vp e vg, Segue que a soma total é zero, isto é, 
00c = 0. E 


Em um n-simplexo s” dado pela expressão (111), os coeficientes t;, cuja 
soma é 1, são denominados de coordenadas baricêntricas de s”. 


*8Note que o i-ésimo ponto desse (n — 1)-simplexo não necessariamente é v;. 
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+ Un] 


Queremos estabelecer uma relação entre os pontos de dois simplexos distin- 
tos usando suas coordenadas baricêntricas. Para isso precisaremos do próximo 
teorema. 


Teorema 4.2 (Unicidade das coordenadas baricêntricas). Seja S o conjunto 
dado por 


n 
S = [tom + ot te ; Bju] 


i=0 
onde X = fpo,P1,.-.,Pn+ é um conjunto de pontos independentes em R*. Então 
existe uma correspondência biunívoca entre os pontos de S e as (n+1)-uplas 
(to,ti,..., tn) 
com 


s t=1 (113) 
i=0 


Demonstração. Suponha que as (n + 1)-uplas (to,ti,...,tn) e (tos t4,...sth) 
restritas às condições (113) representem um mesmo ponto em S, ou seja, 


toPo =: tips +...+ tnPn = toDo + pr +...+ Ea e R* (114) 
Segue disso que 
Dt — ipi =0 
i=0 
De maneira conveniente, reescreveremos essa expressão como 
n 
(to — to)po + BR —t)p=0 (115) 
i=1 


Como os pontos de X são independentes, temos que 


Is 


ai(pi p;)=0 a; = 0Vi 
0 
j 


KI 


Como isso é válido para todo j, sem perda de generalidade faremos j = 0. 
Queremos saber o valor de 


Bt — po) 


pois se for sempre nulo, teremos que t; — t; = 0, que é o que queremos. Expan- 
dindo o produto na somatória, temos 


n n 

> (t; — t)(p; — Do) = o t;D; + t;po — tipo — tp; 

i=1 1=1 
n n is: 
= 5 (ti — pi + po [Bu Se) 
1=1 1=1 1=1 
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Mas de (113) e (115) temos 


RU ti) (p: — po) = — ((to — to)po) + po ((1 — to) — (1 — to) 


= —L(to — to)po + (to — to)po 
=0 


e portanto t; — ti = 0 ti=t Vi=l1,2,..,n. Isso significa que os n 
últimos termos de cada membro da equação (114) são idênticos e então 


topo = topo to=t 


Logo para cada ponto em S há apenas uma única (n+1)-upla (to, ...,tn) restrita 
às condições (113). E 


No caso particular dos simplexos em que t; > O Vi, o teorema (4.2) se 
aplica. Agora sejam dois n-simplexos orientados dados por s” = [po,...,Pn] € 
t” = [go,...,Gn]. Podemos criar um mapa entre os pontos de s” e de t” fazendo 


n n 
tip; = px tiqi 
i=0 i=0 


Como as coordenadas baricêntricas são únicas para cada simplexo, então essa 
relação é biunívoca. Além disso, se t; = 0;;, temos que 


PD; q VI 


Isso significa que esse mapa preserva a ordem dos vértices dos dois simplexos. 

Sabendo disso, podemos definir uma integral de uma p-forma em um sim- 
plexo específico A” e então calcular essa integral em outros simplexos quaisquer 
com base em A”. 


Definição 4.6 (Simplexo padrão). Um n-simplexo cujos vértices inclui a ori- 
gem enquanto que os demais são representados pelos vetores da base canônica do 
R” é denominado de n-simplexo padrão, que é denotado por A”. O n-simplexo 
padrão pode ser escrito na forma 


aro foto | NES > 0ve 
1=1 


Note que o n-simplexo padrão é definido!” apenas no espaço R”, ao contrário 
de outros n-simplexos que podem existir em R*, com k > n. 

O n-simplexo padrão é o simplexo cuja origem é um de seus vértices e é a 
ponta de um cubo n-dimensional. Por exemplo A? = [(0,0); (1,0); (0,1)] 


Ox 


Em outros textos um n-simplexo padrão é definido substituindo diciti <lpordiot = 
1e R” por R?ºt1 na definição (4.6), o que exclui a origem e diminui a dimensão em 1. 
Entretanto não usaremos essa definição aqui. 
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um triângulo retângulo com catetos unitários no R? cujo ângulo reto está na 
origem e portanto é a ponta de um quadrado. Podemos equivalentemente falar 
que A? é o subconjunto fechado do R2 cujas fronteiras são os eixos 1 e y e 
areta ar +y = 1. Analogamente Aé = [(0,0,0); (1,0,0); (0,1,0); (0,0,1)] 
é o subconjunto fechado do Rê delimitado pelos planos xy, yz, zz e o plano 
cty+z=1. 

4.4 Integração de formas diferenciais 

Iremos definir a integração de uma n-forma w = fd! dx2...dr” em um n- 
simplexo padrão A” (no qual w está bem definida) pela equação 


/ tu = E RADAR a o ra > O a 
ú Ra 


Por exemplo, se w = dx dy então 


|, l-a ah 
/ w= dody = [ K dyda = = 
A? A? o Jo 2 


no qual também pode ser escrito na forma 


fas = ML “anjos 


Em geral, para uma (p + 1)-forma do tipo w = Adal...da? daP*! e o simplexo 
padrão AP+1 teremos 


1-SEn o 
/ u=[ E AdaP* | dal... da” 
Ap+ ar Jo 


onde A? é o p-simplexo padrão envolvendo as coordenadas x!,...,a?. 

Agora sejam s” = [po,...,;Pn) et” = [g0,...,gn] dois n-simplexos orienta- 
dos. Queremos uma vizinhança U que inclua o simplexo s” e uma vizinhança 
V de t”. Entretanto os simplexos são conjuntos fechados e necessitamos que U 
e V sejam conjuntos abertos para não haver problemas relacionados à derivada 
nas fronteiras. 

Para resolver esse problema, iremos definir dois espaços planos W e Z dados 


por 
n 
= [tomo citam ; Sua] 


2 = [tomt ot ta Sua] 
i=0 


O que diferencia W e Z dos simplexos é que as coordenadas t; podem ser 
quaisquer reais, não necessariamente não-negativas. Por exemplo, se s” for um 
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segmento de reta, então W é a reta que o inclui. Se s” for um triângulo, então 
W é o plano que o contém e assim por diante. 

Dessa forma, definimos U e V vizinhanças abertas contidas no menor sub- 
conjunto de W e de Z que inclua s” e t” respectivamente. 

Agora seja M uma variedade diferenciável de dimensão n. Podemos construir 
um mapa suave à que leva um ponto de U até M, 


d4:USM 
e um segundo mapa também suave Y dado por 
viVoM 
Com isso definimos dois objetos 


(s”, U, É) e (t”, V, 1) 


que relacionam os dois simplexos com a variedade M. Dizemos que ambos os 
objetos são equivalentes se para as mesmas coordenadas baricêntricas os mapas 
fe levam UeV respectivamente ao mesmo ponto em M, ou seja, 


n n 
(o) a tip; |=U se tiqi 
i=0 i=0 


Todos os objetos do tipo (s”, U, &) que são equivalentes entre si formam uma 
classe de equivalência que denotaremos por o” e chamaremos de n-simplexo em 
M. 

Uma vez que sempre podemos estabelecer uma relação de um para um entre 
dois simplexos quaisquer, então para um n-simplexo o” qualquer em M sempre 
podemos construir um mapa À que mapeará os pontos da vizinhança U de um 
simplexo padrão A” para a variedade M e que simultaneamente (A”, U, q) 
pertença à classe 0”. 

Obtido tal mapa É: U > M, podemos construir uma função qX que leve 
uma n-forma em M para uma n-forma em U, 


Pa: FM) FU) 


de maneira que se w é uma n-forma em M, então dkw é uma n-forma em U. 


Logo temos que 
N w= WI 
or AM 


Esse é o procedimento no qual calculamos uma integral de uma n-forma w em 
um n-simplexo o” em M. 

O próximo passo é aglomerar vários simplexos de maneira que cada um seja 
homeomorfo com um subconjunto aberto de M e englobem toda essa variedade. 

Antes disso, devemos definir o que seriam as faces de um simplexo s”. As 
fronteiras de s” são dadas por 

ôs" = 53 dt; 
i 
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e então denominamos os (n-1)-simplexos t; de faces de s”. 
Seja (s”, U, 4) um simplexo em M pertencente à classe o”. Restringindo 
o mapa à às faces t; de s” e V; sendo uma vizinhança de t; no menor espaço 
plano que contém t;, então montamos um novo (n-1)-simplexo 7; em M dado 
por 
Ti= (t;, Vi, É) 


de maneira que podemos definir as faces 7; de um simplexo o” em M pela 
expressão 
Oo” = > +; 
F 


Iremos definir uma n-cadeia simplicial c em uma variedade M como sendo 
uma combinação linear de vários n-simplexos o? em M dada por 


c= Duo? 
i 


onde a; são constantes. À integral em uma cadeia simplicial é dada pela soma 


das suas partes: 
[o-Da [ 
c % Ci 


Um espaço topológico M é dito triangularizável se existir pelo menos uma 
cadeia simplicial c tal que seja possível construir um homeomorfismo entre M 
ec. 

É um resultado conhecido da topologia que todas as variedades diferenciáveis 
de quaisquer dimensões são triangularizáveis. Um exemplo de triangulação é 
ilustrado na figura 10, onde a superfície de um toro é coberta pela junção de 
vários 2-simplexos. Em cada simplexo o; podemos construir um conjunto aberto 
U em um plano ax + by + cz = d que o contém de maneira que possamos criar 
um mapa biunívoco entre U e uma parte da superfície do toro. E então todos 
os pontos da variedade são cobertos por um homeomorfismo com a cadeia de 
simplexos. 


Exemplo 4.3 (Integrando em um 2-simplexo no R$). Suponha que sº = [po,P1,9] 
seja um simplezxo no Rê e f(x,y) uma função diferenciável definida em s2. O 
simplexo padrão A? é dado por [(0,0), (1,0), (0,1)]. Queremos criar um mapa 
entre esses dois simplexos, simbolizado pela expressão 


togo + tagr + tago <> topo + tipi + topa 
Dados go = (0,0), q = (1,0) e go = (0,1) e denominando ty = u eto = 0, 
impomos t)=1—-u-—vw para que as coordenadas baricêntricas t; sejam únicas 


(teorema 4.2) e então 


(u,v) > (1 = u— v)po + upi + vpo 
(u,v) <> po + u(py — po) + v(p> — po) 
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Figura 10: Triangulação de um toro com raio externo d = 2 e raio interno r = 1 
usando 2-simplexos. 


Segue disso que as transformações de coordenadas são 


T(u,v) = po, + u(p1, — Po, ) + v(pa, — Po,) 
y(u, v) = Po, + u(pr, e Po,) E v(pa, ES Po,) 
z(u,v) = po, + u(pi, — po.) + v(pa, — Po.) 


2 


Agora suponha que o simplexo sº não seja paralelo a algum plano perpendicular 
ao plano xy, de maneira que possamos isolar z em função de x ey. Logo 


dx = (pi, — po,) du + (pa, — po, ) dv 
dy = (pi, — po,) du + (pa, — po, ) dv 


e então para uma 2-forma f(a,y) da dy temos 


ds (Fla, y) de dy) = ye f(a(u,v), y(u, v)) dudo 


onde 


Yf = ((p1, — po. )(pa, — Po,) — (pa, — po, )(p1, — Po,)) 
= ((pi — po) x (p2 — po)) |. 


é a componente z (não-nula) do vetor normal à superfície do simplexo s2. Segue 
disso que a integral 1 da 2-forma f(x,y) da dy em s2 é 


1= ff. ttoardy = ff fta(uc0) uu) dude 
sf fo fetal) dude 
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4.5 Teorema de Stokes generalizado 


Vimos na seção de revisão cinco teoremas fundamentais que relacionam uma 
integral em uma dimensão n < 3 com uma integral em uma dimensão menor ao 
custo de uma derivada. No teorema fundamental do cálculo temos 


b 
[ at= 10-10) 


no qual no membro esquerdo temos uma integral da derivada de f em um 
intervalo [a,b] e do lado direito temos f avaliado nas bordas desse intervalo, b 
e a respectivamente. O teorema do gradiente é similar, pois é a generalização 
do teorema fundamental para uma curva em R”. 

No teorema de Green, temos uma integral unidimensional se relacionando 
com uma integral dupla de derivadas parciais: 


“W(fQ op 
À. Pas) do + Olou) dy = ff (SS de dy 


A sua generalização para o R?, o teorema de Stokes, também é desse tipo, onde 
as derivadas se tornam no rotacional. 

E o teorema da divergência de Gauss temos uma integral de superfície se 
relacionando com uma integral tripla de um divergente: 


dh F-da = [) V-Fdr 
oV Vv 


Esses teoremas citados, que reduzem ou aumentam a dimensão de uma in- 
tegral inclusive fazem parte das demonstrações das integrações por partes. 

Usando tudo o que foi aprendido sobre formas diferenciais e a integração 
destas em variedades diferenciáveis, podemos mostrar que todos esses teoremas 
fundamentais são casos particulares de um único, que é o teorema de Stokes 
generalizado. 


Teorema 4.3 (Teorema de Stokes). Seja w uma p-forma definida em uma 
variedade orientável M e c uma (p + 1)-cadeia simplicial em M. Então 


[o= [au 
de c 


Demonstração. Sabendo que uma cadeia simplicial é uma combinação linear de 
vários simplexos, é suficiente provarmos que 


[= [au 


onde o é um (p + 1)-simplexo de c com uma representação em um simplexo 
padrão dada por 
(AP, q) 


130 


para algum mapa à. Segue disso que 


/ u=[ dou e [do=[ Ppri(do) 
do oAP+ & Apt 


Definindo u = djw e sabendo que q, ;(dw) = d(djw), devemos provar que 


K n= |, du 
CAP+I Ap+1 


A p-forma u pode ser representada em uma base dxf: 


H 
onde Ay = Ay(al,...,aP,xPtl) e H são subconjuntos de p índices distintos 
de (1,2,...,p + 13. Uma vez que a derivada exterior é linear, assim como a 


integral, nos é suficiente demonstrar apenas para o caso particular em que u é 
representado por uma única p-forma, 


u= Ada!...da? 
Logo 
OA 


Ogp+l 


del dy! ,..dr? = (-1P? det ada! 


oa Out 


Segue disso que 
OA 
= =[)P———— 1 p+1 
e au al 1) Out assa 


1-5. x a] 
= (e f ( Sede!) dg!...dz? 
A” AJo = 


p 
l du = (1º [ A (at, custa Biot) ast. cas 
AP+1 A? e 
(e f A(a!,...,x?,0)de!...da? (116) 
AP 


Por outro lado sabemos que A?*! = [0,e;,...,e,41], onde e; é o i-ésimo 
vetor da base canônica do R?P*!. Então 


0A?! = [e,,... epri]— [0,e5,... epa] +... + (—1PH0,61,...,e,] 
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Mas como yu = Adz!...dx” e pelo menos uma coordenada x*, i + p+1 é 
constante? em todas as faces exceto a primeira e a última, temos que as integrais 
em | nessas faces serão nulas. Então temos que 


A n= |, n+ (ape f u 
oAP+ [e1,...,ep+1] [0,...,ep] 


no qual identificamos [0,...,e,] = A?. Como «?*! = Q em todo esse simplexo 
padrão, teremos 


o di E ai O dg q A(x!,...,4?,0) de! ... da? 
0,...,ep Pp 


=—(—1)? A(x!,...,2?,0)da!...da? 


AP 
que é o segundo termo da equação (116). Já na integral em [ej,...,ep+1] pode- 
mos usar um mapa entre esse simplexo e o simplexo padrão A? = [0,e1,...,e,] 


usando uma generalização do exemplo 4.3: 


onde y são as coordenadas em A? e 
Pj; = Oij V9 A 0 
Po, = Oito) 


de maneira que o ponto po seja e,, 1 enquanto que os demais pontos p; são inal- 
terados?!. Podemos ver isso ao substituir os pontos pelos deltas de Kronecker: 


p 
2º = dip+1) + D, (Og — Oigp+1)) 
j=1 
Se 1 p+ 1, temos 
p 
ES poa (117) 


20No caso do tetraedro A?, temos que uma de suas faces é —[0, e», es], que é um triângulo 
retângulo no plano yz com o ângulo reto na origem e então x! = 0 em toda essa face. Portanto 
se houver uma forma, diferencial dx! na integração dessa face, teremos uma integral do tipo 
Te: que é zero. 

21No caso do simplexo [e1,e>,es3], que é um triângulo cujos vértices são (1,0,0), (0,1,0) 
e (0,0,1) respectivamente, temos que po = (0,0,1), enquanto que pi = (1,0,0) e p2 = 
(0,1,0). O simplexo padrão A? terá os mesmos vértices do outro simplexo porém com (0,0,1) 
substituído por (0,0,0). 


132 


Mas sei=p+1, 


Pp 
=1+5 y'ô;-D,y 
j=1 


p 


=1-S,y 


j= 


uma vez que j *i = p + 1 em todo o primeiro somatório. Usando o resul- 
tado (117), ficamos com 


p 
aPti =1— a 
j=1 
p 
=1-— x! 


Além disso, temos dx” = dy' para todo i Lp «1. 

Porém nesse mapa trocamos o simplexo [e1,...,ep,ep+1] por [e1,...,e,,0] 
(já que usamos po como sendo ep,1). O simplexo padrão A? é obtido ao mudar a 
orientação dos vértices permutando o vértice O da última posição até a primeira, 
o que deve carregar um fator (—1)?. Então temos que 


n= |, A qasaçab Ped corda? 
[e1,...,ep41] [e1,...,ep,ep+1] 
Pp “ 
-[ A (et, custa= Sie) as! cao 
[e1,...,ep,0] mr 
Pp is 
(1 [ A (at, custa Sia) ast. cao 
[0,e1,...,ep] o 
p , 
(1 / A (at, custa Sie) at. car 
Ar ; 


que é o primeiro termo da equação (116). Logo 


/ n= |, du 
CAP+I AP+1 


o que completa a demonstração. E 


Corolário (Teorema Fundamental do Cálculo). Embora tenhamos usado o Te- 
orema Fundamental do Cálculo para demonstrar o teorema generalizado de Sto- 
kes, podemos recuperá-lo ao fazer o = [a,b| um intervalo fechado em R. Então 
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do = [b] — [a] e para uma função integrável w = f temos 


he E a f = (6) — fa) 


fue fu 


Então o teorema de Stokes afirma que 


enquanto que 


b 
! af = f(b) — f(a) 


Corolário (Teorema do gradiente). Dado uma função escalar diferenciável 
V(al,...,2”), podemos fazer w = V, de maneira que 


HOV o, 
dw = > 7 dz 
i=1 


Sendo c a cadeia 1-simplicial dada pela junção de vários segmentos de retas no 
R” que une dois pontos P e Q, a curva resultante entre esses dois pontos é 
suave por partes. Denotando essa cadeia como 


c=[PRi]+ Ri, Ro)+...+ Re, Re] + [Re Q] 


onde Ri; denota os pontos comuns entre dois segmentos de reta, temos que a 
fronteira de c se torna uma soma telescópica, no que 


de = —[P] + [0] 
Logo 


Portanto o teorema generalizado de Stokes afirma que 
[vv “de= V(Q)— V(P) 


que é o teorema do gradiente. 


134 


Corolário (Teorema de Green). Sejam P(x,y) e Q(x,y) duas funções dife- 
renciáveis em x ey pelo menos até primeira ordem e definidas em um subcon- 
junto aberto X do R2. Seja D uma região fechada em X delimitada por uma 
curva suave por partes OB. Fazendo 


w=Pdr+Qdy 


temos 
dw = SE india do 
Oy dx 
0Q OP 
Então 


/ = P(x,y) de + Q(x,y) dy 
ox 5)» 


[de= | (5 =) dx dy 
>) sz Ux oy 


Logo o teorema generalizado de Stokes nos garante que 


“MtoQ op 
O Pesar + Aaapdy = ff (E D)aray 


que é o teorema de Green. 


enquanto que 


Corolário (Teorema de Gauss da divergência). Seja F um campo vetorial de- 
finido em um subconjunto aberto X de R$. Seja V uma variedade diferenciável 


volumétrica em X com coordenadas locais ortonormais al = u,z2 = vex) = w. 
O campo F pode ser dado por 
F = Fã + Fé + Fes (118) 


onde fe;) forma uma base ortonormal do Rê, com e; sendo a derivada normali- 
zada do vetor posição em função de x". Esse campo pode ser escrito de maneira 
equivalente pela 2-forma 


w = F,0203 + F50301 + F301093 (119) 


onde F, = Fi(u,v,w) e o elemento o10903 forma uma base ortonormal das 3- 
formas em R$. As formas diferenciais o; são as mesmas calculadas na seção 3.7 
sobre o operador nabla em coordenadas ortonormais: 


o, = Adu 
02 = udu 
03 = vdw 
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Segue da equação (119) que 


w = (uvF,) dv dw + (AvFs) dw du + (AuFs) du do 


RENO dv) GOMES) QUE 
do = ( (um DU a 0. faia 2) dudu day 
ou Ou 
ou 


ou Ou ow 


Também vimos na seção 3.7 que oj0203 é o elemento de volume dr e que 


v.F 1 (Eee e AvFs) m ni) 


E Av ou Fox) dw 


is 1 (Sua 4 (AVE) 4 aa 010203 


Logo 
dvu=V-Fdr 


Por outro lado, vimos na seção 3.8 que o elemento de área infinitesimal da 
quando a coordenada xº é constante é dado por 


da = 0102€3 
Em geral, o elemento de área da é dado por 
da = os03& + 030162 + 010963 
e portanto 


F.da= (Fé + bBbé& + F363) : (co03& + 030163 + 010263) 
= F,0903 + F503014 + F301092 


=W 


Logo 


K a = dh F.da 
oV ov 
fo = | V.Fdr 
V V 


Segue do teorema generalizado de Stokes que 


dh F-da= [) V-Fdr 
ov v 


que é o teorema de Gauss da divergência. 
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Corolário (Teorema de Stokes do rotacional). Seja F um campo vetorial pelo 
menos uma vez diferenciável definida em um subconjunto aberto X de R$ dado 
pela equação (118). Seja D uma variedade diferenciável em X delimitada por 
uma curva suave por partes 05. Esse campo pode ser escrito de maneira equi- 
valente usando a 1-forma 


wW = Fo1 + Fã0s5 + F303 


nas coordenadas curvilíneas ortonormais x* cuja base é (1, é2,€3!. Conforme 
vimos na seção 3.7, os elementos de comprimento são 01, 03 € 03, No que 


dr = cj€& + os€2 + 0363 
e então 


F.dr = (FR é& + bBé + F3€63) . (o1& + 0962 + 0363) 
= Fo, + Fo509 + F303 


=Ww 


Além disso, dado 


o1 = Adu 
02 = udu 
03 = vdw 
temos 
WU = (AF) du + (uF») dv + (vF3) dw 
e então 
AF F, 7) 7) 
dw = dida) du + Bu) dw | du +... + Puks) du + SF) dv ) dw 
Ov Ow Ou Ov 
Já fizemos esse cálculo na seção 3.7 e resulta em 
1/0 9) 1/0 9) 
ss uv fm (vFs) ow (1F3) )ozos ii Av (E E) du (E) Josor+ 


1/0 Õ 
Au (= (utto) Ex (AF1)) 0102 


Por outro lado, sabemos também que 


1/0 0 E 1/0 0 a 
EE uv E (bs) dw (Fa) )6s Es Av (au Uai ou (2F3) )62+ 
1,0 Õ q 
Au E (ufa) Ov (NF) ) 6s 
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e que 
da = 090361 + 0301€2 + 010963 


Logo 
VxF-.da=dw 


fao= | veda 

) õ 
[u-d F.dr 
oz FD) 


Portanto o teorema generalizado de Stokes afirma que 


Jovxr-a-d F.dr 
õ os 


que é o teorema de Stokes no Rº. 


Finalmente temos que 


e também 


Exemplo 4.4 (Integração por partes generalizada). Dado o teorema generali- 


zado de Green 
/ w= [ao 
de c 


se fizermos wu = u A À com | sendo uma p-forma, obteremos 


fiuna= [annx+tapus ar 


fusna= [dunas funds 
fauna= | unaj-(1P fumar (120) 


A expressão (120) é a forma generalizada das seguintes integrações por partes: 


e então 


b b 
= F Hag'(x) de 


2. fa conar= db tada ff fev -ajar 
o) Lico xada-d tacar+ fl IA (Vol da 
à hf Beco x ajar = db (A x B)-da+ JJ A-(V x B)dr 


A primeira é obtida ao fazer u = f,A=gec=l[a,b]. A segunda advém da 
escolha u = f, À = Ay0903 + As0301 + Agoios e c um volume. Já a terceira 
c é uma superfície aberta, u = Ajoy + Asos + Agog e A = f. Por último, a 
quarta integração por partes pode ser obtida usando o mesmo u do terceiro item 
além de fazer À = Bi01 + B509 + B303 ec um volume. 


b 
ê / F(a)g(a) de = f(x)g(2) 
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4.6 Formas exatas e fechadas 


Definição 4.7 (Formas fechadas). Uma forma diferencial w em uma variedade 
diferencial M é dita fechada em um subconjunto XE M sedvu=0 em X. 


Definição 4.8 (Formas exatas). Uma forma diferencial w em uma variedade 
diferencial M é dita ser exata em um subconjunto X E M se existir alguma 
forma a em X tal que da = w. 


Das definições segue imediatamente que toda forma exata é fechada, pois 
d(da) = 0. Porém nem toda forma fechada é exata, apesar de existirem 
condições suficientes para isso ser verdade conforme vimos na inversa do lema 
de Poincaré. 


Teorema 4.4 (Integrais de formas fechadas). Se w é uma p-forma fechada em 
Fº(M), com M sendo uma variedade diferenciável, então dada duas n-cadeias 
simpliciais quaisquer cy eos em M tais que 04 — 02 seja a fronteira de uma 
n + 1-cadeia simplicial em M temos 


oq o2 


Demonstração. Seja D tal que 0, — 05 = 05. Então 


[o [u-[ u= [ u= | du 
o (0) 01—092 os DD) 


Mas como w é uma forma fechada em M, temos dw = 0 em X e então 


[o-[ w=0 
o1 [0 8») 


Exemplo 4.5 (Integrais de caminho em campos conservativos). Seja f uma 
função escalar diferenciável pelo menos até primeira ordem definida em um 
conjunto X € R” assim como dois pontos p eg em X. Sejam c, e c» dois 
caminhos suaves por partes que começam em p e terminam em q. Fazendo 
w = df, temos que w é uma forma diferencial exata e portanto é também uma 
forma fechada. Uma vez que cy — c» é o caminho que vai do ponto p até ele 
próprio passando pelo caminho cy e depois —ca, segue que cy — cs forma um 
caminho fechado e consequentemente é a borda de alguma superfície 5 em X. 
Então pelo teorema (4.4) temos que 


fo hio 
ci Ca 


Mas também vimos que df = Vf-dr e portanto 


/ vi dr= [ Vf-dr 
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Isso significa que para um campo vetorial conservativo F = Vf a integral de 


um ponto p até um ponto q independe do caminho. Esse resultado é equivalente 
ao teorema do gradiente, 


! Vj$-dr= f(a)- f(p) 


onde a integral de caminho só depende dos pontos inicial e final. 
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5 Aplicações de formas diferenciais na Física 


5.1 Eletromagnetismo 


Para obtermos as equações de Maxwell na notação de formas diferenciais, iremos 
começar com a mais fácil, V.B = 0, onde B é o campo magnético. Definindo 
a forma 


B = Adydzdt+ Cdzdxdt+ Ddzdydt+Fdardydz (121) 


temos que 


Fr oA oO OD 
dB & 9 oC ' Z 
dt dx Oy Oz 
Fazendo F=0, 4-B,, C-B,e D= B,.,temos 


) dx dy dz dt 


dB =(V.B)dxdydz dt 


Segue disso que a equação de Maxwell da não-existência de monopolos magnéticos 
é dado por 
dB =0 


onde 
B = B, dydzdt + B, dz dx dt + B, dx dy dt 


Agora iremos verificar o que acontece se eliminarmos dt na equação (121). Seja 
À definida por 
Ao = Br dydz + B, dzda + B, dx dy 


de maneira que B = Ag dt. Então 


oB B, 


X B a) 
do = (V - B) dedydz + TE dydedt + de da dt + Gg dx dy dt 
9 
0 By B OB. 
ddr Cl dndra Carai 
ot ot ot 
: oB ; sr 
Para obtermos a lei de Faraday V x E + Em 0, onde E é o campo elétrico, 


devemos encontrar uma 2-forma À, tal que 


mel Cen Cria o dd 
Oy õz z Fopi dx Oy 
(122) 


de maneira que dAg + dA; = 0, uma vez que essa equação de Maxwell nos dá 


7) e ú 
1) 0E, | 0B, 0 
Õz Fohi ot 

dE, 0E, 0B. 0 


Naturalmente À, só pode ter as componentes do campo elétrico. Como não 
aparece nenhuma derivada temporal em dAy, segue que em todas as 2-formas 
que compõe A, possuem um dt. Além disso, não aparecem derivadas parciais 
de E, em relação a x e idem para as outras coordenadas, o que significa que 
deve ter uma forma dx onde há E, e assim por diante. Então intuitivamente 
chegamos na expressão 


A = E,drdt + E, dydt+ E, dz dt 


que cuja derivada exterior de fato reproduz (122). Portanto, denotando a = 
Ao + As, temos que 
da = 0 


com 
a = (E, dar + E, dy + E. dz) dt + Br dydz + B, dzdr + B, de dy 


nos dá as equações homogêneas de Maxwell no sistema internacional (ST): 


vV.B=0 

oB 
VxE+>— = 
x + 0 


No sistema gaussiano, devemos apenas substituir a coordenada t por ct em todos 
os passos além de considerar a forma diferencial cdt = d(ct) como um elemento 
da base das 1-formas. 

Ainda no sistema gaussiano (que nos será útil no contexto da relatividade 

restrita), as duas outras equações de Maxwell são as não-homogêneas: 
V.D=47p; 

1 oD 

VxH=-—[47],+— 

c ot 


onde py e Jy são respectivamente a densidade de cargas livres e a densidade 
de corrente de cargas livres, além de c ser a velocidade da luz no vácuo, D o 
deslocamento elétrico e H o campo magnético auxiliar. 

Suponha que 


D = Dii+ Dj + D.k 
H-= Hi+Hj+ H.k 
Jp =Jri+Jçj+ Jr 


Então a lei de Ampére-Maxwell nos dá 


0H, 0H, 1 OD o 
Em ua (am, )-o 
0H. “0H, d DV 
E Ea (am + E) Ee 
0H, 0H, 1 9D, 

Es Re | ao = 
Foto Oy c ( tia ot ) É 


Como agora as derivadas parciais em ct aparecem com sinal negativo, iremos 
definir uma 2-forma 1) dada por 


y=-—D, dydz — D, dzdx — D, de dy 
de modo que 


10D, 10D 
Dies D)dadudass É dy dz det — = RR 
c O c Ot 


PE” da dy det 


Queremos encontrar a expressão 


0H. 0H, 1 0D, 
(47py — V- D) dr dy dz + | Ro E E (amo, + mM )| dy dz det+ 


0H, 0H, 1 0D, 
0H, 0H, 1 oD. 
Get CO ldedyde= 
[5 dE ( tJf, + En )| adydet=0 (123) 


Assim como nas equações homogêneas de Maxwell, devemos somar uma 2-forma 
do tipo 
H, da det + H, dydet + H, dz det 


à 2-forma 1) para que apareça V x H. Então se definirmos uma 2-forma 5 dada 
por 


8 = —Ds; dydz — D, dzdr — D, dr dy + H, da det + H, dydet + H, dz det 


então 
0H. 9H 10D, 
dB =(-V- D)drdydz + |É Ciy “Coe | qydzdet+ 
Oy z c Ot 
0H, 0H, 10D, 
EE o Ed fas da det+ 


[Ou 0H, 10D, 


E Em E Õ | dz dy det 


Para obtermos (123), devemos somar uma 3-forma 41% dada por 
1 1 1 
47% = 47py da dy dz — Ele dydz det — — Jg, dz da det — — Jg, da dy det 
c c 
ou equivalentemente 
y = py dardydz — Jy, dydzdt — Jy dzdar dt — Jy, dx dy dt (124) 


de maneira que 
dB + 477 =0 
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resulta em (123). 
Portanto dada as formas diferenciais 


a = (E, dr + E, dy + E, dz) det + B, dydz + B, dzda + B, dr dy 
8 = —Dy; dydz — D, dzdx — D, dr dy + H, de det + H, dydet + H, dz det 
y=pydxdydz-— Jy, dydzdt — Jp dzda dt — J, da dy dt 


as equações 


do = 0 (125) 
d8+47y=0 (126) 


resumem as quatro equações de Maxwell. 
Aplicando a derivada exterior na equação (126), obtemos 


dy =0 


que após recolocar dct em (124) essa equação se torna 


10 10J, 10J, 10J, 
Sparad dct dx dy dz— — É da dy dz det— — E dy dz da det- =“ da dz dx dy dct = O 
c Ot c Or c Oy c Oz 
e 1 dp 0TE- Oi 8] 
CPf fa O ty fe 
t = 
c ( ot dz Oy Õz ) au Apto 

e portanto 

ps 

Ed, Ei 

En +V f 


que é a equação de continuidade ou conservação local de cargas livres. 

Uma vez que da = 0, então pela inversa do lema de Poincaré em toda região 
estrelada do espaço-tempo existe alguma 1-forma À tal que dA = a em toda 
essa região. Denotando 


A= Adr + A, dy + A, dz + Ag det 


A-=Ai+Aj+Ak 


temos que os termos de dA que não envolvem derivadas em dct resultarão em 
2-duas formas do tipo dx dy, dz dx e dy dz, e então concluiremos que 


VxA-B 


Já nos demais termos, obteremos 


Roo da det + aaa dy det + eso dz dct— 
fop Oy Õz 
104, 104, 104, 
E da det Rr dy dct di dz dct 
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ou 


OA, 104, 04, 104, 04, 104, 


que pela equação dA = a essa última expressão é igual a 
E, da det + E, dy det + E, dz det 


e portanto chegamos nas equações 


OA, 104, 
SRD Es e, 
Ox c dt 
0 Ao nd OA, 
dy cdt 
0 A, il OA, =. 
oz edi * 
Denotando Ay = —V como sendo o potencial elétrico, essas expressões se redu- 
zem a E 
—VV-=>>— =E 
c Ot 


A forma diferencial a dada por 
a = (E, de + E, dy + E. dz) det + B, dy dz + B, dzde + B, de dy 


pode ser escrita na base do espaço co-vetorial das 2-formas no espaço tempo da 


seguinte forma: 

4 u 

a= > > Fudut do” 

u=1v=1 
onde du! = dz, du? = dy, du? = dz e du! = det. Então os coeficientes HF, 
tomam a forma da matriz antissimétrica 
O —E, —-E, —E, 
Es ) B. -—B, 
EL, —B, 0 Be 
EL. B, —Bs 0 


que é a representação matricial do tensor eletromagnético em sua forma covari- 
ante. 


5.1.1 Equações no vácuo 


No vácuo temos pp; =0,D - E, H = Be J= 0. Então as formas diferenciais 
a, 8 ey ficam 


a = (E, de + E, dy + E. dz) det + B, dy dz + B, dz de + B, de dy 
8 = -E, dydz — E, dzdz — E, dr dy + B, dr det + B, dy det + B, dz det 
Y=0 
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Escrevendo 5 pela expressão 


temos que 


Eu Bo Ep D E; 
FE 


que é a representação matricial do tensor dual de F,,. 
Definindo a norma das formas diferenciais como 


(da, dz) = <dy, dy) = <dz,dz) = 1 


(det, det) = —1 


e considerando que da dy dz d(ct) tem orientação positiva, então temos que essa 
base das 4-formas tem outras representações com orientação positiva dadas por 


da dy dz d(ct) = da d(ct) dy dz = dy d(ct) dz de = dz d(ct) da dy 
= de dz d(ct) dy 


dentre outras. Relembrando a equação (28) da seção 2.4 sobre o operador es- 
trela, temos que se dx” À dxÉ é uma base das n-formas de um espaço cotangente 
de dimensão n, então 

«(dr) = (dr*, dal doÉ 


Isso significa que 
«(da det) = (dy dz, dy dz) dy dz 


Mas 


(ay dz, dy dz) = der | Rd o | E der E | po 


e então 
«(da det) = dy dz 


Analogamente, temos 


«(dy det) = dz da 
«(dz det) = dx dy 


Por outro lado 


«(da dy) = (dz det, dz det) dy dz 
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a qdades (dede |. 1 0 |. 
(dz det, dz det) = det | tdetides: Ade desse] det asda fe —1 


ou seja, 
«(da dy) = — dz det 


Similarmente, 


«(dy dz) = — de det 
«(dz da) = — dy det 


Com isso, temos que 


sa = E; dydz + E, dzda + E, dedy — By da det — B, dy det — B, dz det 


=-8 
Portanto as equações de Maxwell no vácuo ficam 
da = 0 
d(—- a) = 0 
ou 
da = 
d(za) = 0 
Ágora note que 
0E 0E 0E. 0E, 
d(—8) = 2d “ det | dyd pa “ det | dzd 
(— 8) (5 a ct) dude + (SS a ct) dedo 
E. LE, By Bs 
Ê dz + á dct | dx dy di dy + Ê dz |) dz det— 
õz Oct Oy Õz 


2by da + Ely dz | dy dct qe da + Es dy | dz dct 
dx Õz dx Oy 


a(=B) = 0Es 4 0oB, 0B, ddr des E, a oB. OB, drdos 
Oct 


Oct z Oy dx Õz 
0E,  0B, 0By dE, 0E, 0E 
A = da dy det + * | dedyd 
(a Oy 2) do dydet + (Si 4 E) dedyas 
Então d(xa) = O nos leva às equações 
e VxB=0 
c Ot 
V.E-0 (127) 


Mas também 


a d(xa) = (5 VxB) d+ (CS -VxB) dy+ 
x y 


(25 V=B) de —(V- E) det 


e portanto 


d(* d(xa)) = [v = (2% VxB)) dyde+[V =» (5 -VxB)] dz dx+ 
y 
uva B)| dx dy— (5 
c Ot x hi 
1 
c 


7 7 
RO (Roo, 42 170 MAE DS 
E +a( 


c Ot 
Õ o (10E 


Usando a equação (127) e sabendo que 


Vx(VxB)=V(V-B)- V?B=-VêB 


VxE= sen 
ot 
7) oB 3 
——(V x B) V x — Vx(VxE)=-VÍE 
oct t 
temos 
do 2B 1 02B 5 
d(x d(z0)) = [-5 22 Ea B|, dz dr+ 
102B IPÊ a 1 da E 
0? 
— nes a Eee dz det 
c2 Ot? x 
Porém como d(xa) = 0 e o operador estrela é linear, segue que d(x d(xa)) 


também é zero. Então, pela independência linear das duas-formas temos 


1 02B 


2 
-352 + VB-0 
12E 

c2 Ot? io 


que são as equações de onda para E e B. 
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5.1.2 Teorema de Poynting 


Agora sejam as formas diferenciais de um espaço vetorial tridimensional dadas 
por 

w; = Erdr + E, dy+ E, dz 

wa = By dydz + B, dz dar + B, dr dy 


Denotando d” como a derivada exterior desconsiderando termos do tipo dct, 
temos 


dw =(V x E).dydz+(V x E),dzde+(V x E), dedy 
d'wo, = (V - B) de dy dz 


Também denotando a derivada temporal de uma forma diferencial pela ex- 
pressão 


Õ oB, oB oB 
= wo = E dyd = dzd 2 de d 
a? wa ot y dz + En Rb x dy 
segue que as equações homogêneas de Maxwell se tornam 
d'w = ——wa 
c 
d'wo =0 


Para obtermos as outras duas equações, que são as não-homogêneas 


V.D=4rp; 


1 oD 


precisaremos calcular o divergente de D e o rotacional de H, no que nos motiva 
a definir 


wy = H;dr+ H,dy-+ H, dz 
wa = Dy dydz + D, dzdax + D, dx dy 


de modo que 


dws=(V x H)dydz+(V x H),dzdr+(V x H), de dy 
d'ws =(V- D) drdydz 


Como a fórmula para gerar o rotacional resulta em duas-formas, também iremos 
definir 
ws = Jy dydz + Jy dz da + J, de dy 


Então as equações não-homogêneas de Maxwell ficam 
T 1. 
d'ws = —ws + >w4 
c c 


d'ws = *(4mpf) 
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onde as normas (dx, dx>, (dy, dy) e (dz, dz> são todas iguais a +1. 
O vetor de Poynting na notação vetorial é dado por 


c 
=—ExH 
Ss E x 


Vimos que o produto vetorial pode ser representado pelo produto exterior das 
1-formas. Nesse caso temos 


Sx dydz + Sy dz da + S. drdy = qt Aus 
T 
Ágora 
d'(wy A w3) = d'w1 A ws + (=1) Hc A d'ws) 


1. (= 1 1) 
= ——Wo AW3s — Wi A | —ws + —wa4 
c c 


c 
; T é 
=——W9o AWs — — WI A W5 — —Wj A W4 
c c c 
Mas por outro lado 
4 
d'(w1 A wa) = * (Ev . s) 
É 
e então 
4 1 4 1 
«(EVS) inn + Cunnar + Cu adg=o (128) 
c c c c 


Também vimos anteriormente que o produto escalar pode ser representado pelo 
produto exterior de uma 1-forma com uma 2-forma. Então a equação (128) nos 
fornece 

41V .S+B.H+47E.J+E-D=0 (129) 


que é o teorema de Poynting. 
Em um material em repouso com comportamento dielétrico e de magne- 
tização linear temos 
D =eE 
B = uH 
onde e é a permissividade elétrica e yu a permeabilidade magnética do material, 


ambos constantes no tempo. Substituindo essas expressões em (129) e dividindo 
por 47 iremos obter 


v.s+2H.H+E.J+ ÍE.E-=0 
47 47 


Mas 
- lo 10B? 
E-.E=--—(E.E)=-— 
> ul 2 
. 2 
gm 1H 
2 Ot 
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e então 


o(1 
E.J+—([>(uHº +eE?) |) = 
V:S+ + a (ara +e ) 0 
Denotando 


1 
Eae E H? E2 
a (uH" + eBº) 


como sendo a densidade de energia nos campos eletromagnéticos, concluímos 


que 


ou 
a AR = e 
A V.S+E.J 


5.2 Termodinâmica 
5.2.1 Relações de Maxwell 


Os potenciais termodinâmicos molares são dados por 


u(s,v) - fTas-pas 
HMTv=u—Ts 
h(s,p) = u+ po 
Tv) =u+pv—Ts 


onde f é a energia livre de Helmholtz, h a entalpia e g a energia livre de Gibbs. 
A derivada exterior da energia interna molar é dada por 


du=Tds-—pdv 
enquanto que para os demais temos 
df=du-sdT-—-Tds 
df=-sdT — pdv 


dh = du +vdp+pduv 
dh=Tds+vdp 


dg=du+vdp+pdv— sd” —- Tds 
dg=vdp— sd” 


Pelo lema de Poincaré temos d(du) = 0, d(df) = 0, d(dh) = 0, e d(dg) = 0. 
Mas por outro lado, 


(97 oT Op op 
d(du) = (Sel, ds + So], do) ds (2), ds + 52], du) do 
= (5 4. 0P )avas 
Ou ls Oslv 
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e consequentemente 
(130) 


onde o traço vertical indica a variável que é constante durante a derivação 
parcial. A equação (130) é uma das quatro relações de Maxwell. As demais se 
seguem da mesma maneira: 


d(df) = — (5 AT + : ndo) dT — (5 AT + o pão) dv 
— (ô0s op 
= (= = oT ) gem 
oT Õ 
d(dh) = (|. =) ras 
Ou Os 
e então 
ôs| Op 
Fole” Ph, Fo) 
oT Ou 
oi = el, (132) 
Ou Os 
ode 8 (133) 


Pela distribuição da derivada exterior, também temos 


dT ds = dp dv 


5.2.2 Propriedades dos materiais 
Dado um material, na termodinâmica definimos as grandezas 


Os Tê 1ldv “Id 


Gps DT = 0T lp “TE vôplr “Cuório 


onde cy, Cp, KT € Q são respectivamente a capacitância térmica molar a volume 
constante, à pressão constante, compressibilidade à temperatura constante e o 
coeficiente de expansão térmica. 

A forma ds em termos da temperatura e do volume molar é dada por 


Os Os 


ds v OvIT 
-2a+|d 
“T qo” 


Os 
Ou 
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Mas escrevendo dv em termos da temperatura e da pressão temos 


Ou Ou 
do = Sa) T+ sobr? 
e então ã E - E 
Cy s v 5 v 
di (2+5 el) ar+ (Gl, er) do 


Por outro lado, escrevendo ds diretamente em termos da temperatura e da 
pressão temos 


Os Os 
SE dra bd 

Tl, * gplr 
Os 


odpIr 


ds 


Cp 
= E dT 
T + 


dp 


Como as formas dT e dp são independentes, comparando as expressões chegamos 
nas equações 


Cp CG dd Ou 
= 134 
Tr T Ovlr dTlp ess 
ú ôs| O O 
da Re (135) 
opir - Odvlr Oplr 
Na equação (134) identificamos 
Ov 
ES SA 136 
Orla E too) 
Já na equação (135), temos que 
o E oo T 
Oulr o 
ôp lr 
Usando uma das relações de Maxwell (equação 133), obtemos 
: du 10 
s oT o voT p a 
= = 137 
Gu ou — 1d RT sm 
op T vop T 
Substituindo (136) e (137) em (134), concluiremos que 
Sopas Rc 
FRRNRETA KT 
ou equivalentemente 
2 
cp = Cy + Tu— (138) 
RT 


A equação (138) traz uma relação entre cp, Cy, Q € KT, O que indica que pelo 

menos um desses parâmetros não é independente dos demais. Usando a segunda 

lei da termodinâmica é possível mostrar que cp, Cy € KT São necessariamente 

positivos. Portanto da relação (138) também podemos concluir que cp > cy. 
De maneira análoga, começamos com 


Ov Fox) 
dv = dp+ —| ds 
op s P Os p 
Ov 
=-—uv-nsdp+ —| ds 
Os p 
onde «5 é a compressibilidade à entropia constante, dada por 
10v 
RR 
é v Ops 


Agora fazemos uma mudança para as coordenadas termodinâmicas para p e T, 
obteremos 


do =v-msdp+ 2] (5 aa e ar) 
p T 


Os op 0T Ip 
Ou Os Ou Os 
= (ut lo br) do+ (5: pl, ) ar 
Mas também 
Ou Ou 
dv = o dp + Rs 


=-v-nrdp+ did dT 
p 


Comparando as duas últimas equações, temos 


Ou Os 
Emis E ape das é — 139 
RE AR p Oplr fi 
ú du, do O 
==>) 5 (140) 
oT p Os p oT p 
De (140) temos 
R du 10 
OU To mp To my! 
Oslp ds Ts, Cp 
oT ã oT p 
Usando (133) em (139), obteremos 
a Ou 
: = DR E ; 
V-kr V-K ço |, 
a 
vinrp=—v-ks— To— - au 
Cp 
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e então 
2 


fer ms Tv — (141) 
Cp 
Isso significa que, pela segunda lei da termodinâmica, temos K7 > ks. 
Por fim, de (141) temos 
(kr — Ko)cp = Tua? 
Por outro lado, a equação (138) nos leva a 


(Cp — Co)kr = Tva? 


e consequentemente 
(Kar — Ks)Cp = (Cp — Co) 


—KsCp = —CykT 

ou então je a 
> (142) 
KT Cp 


Ás grandezas em que temos mais facilidade de medir são a expansão térmica 
a, a capacidade térmica molar à pressão constante c, e a compressibilidade à 
temperatura constante «7. Usando 


a? 


Cy pd 
T 


obtemos a capacidade térmica molar com volume constante. E então, utilizando 
a equação (142) obtemos a compressibilidade à entropia constante «,. 
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